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FOREWORD 


I This  Summary  Report  by  Arthur  L.  Smirl,  Department  of  Physics,  North 

Texas  State  University,  Denton,  Texas,  covers  research  progress  for  the 
! period  September  1,  1976  to  August  31,  1977,  under  Office  of  Naval  Research 

a Contract  NR  318-048. 

The  use  of  trade  names  in  this  report  does  not  constitute  an  official 
I endorsement  or  approval  of  the  use  of  such  hardware  or  software.  This 

report  may  not  be  cited  for  the  purposes  of  advertisement,  and  its  publi- 
| cation  does  not  constitute  Navy  approval  of  findings  or  conclusions. 
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SUMMARY 


In  this  work,  we  report  on  our  progress  in  attempting  to  use  the  high 
electric  fields  and  the  ultrashort  pulses  that  can  be  derived  from  Nd:glass 
and  CC>2  lasers  to  investigate  hot  electron  effects  in  semiconductors.  In- 
tense ultrashort  optical  pulses  having  durations  of  a few  picoseconds  and 

3 

peak  powers  of  10  watts  can  be  readily  generated  by  mode-locking  a Nd: glass 
laser.  The  extremely  high  power  and  short  duration  of  pulses  from  these 
lasers  make  possible  the  study  of  saturable  optical  transmission  properties 
and  hot  electron  dynamics  of  semi conduc tors  on  a picosecond  time  scale.  In 
addition  to  this  study  of  carrier  dynamics,  we  Me  also  attempting  to  inves- 
tigate stress  and  electric  field  biasing  of  the  absorption  properties  of 
semiconductors  under  intense,  short-pulse  C02  laser  radiation. 

Hence,  we  report  on  the  development  (with  coworkers  at  the  University 
of  Arizona)  of  an  initial  model  (see  Appendix  A)  that  describes  the  gener- 
ation and  the  temporal  evolution  of  hot  electron-hole  plasmas  produced  in 
the  semiconductor  germanium  by  the  absorption  of  intense,  ultrashort  optical 
(1  “ 1.06  pm)  pulses.  *niis  model  is  shown  to  be  in  agreement  with  the 
experimentally-observed  enhanced  transmission  of  single,  energetic  pulses 
through  thin  germanium  samples  as  a function  of  incident  pulse  energy  for 
two  sample  temperatures,  105  K and  300  K.  It  has  also  successfully  des- 
cribed the  transmission  of  a weak  probe  pulse  through  the  germanium  as  a 
function  of  time  delay  after  an  energetic  excitation  pulse  for  different 
excitation  energies  and  for  sample  temperatures  of  105  K and  300  K.  Subse- 
quent to  the  development  of  this  model,  van  Driel  and  coworkers  conducted 
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further  tr.in.smi  ss  ion  studios  (soo  Appendix  B),  in  which  tho  energy  band  gap 
'.I  i h<-  germanium  was  tuned  by  hydrostatic  pressure,  that  corroborated  this 
proposed  model. 

Recently,  however,  the  optical  pulsewidth  dependence  of  thin  qermanium 
samples  was  experimentally  measured,  as  suggested  in  our  earlier  work  (see 
Appendix  C) . We  found  these  measurements  to  be  in  substantial  disagreement 
with  the  predictions  of  our  model.  The  question  then  arose  as  to  whether 
this  disagreement  with  the  proposed  model  was  due  to  assumptions  made  in 
deriving  the  model,  experimental  uncertainties  in  physical  constants  used 
in  the  calculations,  or  to  relevant  physical  processes  omitted  from  the 
model.  A numerical  investigation  (Appendix  D)  of  the  effect  of  the  electron- 
phonon  coupling  constants  on  the  model  revealed  that,  given  the  uncertainties 
in  these  and  other  constants,  the  pulsewidth  experiment  was  not  a definitive 
test  for  the  model.  In  fact,  given  the  uncertainties  in  certain  of  the 
physical  constants  contained  in  the  calculations,  many  of  the  questions  con- 
cerning the  physics  of  ultrafast  phenomena  in  germanium  can  only  be  answered 
by  further  experimental  investigation. 
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ULTRAFAST  TRANSIENT  OPTICAL  RESPONSE  OF  ELECTRON-HOLE 


PLASMAS  IN  SEMICONDUCTORS 

Theory 

During  the  past  year,  the  theoretical  activity  in  this  project  has 
been  a detailed  study  of  the  theory  for  the  transmission  of  ultrashort 
optical  pulses  through  germanium.  At  the  beginning  of  the  year,  an  initial 
model  had  been  developed  by  A.  Elci,  A.  L.  Smirl,  J.  C.  Matter  and  M.  0. 
Scully.  Smirl  and  Matter  had  also  completed  single  pulse  and  excitation- 
probe  pulse  experiments  with  germanium.  In  addition,  an  attempt  had  been 
made  to  generate  numerical  solutions  to  the  complicated  coupled  set  of 
integro-dif ferential  equations  in  the  theory. 

However,  there  were  several  problems  in  the  original  computer  codes. 

The  computer  program  was  lengthy  and  expensive  to  execute,  and  the  solution 
to  the  differential  equations  was  crude.  Several  modifications  were 
needed  to  improve  the  accuracy  and  execution  time. 

There  were  four  areas  in  which  the  numerical  techniques  were  improved. 
First,  the  execution  time  was  considerably  lowered  by  using  Gauss-Laguerre 
quadrature  formulas  to  evaluate  the  integrals.  Next,  the  accuracy  of  the 
solution  to  the  differential  equations  was  improved  by  using  a fourth  order 
Rungo-Kutta  solution.  The  third  major  area  of  improvement  was  the  intiali- 
zation  routine  which  starts  the  program.  By  including  limiting  valves  calcu- 
lated in  the  high  temperature  limit,  the  executiom  time  was  improved  by 
nearly  an  order  of  magnitude.  Finally,  a correction  routine  was  included 
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in  the  program  to  check  for  and  correct  for  roundoff  error  due  to  the  large 
number  of  iterations  required  to  complete  the  calculation. 
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Having  developed  an  efficient  and  accurate  computer  routine,  W.  P. 

Latham  and  A.  L.  Smirl  generated  the  final  theoretical  curves  comparing  the 

experiments  of  Smirl  and  Matter  to  the  theory  of  A.  Elci,  et  al.  The  results 

were  presented  in  the  paper  by  A.  Elci,  A.  L.  Smirl,  J.  C.  Matter  and  M.  0. 

Scully.1  A copy  of  this  paper  is  included  as  Appendix  A of  this  report. 

By  the  time  the  paper  by  Elci,  et  al.1  was  completed,  additional  optical 

experiments,  in  which  the  band  gap  was  tuned  by  pressure,  had  been  performed 
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by  H.  M.  van  Driel,  J.  S.  Bessey  and  R.  C.  Hanson.  The  computer  code  was 

modified  to  include  these  pressure  effects  and  the  corresponding  theoretical 

curves  were  generated  by  Latham  and  Smirl.  The  agreement  between  these  new 

experiments  and  the  proposed  model  was  excellent.  The  pressure  experiment 

is  reviewed  in  the  preprint  contained  in  Appendix  B of  this  report. 

Subsequently,  J.  S.  Bessey,  B.  Bossachi,  H.  M.  van  Driel  and  A.  L.  Smirl3 

completed  single  pulse  transmission  studies  in  which  the  incident  optical 

pulsewidth  was  varied.  (A  preprint  of  this  paper  is  included  in  Appendix  C 

of  this  report.)  The  pulsewidth  experiment  indicated  the  first  disagreement 

between  theory  and  experiment.  Whether  this  disagreement  was  due  to  the 

simplifying  approximations  made  in  the  model,  the  uncertainties  in  physical 

constants,  or  the  details  of  the  physical  processes  included  in  the  model 

was  investigated.  By  numerically  solving  the  equations  for  several  valves 

of  the  parameters  in  the  theory,  Latham  and  Smirl  generated  a large  library 

of  figures  in  which  theory  and  experiment  were  compared.  The  culmination  of 

this  theoretical  study  was  the  paper  by  W.  P.  Latham,  A.  L.  Smirl,  A.  Elci, 

4 

and  J.  S.  Bessey  , which  is  presented  as  Appendix  D of  the  report.  Finally, 
a review  paper,  in  which  diffusion  is  included  in  the  model  for  the  first 


time,  was  presented  by  A.  Elci,  A.  L.  Smirl,  M.  O.  Scully,  and  C.  Y.  Leung.5 
A copy  of  this  paper  is  also  included  as  Appendix  B of  this  report. 

In  summary,  throughout  the  year,  this  project  provided  the  theoretical 
and  numerical  support  needed  to  compare  theory  with  experiment.  Members  of 
this  project  submitted  several  publications.  The  initial  paper  by  Elci,  et 

al.1  has  been  published  in  Physical  Review  B.  The  paper  by  Bessey,  et  al.3 
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has  been  submitted  for  publication.  Both  the  paper  by  Latham,  et  al.  and  Elci, 
et  al.5  were  presented  at  the  International  Conference  on  Hot  Electrons  in  Semi- 
conductors, which  was  held  at  North  Texas  State  University,  and  have  been  ac- 
cepted for  publication  by  Solid  State  Electronics.  All  of  the  members  of  this 
project  actively  participated  in  the  conference.  A.  L.  Smirl,  principal  inves- 
tigator, served  as  conference  vice-chairman. 

Exper iment 

In  the  first  year  of  this  grant,  the  experimental  effort  has  been  devoted 
to  designing,  constructing,  and  assembling  a system  that  enables  us  to  directly 
monitor  ultrafast  processes.  This  system  utilizes  an  excitation-probe  technique 
first  used  by  Shelton  and  Armstrong6  to  monitor  the  evolution  of  nonequilibrium 
processes  on  a picosecond  time  scale. 

This  excitation-probe  technique  requires  the  use  of  the  following  equip- 
ment (see  Figure  1) s 

(1)  A mode-locked  laser  producing  trains  of  ultrashort  pulses. 

(2)  An  electro-optic  shutter  to  separate  one  of  the  ultrashort  pulses  from 
the  train  of  pulses  produced  by  the  mode-locked  laser. 

(3)  A detection  system  to  monitor  the  pulse  train  and  the  shutter. 

(4)  An  apparatus  to  split  the  "switched-out"  pulse  and  temporally  delay  one 
part  with  respect  to  the  other. 

(5)  An  optical  dewar  to  mount  the  sample  and  to  control  its  temperature. 
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((>)  A detection  system  to  monitor  the  "switched-out"  pulse  both  before  and 
after  the  sample. 

This  last  year  has  been  spent  in  designing,  constructing,  assembling, 
and  testing  every  component  of  the  system  described  above  as  well  as 
other  pieces  of  peripheral  equipment. 

The  mode-locked  laser  system  consists  of  an  end  mirror,  a laser  head 
containing  a xenon  flashlamp,  a Nd:glass  rod,  a cooling  system,  an  output 
coupler,  a dye  cell  in  contact  with  the  output  coupler,  a power  supply  for 
the  flashlamp,  and  an  alignment  laser  to  ensure  the  alignment  of  the  laser 
cavity. 

This  was  the  first  part  of  the  system  to  be  assembled  and  tested.  Mode- 
locking the  laser  requires  that  the  dye  concentration  be  balanced  correctly 
with  the  power  applied  to  the  xenon  flashlamp  to  ensure  that  the  rod  lases 
only  once  during  the  symmetrically  shaped  flashlamp  pulse.  As  a consequence 
the  dye  concentration  (i.e.  the  volume  to  volume  ratios  of  1-2  dichloroethane 
and  an  Eastman  Kodak  dye,  A9860)  was  varied  from  10:1  to  3:1  while  the  laser 
output  was  monitored  at  each  concentration  as  a function  of  voltage  applied 
to  the  flashlamp.  Two  monitoring  processes  were  used.  The  output  of  the 
laser  was  monitored  directly  using  a fast  risetime  biplanar  photodiode- 
oscilloscope  combination.  Also  the  output  was  monitored  indirectly  using 
a silicon  PIN-photodiode  in  the  photoconduct ive  mode  with  an  oscilloscope, 
and  the  laser  radiation  was  observed  as  a spike  riding  on  top  of  the  sym- 
metric flashlamp  pulse.  As  a result  of  these  measurements  a 5:1  dye  solution 
was  chosen  as  the  best  match  to  our  needs.  It  is  necessary  to  point  out,  how- 
ever , that  mode-locking  is  such  a complicated  statistical  process  that  we 
have  from  time  to  time  run  into  difficulties  in  obtaining  clearly  mode-locked 
output  from  the  laser  system.  At  such  a time,  one  must  check  all  of  the 


I’OBfllbl*  causes  of  problems  (e.g.  bad  dye,  dirty  dye  cell)  until  the  problem 
is  corrected. 

The  next  piece  of  equipment  that  was  installed  was  the  electro-optic 
shutter  (see  Fig.  2)  consisting  of  two  Glan-Foucault  polarizers,  a Pockel's 
cell,  a spark  gap,  a d.c.  power  supply,  and  coaxial  cable  to  connect  the 
power  supply,  spark  gap,  and  Pockel's  cell.  Among  the  important  problems 
associated  with  the  set-up  of  the  shutter  were: 

(1)  One  would  like  the  ratio  of  light  passed  by  the  shutter  with  voltage 
applied  to  the  Pockel's  cell  compared  with  the  light  passed  by  the 
shutter  with  no  voltage  on  the  cell  to  bt  as  large  as  possible.  In 
order  to  maximize  this  ratio  (the  switch-out  ratio)  the  polarizers  must 
be  crossed  precisely  before  the  Pockel's  cell  is  placed  between  them. 

The  switch-out  ratio  was  tested  by  monitoring  the  light  passed  with  no 
voltage  applied  to  the  Pockel's  cell  and  comparing  it  to  the  light 
passed  when  a single  pulse  was  switched  from  the  pulse  train  by  the 
shutter.  This  ratio  was  measured  to  be  approximately  100:1  and  should 
not  be  confused  with  the  usual  extinction  ratio. 

(2)  The  KD*P  crystal  in  the  Pockel's  cell  must  be  aligned  very  precisely 
with  its  optical  axis  parallel  to  the  direction  of  propagation  of  the 
pulse  train.  Alignment  of  *-he  Pockel's  cell  was  ensured  by  viewing  the 
crystal  between  the  crossed  polarizers  while  it  was  illuminated  by  an 
extended,  diffuse  source  and  by  orienting  the  crystal  in  such  a manner 
that  a dark  or  opaque  area  shaped  like  a Maltese  cross  appeared  centered 
in  the  aperture  of  the  polarizers. 

(3)  Breakdown  of  the  spark  gap  must  be  controlled  in  one  of  two  ways:  either 
by  adjusting  the  gap  spacing  or  by  adjusting  the  gas  pressure.  The  pro- 
duct of  gas  pressure  and  gap  spacing  must  be  chosen  so  that  breakdown 
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Schematic  Diagram  of  the  Electro-Optic  Shutter  System. 
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occurs  far ly  in  the  pulse  train.  One  of  the  first  few  pulses  must  be 
switched  from  the  pulse  train  since  the  pulsewidth  increases  at  the  back 
of  the  train  due  to  the  various  nonlinear  processes  associated  with  the 
pulse  production. 

(4)  The  length  of  the  coaxial  cable-connector  system  between  the  d.c.  power 
supply  and  the  spark  gap  must  be  chosen  carefully  since  it  controls  how 
long  the  shutter  will  be  "open".  A problem,  which  still  has  not  been 
solved,  is  that  the  shutter  occasionally  switches  out  parts  of  two  pulses 
instead  of  one.  This  problem  is  associated  with  the  "jitter"  in  the 
overlap  of  the  voltage  applied  to  the  KD*P  crystal  and  the  pulse  train. 

When  the  voltage  overlaps  a single  pulse,  it  is  switched  out;  whereas, 
when  the  voltage  overlaps  two  pulses,  parts  of  both  are  switched  out. 

This  problem  has  been  approached  in  two  ways.  First,  the  frequency  of 
double  switchout  was  reduced  by  varying  the  spark  gap  spacing  and  the 
gas  pressure.  We  are  now  in  the  process  of  trying  a shorter  length  of 
coaxial  cable  between  the  d.c.  power  supply  and  the  spark  gap  (see  Fig.  2) . 
We  believe  that  this  will  solve  the  problem. 

(5)  The  length  of  the  coaxial  cable-connector  system  between  the  spark  gap  and 
the  Pockel's  cell  must  be  chosen  carefully  as  it  controls  the  time  delay 
between  breakdown  in  the  spark  gap  and  turn  on  of  the  Pockel's  cell. 

(6)  The  coaxial  cable  on  the  other  side  of  the  Pockel's  cell  must  be  termi- 
nated properly  or  ringing  will  occur  allowing  more  than  one  pulse  to  be 
switched  out. 

Once  we  had  the  laser  and  shutter  assembled  and  operating,  we  calibrated 
the  silicon  PIN-photodiode  detection  systems  for  monitoring  the  energy  of  the 
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switched-out  pulse.  These  energy  meters  were  constructed  "in  house"  at  a 
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cost  of  approximately  $50.00  apiece.  The  detectors  were  checked  for  linearity 
and  found  to  linear  over  a 1-9  volt  range.  The  absolute  calibration  of  the 
detectors  was  obtained  by  comparison  of  their  readings  with  those  obtained 
from  a Scientech  Laser  Energy  Meter  and  sensitivities  were  obtained  for 
each  of  the  detectors. 

In  addition,  two  detector  systems  that  monitor  the  pulse  trains  and  the 
electro-optic  shutter  have  been  installed  and  tested.  The  first  consisted 
of  a Hamamatsu  R617U  biplanar  photodiode  in  combination  with  a Tektronix  7904 
oscilloscope  with  a fast  write  option.  This  combination  (mainly  the  oscillo- 
scope) proved  to  be  less  than  optimal  for  our  purposes  and  the  oscilloscope 
is  presently  being  used  on  the  other  project  funded  by  this  grant.  The 
second  system  consists  of  an  ITT  F4000  biplanar  photodiode  in  combination 
with  a Tektronix  519  oscilloscope.  This  system  works  satisfactorily  although 
we  are  now  redesigning  the  housing  on  the  biplanar  photodiode  to  reduce  its 
ringing. 

The  pulse  splitting  and  delay  optics  (see  Fig.  3)  consist  of  a beam 
splitter,  a stationary  90°  prism  (mounted  on  a very  accurate  Klinger  Scientific 
translation  stage  providing  a resolution  of  0.01  mm  and  an  accuracy  of  3 mm 
over  the  40  mm  displacement  allowed)  and  two  mirrors  on  orthoganally  adjust- 
able mounts.  Two  of  the  problems  associated  with  this  piece  of  equipment  are: 
(1)  The  paths  traveled  by  the  different  pulses  separated  by  the  beamsplitter 
must  be  measured  precisely  so  that  the  path  containing  the  translation 
stage  can  be  adjusted  from  a few  picoseconds  shorter  than  the  fixed  path 
to  many  picoseconds  longer  than  the  fixed  path.  This  can  be  done  only 
coarsely  by  mechanical  measurement  techniques,  and  zero  delay  must  be 
determined  by  more  accurate  optical  techniques. 
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(2)  The  two  pulses  must  overlap  at  the  sample  surface  requiring  tests  to  be 
made  before  and  after  each  experiment  utilizing  burn  patterns  on  a piece 
of  photographic  emulsion  positional  in  the  plane  of  the  sample. 

The  optical  dewar  itself  consists  of  a relatively  inexpensive  brass  con- 
tainer with  cold  finger  and  a brass  jacket  with  four  optical  ports  covered 
by  3 inch  quartz  optical  flats.  The  dewar  is  evacuated  by  a "roughing"  pump 
and  is  capable  of  lowering  the  sample  temperature  to  approximately  100  K. 

The  dewar  sits  on  a NRC  lab  jack  and  is  capable  of  being  moved  vertically 
through  3 inches. 

Finally,  the  latest  piece  of  equipment  to  be  installed  and  tested  is  a 
flowing  dye  cell  for  the  mode-locked  laser.  The  old  arrangement  consisted 
of  a dye  reservoir,  a dye  cell,  and  a hypodermic  syringe  all  connected  by 
teflon  tubing.  Every  few  laser  firings  new  dye  had  to  be  pulled  into  the 
dye  .ell  by  adjusting  the  syringe.  The  new  flowing  system  consists  of  the 
dye  reservoir,  dye  cell,  and  a peristaltic  pump  which  pulls  the  dye  solution 
through  the  cell  at  a rate  of  0.2  cc /minute. 

Once  the  entire  system  was  assembled  and  each  of  the  individual  components 

had  been  shown  to  be  functioning  properly,  it  was  necessary  to  test  the  system 

as  a unit.  This  was  done  by  repeating  the  single-pulse  and  excitation-probe 

7 8 9 

experiments  performed  on  germanium  by  Kennedy,  et  al.  ' and  Shank  and  Auston  . 
These  experiments  have  been  completed  and  the  results  are  shown  in  Fig.  4 and 
Fig.  5,  and  the  results  are  in  good  agreement  with  previous  experiments. 

In  conclusion,  experimentally,  the  first  year  of  this  grant  has  been 
devoted  to  developing  and  ensuring  the  operational  integrity  of  the  basic 
system.  Each  component  of  the  system  has  been  installed  and  tested  rigor- 
ously. The  system  as  a whole  has  been  tested  by  comparing  data  obtained  using 
the  present  system  with  data  obtained  at  other  laboratories. 


INCIDENT  ENERGY  (QUANTA  @ 1.06  ^m) 


Versus  Time  Delay  of  Probe  Pulse  with  Respect  to  Excitation  Pulse  at  Temperature 


STRESS  AND  ELECTRIC  FIELD  BIASING  OF  OPTICAL 


PROPERTIES  OF  SEMICONDUCTORS 

For  p-type  GaAs  and  Ge,  the  optical  absorption  at  10.6  pm  is  due  pri- 
marily to  direct  intervalence  band  transitions  between  the  heavy  hole  and 
light  hole  valence  bands.  These  direct  transitions  depend  on  the  effective 
mass  (energy  band  structure) , the  occupancy  of  the  various  energy  levels 
(Fermi-Dirac  distribution  function)  and  the  value  of  the  matrix  element  of 
the  crystal  momentum  between  Block  states  from  the  two  bands. 

A low-intensity,  cw  C02  laser  is  being  used  to  monitor  the  change  in 
the  optical  absorption  of  these  semiconductors  as  a function  of  applied 
stress.  Since  the  effect  of  uniaxial  stress  on  the  band  structure  of  these 
materials  is  well  studied,  these  low- intensity  experiments  should  isolate 
and  calibrate  the  effects  of  stress  on  the  optical  absorption  at  10.6  pm 
without  altering  the  carrier  temperature  in  either  band.  The  equipment  for 
these  measurements,  as  shown  in  Fig.  6,  has  been  designed  and  constructed 
in-house. 

Much  of  the  interest  in  p-type  Ge,  however,  is  for  devices  that  operate 
in  the  region  of  saturated  absorptivity . At  high  illumination  intensities 
two  significant  perturbations  to  the  Fermi  statistics  of  the  two  bands  should 
result:  (1)  A large  number  of  direct  absorption  events  will  occur,  leading 

to  a redistribution  of  the  electrons  and  holes  between  the  heavy  hole  and 
light  hole  bands  and  (2)  the  carriers  in  both  bands  will  be  heated.  These 
changes  will  be  monitored  by  measuring  the  optical  transmission  and  elec- 
trical conductivity.  In  order  to  provide  high  intensity  illumination  of 
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Fig.6. Small  signal  transmission  scheme. 
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the  semiconductors  for  saturated  piezoabsorptivity  studies,  a pulsed  CC>2 
laser  was  designed  and  built.  This  double-discharge  TEA  laser  has  a 60  cm 
active  region  using  Rogowski  profile  electrodes.  The  flowing  gas  is  mixed 
in  a gas  manifold  that  allows  for  variations  and  innovations  in  the  gas  mix- 
ture. The  gas  is  pre-ionized  by  two  flashboards  about  1.5  usee  prior  to  the 
main  discharge  of  about  35 j.  The  peak  optical  power  output  is  in  the  range 
of  1-5  Mw.  The  width  of  the  optical  pulse  can  be  chosen  from  one  of  three 
ranges.  First,  by  adding  a cw  gain  medium  within  the  optical  cavity,  one 
obtains  the  gain-switched  output  from  the  laser,  as  shown  in  Fig.  7.  The 
FWHM  pulse  width  is  about  80-90  nsec.  Secondly,  if  the  cw  section  is  replaced 
by  a p-type  Ge  passive  mode-locker,  a mode  locked  train  can  be  obtained  with 
a single  pulse  width  of  about  1 nsec.  Also,  pulses  with  continuously  variable 
pulse  width  in  the  range  of  70-150  psec  can  be  obtained  from  the  free-induction 
decay  scheme  shown  in  Fig.  8.  Here  the  smoothed,  gain-switched  pulse  is 
truncated  suddenly  by  optically  induced  plasma  in  N2  gas.  The  beam  is  then 
spatially  filtered  and  passed  through  a cell  of  hot  C02  gas.  The  pulse  is 
absorbed  and  a short  coherent  pulse  is  emitted.  The  width  of  the  pulse  is 
determined  by  the  phase  relaxation  time,  which  can  be  varied  as  a function 
of  the  gas  pressure  in  the  hot  (X>2  cell. 

Other  pieces  of  equipment  which  have  been  designed  and  built  in-house 
include  optical  tables,  LS 2 Dewars,  a 10  w cw  (X>2  laser,  stress  apparatus, 
a variable  gas-filled  optical  attenuator,  and  many  pieces  of  optical  mounting 
and  translating  equipment.  Also,  a ratiometer  for  small  signal  transmission 
measurements  and  a dual-channel,  gated  peak-sample-and-hold  for  the  acquisi- 
tion of  pulse  data  were  designed  and  built  in-house  by  Mr.  D.  Maxson. 

The  saturated  pulse  transmission  experiment  is  shown  schematically  in 
Fig.  9.  Preliminary  data  shows  that  amplified  short  pulses  will  be  needed 


g.7  Typical  smoothed  gain- switched  CCX,  pulse. 


Fig.  8.  Short  pulse  generation  by  free  - induction  decay. 


Fig.  9.  Saturated  pulse  transmission  experiment. 
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to  obtain  high  enough  power  density  to  saturate  the  sample  without  optical 
surface  damage. 

In  sumnary,  the  equipment  for  both  the  small  signal  transmission  and 
the  saturated  pulse  transmission  has  been  designed  and  constructed  and  the 
system  is  essentially  complete.  The  techniques  are  currently  being  refined 
to  produce  more  accurate  data  and  results  are  expected  in  the  near  future. 
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A fir»(-pnrKiple»  theory  n developed  for  the  feneration  and  the  lutneqoent  transient  behavior  of  dense 
electron-hole  plasmas  produced  in  germanium  by  intense  picosecond  optical  pulses  Espenmental  dau  are 
discusacd  and  compared  with  the  theory  It  is  shown  that  the  valley  structure  of  the  Ce  conduction  bard, 
phonons,  and  plasmons  play  significant  roles  in  the  generation  and  the  temporal  evolution  of  the  plasma,  t he 
agreement  between  the  theory  and  eapenments  is  good  Some  predictions  of  the  theory  are  discussed 


l INTRODUCTION 

The  development  of  mode- locked  lasers'  that 
produce  optical  pulses  of  a tew  picoseconds  dura- 
tion and  peak  powers  of  2l0>  W has  made  possible 
recent  experimental  investigations  of  ultrafast 
electronic  processes  in  semiconductors.1"’  In 
these  experiments,  the  general  method  is  to  sub- 
ject an  intrinsic  semiconductor  first  to  an  intense 
picosecond  (excitation)  pulse2  of  frequency  w0  such 
that  ft w0  is  greater  than  the  band-gap  energy  of  the 
material  and  to  measure  its  transmission;  then 
after  some  delay,  to  subject  the  same  sample  to 
a second  (probe)  pulse  of  the  same  frequency  but 
lower  Intensity  and  measure  its  transmission. 

The  transmission  properties  of  the  semiconductor 
are  significantly  altered  when  the  excitation  pulse 
creates  a large  (ree-carrier  density,  and  the 
transmission  of  the  probe  pulse  yields  information 
concerning  the  temporal  evolution  of  the  free-car- 
rler  distribution.  The  ultrashort  duration  of  the 
laser  pulses  makes  it  possible  to  follow  this  tem- 
poral evolution  on  a picosecond  time  scale. 

In  this  paper  we  present  a first- principles, 
theoretical  treatment  accounting  for  the  observed 
behavior  in  the  ex  it  ion  and  probe  pulse  meas- 
urements in  terms  o.  the  energy-band  structure 
and  other  well-known  parameters  of  germanium 
(Ce)  We  note  at  the  outset,  however,  that  only 
certain  of  the  myriad  of  possible  electronic  inter- 
actions are  included  in  the  present  work.  Future 
papers  in  this  series  will  deal  with  other  possible 
electron  processes  such  as  diffusion  and  Kane- like 
excitations. 

In  the  remainder  of  this  introductory  section, 
we  briefly  review  the  experimental  method  and 
results  (8ec.  I A),  summarize  the  pertinent  ger- 
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man  turn  energy-band  structure  (Sec.  IB),  and 
qualitatively  discuss  the  physical  processes  which 
occur  during  and  after  the  nonlinear  light  absorp- 
tion by  germanium  (Sec.  IC).  We  expound  on  the 
theory  in  Sec.  H.  This  elaboration  proceeds 
through  a discussion  of  electron-hole  distributions 
(Sec.  QA),  direct  absorption  (Sec.  IIB),  free-car- 
rler  absorption  (Sec.  UC),  phonon- assisted  relax- 
ation (Sec.  IID),  plasmon- assisted  recombination 
(Sec.  IIE),  rate  equations  (or  absorption  and  trans- 
mission (Sec.  QF),  and  integration  of  the  dynami- 
cal equations  (Sec.  DG).  The  final  section  (Sec. 

HI)  of  the  main  text  contains  the  comparison  of 
theory  with  experiment  and  our  conclusions.  Sev- 
eral appendixes  include  detailed  calculations  which 
allow  us  to  make  critical  simplifying  assumptions 
in  our  theoretical  model.  Appendix  A considers 
intervalley  transitions  via  phonon  emission;  Appen- 
dix B discusses  Coulomb  thermalization;  Appendix 
C derives  radiation  absorption  rates;  Appendix  D 
expands  on  phonon-assisted  relaxation.  Finally, 
in  Appendix  E,  we  present  some  details  concern- 
ing the  experimental  method  and  apparatus. 

A.  Experimental  method  and  data 

A mode- locked  Nd  .glass  (JTwu*1.17  eV)  laser  pro- 
duces a train  of  pulses;  each  pulse  is  typically 
5-10  psec  in  duration  and  has  -10”  quanta.  Using 
a laser- triggered  spark  gap  and  electro-optical 
shutter,  a single  pulse  can  be  Isolated  from  this 
train  of  pulses* 

In  the  first  experiment,  this  single  pulse  is 
passed  through  a variable  attenuator  and  focused 
on  an  intrinsic  5.2- pm-thick  Ce  sample  [Fig.  1(a) j. 
The  focused  spot  is  approximately  250  pm  in  diam- 
eter. The  transmission  of  this  pulse  is  then  mon- 
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FIO  1.  Schematic  dia- 
gram of  picosecond  ex- 
periments for  (■)  ssturs- 
tton  and  (b)  excitation 
pulse/probe  studies. 


itored.  The  transmission  data  as  a function  of  in- 
cident number  of  quanta,  for  two  initial  sample 
temperatures  (77  and  297  "K),  are  shown  in  Fig.  2. 
At  low  pulse  energies,  transmission  follows  Beer’s 
law.’  As  the  incident  number  of  quanta  increases, 
transmission  increases  and  eventually  saturates. 
Near  5 x 10“  quanta  samples  are  damaged. 

In  the  second  experiment  [Fig.  1(b)],  the  sample 
is  first  irradiated  by  an  excitation  pulse  of  suffic- 
ient energy  to  cause  enhanced  transmission.  It  is 
then  followed,  at  a later  time  /,,  by  a probe  pulse 
which  is  about  5^  as  intense  as  the  excitation  pulse. 
The  excitation  and  probe  pulses  are  derived  from 
a single  pulse  as  described  in  Appendix  E.  Probe 
pulse  transmission  data  versus  delay  time  tt  are 
shown  in  Fig.  3 for  initial  sample  temperatures 
of  77  and  297  “K  and  for  a fixed  excitation  pulse  of 
about  10*’  quanta.  There  is  a sharp  Increase  in 
transmission,  a spike,  in  the  region  of  temporal 
coincidence  of  the  probe  and  excitation  pulses  (at 
0).  The  spike  is  due  to  the  scattering  of  the 
excitation  pulse  into  the  probe  pulse  by  an  index 
grating  formed  by  the  interference  of  these  two 
coherent  pulses  in  the  crystal.’’*  In  this  paper, 


FIG.  3.  Probe  pulse  transmission  data.  Relative 
transmission  units  are  given  by  the  ratios  of  the  ex- 
citation transmission  to  probe  transmission  for  each 
delay,  normalized  such  that  the  peak  of  the  77 *K  curve 
la  unity. 


we  will  be  concerned  with  the  probe  transmissions 
at  delay  times  longer  than  the  pulse  width.  In  this 
region,  probe  transmission  increases  relatively 
rapidly,  reaching  a maximum  near  t4~  100  peec, 
thereafter  decreasing  slowly.  Figure  4 shows  the 
probe  pulse  data  for  different  excitation  pulse  in- 
tensities and  for  an  initial  sample  temperature 
of  77  LK.  Ratios  labeling  these  curves  are  the  nor- 
malized relative  transmissions  of  excitation  pulse 
to  probe  pulse.  The  details  of  the  experiments  are 
further  discussed  in  Appendix  E. 

B Review  of  Ge  energy  (fracture 

The  energy-band  structure  of  Ge  is  well  known ,* 
and  is  shown  in  Fig.  5.  The  significant  features  of 
the  conduction  band  are  the  locations  of  the  con- 
duction-band valleys.  The  minima  of  these  valleys 
are  quite  close  in  energy.  The  minimum  located 
at  T is  separated  from  the  top  of  the  valence  band 
by  0.805  eV  at  300 °K.  This  separation  is  0.889  eV 
at  77  ‘K.  The  indirect  gap  located  at  L has  a sep- 
aration of  0.664  eV  at  300  rK  and  a separation  of 
0.734  eV  at  77 “K.  The  band  minimum  located  near 
X is  0.18  eV  higher  than  the  minimum  at  L.  There 
are  three  ether  valleys  such  as  the  one  in  the  [111] 
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FIG.  2.  Excitation  pulse  transmission  (saturation) 
data. 


FIG.  4.  Probe  pulse  transmission  data  for  different 
excitation  pulse  Intensities. 
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FIG.  5.  Ge  energy-band  (eV)  at  rue  tu  re  at  300  *K. 


direction  whose  minima  are  located  exactly  at  the 
Brlllouln-zone  boundary.  There  are  live  other 
valleys  like  the  one  In  the  [100]  direction. 10  The 
constant  energy  surfaces  In  these  valleys  are 
elongated  ellipsoids.  The  density-of-states  ef- 
fective mass"  for  the  L valleys  is  me*  0.22m.* 

The  density-of-states  effective  mass  of  the  cen- 
tral conduction- band  valley  Is  taken  to  be*  m0 
-0.1  m at  300 °K  and  0.04m  at  77 “K. 

The  valence  band  maximum  is  at  T.  Two  of  the 
valence  bands,  which  are  degenerate  throughout  the 
Brlllouln  zone  when  spin- orbit  coupling  is  ne- 
glected, we  call  heavy-hole  bands,  while  the  third 
Is  termed  the  light- hole  band.  At  r,  the  light- 
hole  band  is  lowered  by  spin-orbit  coupling  with 
respect  to  the  heavy-hole  bands  by  about  0.28  eV. 
Except  for  a small  region  near  the  center  of  the 
Brlllouln  zone,  the  structure  Introduced  by  the 
spln-orblt  coupling  in  the  heavy-hole  bands  is 
relatively  minor.  We  will  treat  the  two  heavy- hole 
bands  as  degenerate.  The  effective  mass  of  the 
two  heavy- hole  bands  in  the  central  region  of  the 
Brlllouln  zone  is  mk  = 0.34m.*  We  will  ignore  the 
light-hole  band  owing  to  its  small  contribution  to 
the  hole  density  of  states  and  its  large  separation 
from  the  T conduction-band  edge,  which  makes  its 
contribution  to  light  absorption  negligible  In  the 
current  experiments. 

Let  us  remark  that  an  enormous  simplification 
results  In  the  algebra  of  the  subsequent  sections  If 
we  treat  L and  X valleys  in  an  equivalent  manner 
and  assume  parabolic  energy  bands  with  the  effec- 
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tlve  masses  mentioned  above.  The  quantitative 
results  are  not  affected  much  by  such  an  approxi- 
mation because  these  valleys  enter  into  the  absorp- 
tion properties  mainly  through  their  density  of 
states.  This  approximation  results  in  the  simpli- 
fied band  structure  shown  In  hig.  6,  in  which  the 
four  "effective”  L valleys  and  six  X valleys  are 
replaced  by  ten  equivalent  parabolic  valleys.  In 
each  of  these  new  valleys  the  effective  mass  (or 
electrons  is  me  = 0.22m.  These  ten  equivalent 
valleys  are  at  some  indirect  energy  gap  £c,  which 
Is  some  average  of  the  actual  indirect  gaps  EL  (at 
L)  and  Ex  (near  X)  The  average  we  use  is  £- 
■ ♦ 8£z).  The  new  valleys  are  assumed  to  be 

centered  In  the  Brlllouln- zone  boundaries.  In  this 
approximate  band  model,  the  T valley  and  the  two 
degenerate  valence  bands  are  also  taken  as  para- 
bolic, with  the  effective  masses  mo«0.1m  and 
m,*0.34m,  respectively. 

C.  Phynici  of  nonlinear  light  abaorption  in  Gc 

When  an  excitation  pulse  is  incident  on  a Ge 
sample,  a portion  is  reflected.  The  unreflected 
portion  of  the  excitation  pulse  enters  the  bulk  of 
the  crystal  where  most  of  it  is  absorbed.  The 
light  entering  the  bulk  of  the  crystal  is  absorbed 
primarily  by  two  processes.  In  the  first  process, 
“direct  optical  absorption,"  a quanta  of  light  from 
the  excitation  pulse  is  absorbed  by  an  electron 
which  makes  a transition  from  near  the  top  of  the 
valence  band  to  the  conduction- band  valley  near  T, 
leaving  behind  a hole  in  the  valence  band.  Such  a 
process  Is  allowed  since  the  energy  of  the  light 
quanta  is  greater  than  the  direct  band  gap  £„. 
Once  there  are  electrons  In  the  conduction  band 
and  holes  in  the  valence  band,  "tree-carrier  ab- 
sorption” is  possible.  Free-carrier  absorption  is 
the  process  whereby  an  electron  in  any  of  the  con- 
duction-band valleys  makes  a transition  to  a state 
higher  in  that  same  valley  by  means  of  the  simul- 
taneous absorption  of  a quanta  of  light  and  the  ab- 
sorption or  emission  of  a phonon  (optical  or  acous- 
tical). An  identical  process  occurs  for  holes  in  the 
valence  band.  The  rate  for  direct  absorption  is 
usually  larger  than  that  for  free-carrier  absorp- 
tion; however,  the  rate  at  which  direct  absorp- 
tion events  can  occur  decreases  as  the  number  of 
occupied  electron  states  in  the  T valley  increases. 
On  the  other  hand,  the  rate  for  free-carrier  ab- 
sorption events  increases  as  the  number  of  elec- 
trons (holes)  In  the  conduction  (valence)  band  in- 
creases. 

Indirect  phonon- assisted  interband  absorption 
processes,  which  Involve  the  transition  of  an  elec- 
tron from  the  valence  band  near  T to  either  the  L 
or  X conduction-band  valleys  by  the  simultaneous 
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absorption  of  a photon  and  Die  absorption  or  emis- 
sion of  a phonon,  arc  not  Important  In  our  problem. 
The  probability  that  an  electron  will  reach  the  L 
or  X valley  by  means  of  a real  optical  transition 
to  the  r valley  followed  by  a phonon-assisted  scat- 
tering to  one  of  these  side  valleys  is  much  greater 
than  the  probability  that  the  electron  will  reach  the 
same  valley  by  means  of  a second-order  phonon- 
assisted  transition.  For  this  reason,  Indirect  op- 
tical-absorption events  may  be  ignored.11 

Because  of  the  small  number  of  electron  states 
available  in  the  T valley,  one  might  conclude  that 
the  direct  transitions  are  saturated  at  relatively 
low  pulse  energies.  In  fact,  this  is  not  the  case 
since  an  electron  in  the  T valley  will  rapidly  emit 
or  absorb  a phonon  and  make  a transition  to  one 
of  the  valleys  at  X or  L . A rate  calculation  in  Ap- 
pendix A shows  that  the  rate  for  electron  transi- 
tions from  the  T valley  to  the  L valley  is  greater 
than  10M  sec'1  since  fiu'0>EttI+  max(ffnB<)  and 
since  the  electron- phonon  coupling  constants  in 
Ge  are  relatively  large  Here  ELtX  refers  to  the 
indirect  gaps  and  flu<  to  the  phonon  frequency  of 
mode  p and  momentum  3 Thus  electrons  are 
emptied  from  the  central  valley  to  the  side  valleys 
at  a rate  that  is  larger  than  the  direct  absorption 
rate,  and  any  decrease  in  the  number  of  states 
available  for  direct  absorption  is  ultimately  de- 
termined by  the  buildup  of  the  electron  populations 
in  the  X and  L valleys. 

Owing  to  the  narrow  bandwidth  of  the  excitation 
pulse,  only  narrow  regions  of  states  in  the  central 
valley  of  the  conduction  and  valence  bands  are  opti- 
cally coupled  by  direct  transitions.  Thus  nonequil- 
ibrium  carriers  might  be  thought  to  occupy  very 
localized  regions  within  the  conduction  bands  and 
valence  band,  respectively.  However,  because 
each  carrier  moves  in  the  screened  Coulomb  field 
of  the  other  carriers,  carrier- carrier  scattering 
events  will  occur.  These  events  Include  electron- 
electron,  electron-hole,  and  hole-hole  collisions. 
The  rate  for  such  events  is  large.  These  collisions 
ensure  that  the  electron  and  hole  distributions  will 
be  Fermi- like.  They  also  ensure  that  the  Fermi 
distribution  for  holes  and  the  Fermi  distribution 
for  the  electrons  will  reach  a common  tempera- 
ture which  is  in  general  different  from  the  lattice 
temperature.  Rate  calculations  in  Appendix  B 
using  two  non-Ferml-like  distributions  indicate 
that  if  they  are  not  Fermi-like,  their  lifetimes 
are  quite  short  compared  to  the  inverse  of  the 
direct  optical -absorption  rate.  For  this  reason 
we  take  the  electron  and  hole  carrier  distributions 
to  be  Fermi- like  at  all  times. 

Electrons  located  high  in  a conduction- band  valley 
mav  relax  within  that  valley  by  phonon  emission. 
Holes,  of  course,  may  also  relax  by  emitting  pho- 


nons. The  effect  of  this  relaxation  is  to  reduce  the 
electron  temperature  and  increase  the  lattice  tem- 
perature. The  increase  in  lattice  temperature  due 
to  these  phonon-assisted  intraband  relaxation  pro- 
cesses is  significant  only  at  very  large  pulse  en- 
ergies. 

It  is  important  to  notice  that  the  total  number  of 
carriers  is  unchanged  by  free- carrier  absorption 
or  phonon-assisted  relaxation  within  the  valleys. 
These  processes  serve  only  to  elevate  or  reduce, 
respectively,  the  distribution  temperature.  Only 
direct  absorption  will  Increase  the  number  of 
carriers.  Recombination  processes  serve  to  re- 
duce the  carrier  number. 

The  recombination  processes  can  be  divided  into 
two  general  categories:  radiative  and  nonradiative. 
In  the  present  work,  nonradiative  recombination 
is  much  more  important  in  reducing  the  carrier 
number  than  radiative  recombination.  As  the  car- 
rier density  builds  up  (as  a result  of  direct  absorp- 
tions) the  plasma  frequency  of  these  carriers  in- 
creases. At  sufficiently  large  plasma  frequencies, 
an  electron  in  the  T valley  can  recombine  with  a 
hole  near  the  top  of  the  valence  bands  via  emission 
of  a plasmon.  Normally,  electrons  near  the  con- 
duction-band edge  (at  T)  can  recombine  with  holes 
via  emission  of  plasmons  only  if  the  plasma  fre- 
quency w,  is  larger  than  the  direct  gap  frequency 
E0/H  [see  Eq.  (53)].  However,  in  our  case,  the 
plasma  resonance  is  considerably  broadened  due 
to  the  nature  of  the  direct  absorption  and  subse- 
quent scatterings  between  and  within  the  T valley 
and  the  L-X  valleys.  Therefore,  plasmon- assisted 
recombinations  can  occur  at  plasma  frequencies 
lower  than  E0/lf.  As  the  electron  and  hole  popula- 
tions build  up,  the  plasmon- assisted  recombina- 
tion rate  becomes  comparable  to  the  direct  absorp- 
tion rate,  causing  the  carrier  number  to  saturate 
at  plasma  frequencies  much  less  than  E0/l f.  These 
plasmon- assisted  recombinations  can  be  as  im- 
portant in  determining  the  final  number  of  conduc- 
tion electrons  as  the  saturation  in  the  total  number 
of  available  electron  states  in  the  conduction  band 
due  to  the  buildup  of  the  electron  population. 

The  collective  plasma  oscillations  have  a life- 
time that  is  short  compared  to  a picosecond.  The 
energy  lo6t  in  the  decay  of  the  plasma  oscillations 
is  rapidly  transferred  to  single  electron  and  hole 
states  and,  thus,  ultimately  increases  the  temp- 
erature of  the  carrier  distribution.  Thus  the  plas- 
mon-assisted  recombinations  both  limit  the  num- 
ber of  carriers  and  raise  the  electron-hole  temp- 
erature . 

Radiative  recombination  may  be  of  two  types: 
direct  or  indirect.  The  recombination  of  an  elec- 
tron in  the  T valley  of  the  conduction  band  with  a 
hole  in  the  valence  band  by  means  of  emission  e i 
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a photon  Is  termed  direct;  the  recombination  of  an 
electron  in  the  L or  X valley  with  a hole  In  the  val- 
ence band  by  means  of  simultaneous  emission  of  a 
photon  and  emission  (or  absorption)  of  a phonon  is 
termed  indirect.  Both  processes  occur  on  time 
scales  larger  than  nanosecond,  and  are  not  im- 
portant in  our  problem. 

The  diffusion*  of  carriers  from  the  interaction 
region  (focused  spot  size  times  sample  thickness) 
also  reduces  the  carrier  number.  Preliminary  dif- 
fusion (Boltzmann  equation)  calculations  (to  be  pub- 
lished) indicate  that  the  number  of  carriers  leaving 
the  interaction  region  on  a picosecond  time  scale 
Is  small.  In  any  case  we  are  able  to  account  (or 
our  experimental  results  without  invoking  diffusion. 

In  summary,  the  following  rapid  processes  are 
invoked  In  explaining  the  passage  of  our  1.0G- 
(irn  psec  pulses  through  Ge:  particle-particle  scat- 
tering, electron-phonon  scattering,  direct  and  free- 
carrler  absorption,  ard  plasmon-assisted  inter- 
band electron- hole  recombinations.  In  the  cal- 
culation, radiative  recombination  and  diffusion 
are  ignored  as  discussed  above. 

Physically,  the  single  pulse  transmission  of  Ge 
as  a function  of  incident  pulse  energy  (see  Fig.  2) 
can  be  accounted  for  by  direct  interband  transitions 
followed  by  heating  of  the  electron  distribution  as 
follows:  When  the  excitation  pulse  enters  the  crys- 
tal, it  is  absorbed  by  direct  optical  transitions, 
creating  a large  number  of  electrons  in  the  central 
valley  of  the  conduction  band.  These  electrons  are 
quickly  scattered  to  the  side  valleys  by  phonons. 
Particle- particle  scattering  events  ensure  that  the 
carrier  distributions  are  Fermi- like  and  that  both 
electron  and  hole  distributions  have  the  same  temp- 
erature. The  transmission  Initially  rises  due  to 
the  saturation  of  the  available  optically  coupled 
states  as  a result  of  direct  absorption.  Further 
Increase  In  the  transmission  for  intensities  larger 
than  10“  quanta  is  hindered  as  the  electrons  are 
heated  and  removed  from  these  states  by  free- 
carrler  absorption  and  plasmon-assisted  recom- 
bination For  these  excitation  pulse  intensities, 
electrons  m the  T valley  fall  to  the  hole  pockets 
in  the  valer-.e  bands  via  plasmon  emissions  at  a 
rate  comparable  to  light  absorption  thus  limiting 
the  carrier  density. 

After  the  passage  of  the  excitation  pulse,  the 
interaction  region  of  the  sample  contains  a large 
number  of  carriers  with  a high  electron- hole  temp- 
erature. Plasmon-assisted  recombinations  are 
essentially  turned  off  when  the  excitation  pulse 
has  passed  This  is  discussed  in  detail  in  Sec. 

II  C.  As  time  progresses,  the  distribution  temp- 
erature la  reduced  by  phonon- assisted  intraband 
and  Intravalley  relaxation;  the  probe  pulse  inter- 
rogates this  evolution  of  the  distribution,  since  It 


is  a sensitive  measure  of  whether  the  optically 
coupled  states  are  available  (or  absorption  or  are 
occupied. 

The  probe  pulse  transmission  (see  Fig.  3)  can  be 
understood  in  the  following  way  After  the  pass- 
age of  the  excitation  pulse,  the  electrons  (holes) 
are  located  high  (low)  in  the  valleys  because  of 
the  high  distribution  temperature,  leaving  the 
states  that  arc  optically  coupled  available  for 
absorption.  Hence  probe  transmission  is  small. 
Later,  as  the  distribution  cools  by  means  of  intra- 
band phonon- assisted  transitions  and  carriers  fill 
the  states  needed  for  absorption,  the  probe  trans- 
mission increases.  Finally,  as  the  distribution 
temperature  cools  to  near  the  lattice  temperature, 
the  electrons  (holes)  occupy  states  near  the  bot- 
tom (top)  of  the  valleys,  and  states  needed  for  ab- 
sorption are  once  again  available;  the  probe  trans- 
mission then  decreases.  This  is  summarized  in 
pictorial  terms  in  Fig.  7. 


II.  THEORY 

A.  Electron  and  hole  dutributiuru 

In  the  Intioduction,  we  argued  that  the  electron 
and  hole  populations  in  our  problem  can  be  de- 
scribed by  two  Fermi  distributions  which  can  be 
assigned  the  same  temperature  because  of  -he 
rapid  Coulomb  collisions.  Electrons  and  holes 
have  distinct  Fermi  energies.  Thus 


1 + exp( 
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FIO.  7 Schematic  diagram  for  the  temporal  evolu- 
tion of  the  (roe  carriers  created  by  the  excitation  pulse. 
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fcljft  CXP  (Sfi^  * )]  > 

«cs **.*<£><♦«.  (lb) 

and  likewise  for  the  holes 

/„(E)  = |^1  4 exp^J^-—  , 0s5Jf(t)<  + «> . 

(lc) 

Here,  all  energies  are  normalized  with  respect  to 
ha'u,  where  u>„  is  the  circular  frequency  of  the 
optical  pulse.  The  subscript  rT  denotes  the  r 
valley  of  the  conduction  band,  cL,X  the  L-X  valleys 
of  the  conduction  band,  and  H the  holes  In  the  valence 
bands,  t and  li  are  the  normalized  Fermi  energies 
of  the  electrons  and  holes,  respectively;  and  S0 
and  Sc  are  the  normalized  energies  of  the  direct 
and  indirect  energy  gaps,  respectively.  In  the  re- 
duced Ge  energy-band  structure,  the  indirect  en- 
ergy gap  is  taken  to  be  <SC  ^ + 6<yx),  where 

<5  L is  the  normalized  energy  band  gap  at  L and  Sx 
is  the  normalized  energy  band  gap  at  X.  As  dis- 
cussed earlier  we  will  treat  the  L and  the  X val- 
leys as  equivalent  and  assume  that  they  have  sphe- 
roid parabolic  band  structure.  This  simplifies 
the  algebra  considerably,  but  does  not  make  much 
difference  in  the  end  results.  Also  note  that,  due 
to  the  rapid  phonon- induced  scattering  of  electrons 
between  the  central  and  the  side  conduct  ion- band 
valleys  (see  Appendix  A),  we  have  taken  the  T 
valley  and  the  L-  and  X-valley  Fermi  energies  to 
be  identical  In  the  reduced  band  structure,  elec- 
tron and  hole  energies  are  given  by 

*(i/ku„)[Ee  ***(£-  £t.j)V2m,] 


*SC  *n2(i-tL'tY/2Kv0mc  , 

(2a) 

Ser(  E)*<  1 /* u>0)(E0  * k V/2 m0) 

*«S0  + /rV/2»u)0m0  , 

(2b) 

i5„(E)  -k  *<f*/2»w0m, , 

(2c) 

and 


fined  by 

N0  = t.m<i>$/4we*  , (3) 

where  c.  is  the  high-frequency  dielectric  con- 
stant, e the  electron  charge,  and  m the  bare  elec- 
tron mass.  Let  us  note  that  due  to  its  smaller  ef- 
lecttve  mass  and  higher  band  gap,  the  T valley 
has  a negligible  density  of  electrons  compared  to 
the  L-X  valleys.  Thus,  neglecting  the  T valley 
and  considering  the  effective  mass  of  each  of  the 
four  L valleys  to  be  equal  tc  that  of  each  of  the  six 
X valleys,  the  total  electron  density  is  calculated 
as  ten  times  the  density  in  any  one  valley.  The 
electron  density  is  given  by 

n =r5/,p,(zl)  , (4a) 

where  the  p’s  and  the  z’s  are  defined  In  Table  I. 
Likewise,  the  hole  density  is  given  by 

«-f’,/,*(*1)  (4b) 

and  p,(z2)  = p,(z,)  gives  the  integral  relation  between 
( and  h. 

By  taking  the  electron  and  hole  distributions  to 
be  time  dependent  only,  rather  than  dependent  on 


TABLE  I.  Definitions  of  integrals 
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*.  r/Hu>u  , TmklT/Hu>0.  (2d) 

Here  E is  a wave  vector,  kt  Boltzmann’s  constant, 
>»n  the  effective  mass  for  the  F valley,  >ne  the  ef- 
fective mass  for  both  the  L and  the  X valleys,  and 
mt  the  effective  mass  for  the  heavy-hole  bands. 
EtiI  refer  to  the  wave  vectors  for  the  minima  of 
the  conduction  bands.  The  distributions  defined 
above  will  depend  on  time  through  c,  h,  and  T. 

The  Fermi  energies  li  and  e are  related  since  the 
total  number  of  electrons  is  equal  to  the  total  num- 
ber of  holes  at  all  times. 

In  the  present  work,  tt  is  convenient  to  measure 
the  electron  density  in  terms  of  N0,  which  Is  de- 


Ai ■ 

A j - (t/tr’\,>(2 mt  oV*)1  /l 
A,»<t/2lr,^0><2m,u),/»>,/, 

ty’tSc-tl/f 

*i'h/T 

*4  * l St  * mt)(t  - <£,)-  c|/f 

*l  • m,)(l  - So)  *M1/T 

«,*(  So-  (SG/i)  * (11/2)  InU  ,/A,))/rl 
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III.  TBAFAAT  TlAISItNT  RESPONSE  Of  SOI.  ID  STATE 


1*7 


t<o«h  lima  and  »pac«.  »<■  jr»  Ignoring  the  pulse 
pr  opag  At  Ion  pr'iblem  within  the  Interaction  region 
(■pot  area  Ilmen  the  •ample  thickness)  as  (ar  as 
the  electronic  structure  Is  concerned.  This  physi- 
cally reasonable  assumption  will  be  Justified  and 
its  consequences  will  be  detailed  in  a subsequent 
paper  dealing  with  the  pulse  propagation  problem. 
The  distributions  defined  in  Eqs.  (la)—  (lc)  should 
be  viewed  as  spatial  averages  of  the  actual  distri- 
butions over  the  Interaction  region. 

We  also  take  (he  photon  density  within  the  inter- 
action region  to  be  time  (but  not  space)  dependent: 
N(j— ’ »Mlr— ’(/).  In  addition,  we  assume  that  all 
quanta  have  the  same  wave  vector.  The  radiation 
energy  density  at  time  I is  then  given  by 

W(l)  -£flu>,  0N'0T*\t)  * HwaNaNU)  , 


The  behavior  of  the  Fermi  distributions  for  both 
electrons  and  holes  Is  determined  if  «,  h,  and  T 
are  known  at  all  times.  Rather  than  dealing  di- 
rectly with  the  distribution  functions,  we  will  de- 
rive another  set  of  equations  from  (5a)-(5c), 
which  will  describe  the  behavior  of  < , h,  and  T 
with  time.  For  this  purpose,  we  consider  the 
time  evolution  of  the  electron  density,  which  is 
given  by 


t i > 


(6a) 


where  j denotes  a summation  over  all  conduction- 
band  side  valleys.  includes  the  spin  summa- 
tion and  the  normalization  (actor,  l.e.,1’ 


PS 


2d£ 

(2  i?N0 


where  N(t)  is  a normalized  “’photon”  number  den- 
sity. Just  as  with  the  electron  density,  N[l ) should 
be  viewed  as  a soatlal  average  of  the  actual  radia- 
tion density  over  the  interaction  region. 

The  following  equations  (for  the  electron- hole 
distributions)  summarize  the  discussion  of  Sec.  I: 


We  substitute  Eqs.  (5a)  and  (5b)  into  Eq.  (6a),  and 
recognizing  that  the  free-carrier,  T-j  scattering, 
and  the  relaxation  terms  do  not  affect  the  electron 
density,  we  obtain  the  result 


dn 

dt 


(6b) 


Inspection  of  Eq.  (4a)  shows  that  n is  a (unction  of 
time  via  c and  T.  Thus  we  can  also  write 


dn  /9i 

s\  dt  J 

'»m\  dr 

dt 

‘)Td‘  ' 

\*r) ,dt  ' 

We  define 

(in 

), 

/8  n\ 
yr  *(«?). 

In  terms  of  the  integrals  defined  in  Table  I,  we 
obtain 


(5c) 

where  DA  represents  the  direct  absorptions;  FCA, 
the  free-carrier  absorptions;  T-j,  the  phonon- 
assisted  scattering  from  the  T valley  to  the  j val- 
ley; REL,  the  Intravalley  relaxation  of  electrons 
and  holes  by  means  of  emission  or  absorption  of 
phonons;  and  R,  the  direct  plasmon- assisted  re- 
combination. We  note  that  the  direct  absorption 
terms  mainly  contribute  to  an  increase  in  the  car- 
rier number.  The  free- carrier  absorption  contri- 
butions serve  to  increase  the  distribution  tempera- 
ture. The  intravalley  phonon- assisted  relaxation 
terms  decrease  the  distribution  temperature.  The 
direct  recombination  terms  decrease  the  electron 
number  and  increase  the  distribution  temperature. 
The  F-)  scattering  terms  do  not  change  the  elec- 
tron number  or  significantly  change  the  electron 
temperature;  rather,  they  serve  to  populate  the 
side  valleys  In  the  conduction  band. 


v,(<)  *rl/,pl(sl)  - rl/J<Ti(*i)  . (8) 

and 

y1(0  = iT,/Jp1U1)  + zlT4/,(pl(zl)-cr1(zl)]  . (9) 

Therefore,  combining  Eqs.  (6b)  and  (6c)  we  write 


dt 


‘H  --Y 

dt  v 3/ 


(10) 

The  fact  that  the  number  of  electrons  equals  the 
number  of  holes  requires  that 


We  define 
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-(cl/(Sl  > 

r/8„N  /8m\  l /(*n\ 

L \8T  / t ~ \8T  )t  } / \8h)r  • 


Then,  from  (11)  and  (12)  we  have 

<“> 

Expressions  for  ylk  and  yfk  are  evaluated  by  dif- 
ferentiating Eqs.  (4a)  and  (4b).  In  terms  of  the  in- 
tegrals defined  in  Table  I,  these  expressions  are 

r„U) 3 [pMi)  ~ /[Pi(*i>  " ca(*j)]  U4) 


yf,(f)«*,-  z,y,»(f)  . (15) 

To  this  point,  we  have  obtained  two  differential 
equations,  (10)  and  (13),  relating  three  unknowns: 
«,  h,  7 . We  wish  to  find  still  a third  independent 
differential  equation  relating  these  variables.  We 
accomplish  this  by  considering  the  energy  density 
of  the  electrons  and  holes.  The  energy  density  of 
the  electrons  and  holes,  normalized  by  the  photon 
energy  flu),,  and  the  number  density  N0,  is  given  by 

(16> 

and  du/dt  by 


The  energy  lost  to  the  lattice  due  to  T-)  phonon- 
assisted  scattering  events  is  negligible.  As  dis- 
cussed earlier,  plasmon-assisted  recombinations 
do  not  change  u.  Radiative  recombinations  would 
decrease  the  total  energy,  but  they  are  negligible. 
The  free-carrier  absorption  and  phonon-assisted 
tntravalley  relaxation  contributions  of  the  I"  valley 
to  the  energy  are  negligible  compared  to  those  of 
the  L and  X valleys.  The  following  equations  sum- 
marize these  statements; 


du  /•«!  /RiA  />m\ 
dt  dt  )0A  91  ) rCA  \9t  ) „gl_ 


1 (18b) 

f \ 01  'DA 

~|  ■ £ £ a,®^) 

9t  i rCA  £ j \ 91  'rcA 

«*» 

Iril  £ J \ 91  /rel 

*2  (18d) 

[ \ 9t  / HI L 

Since  u =u(« , h,  T)  we  can  write 

du  (9 u\  dt  / 9u\  dh  19 u\  dT  ...  . 

sr-(s;)„sr(n)„.«,(irL  (19*’ 

Making  use  of  the  definitions  of  y„  and  yr,,  we 
obtain 

du  [/9u\  /9u\  Idt 

5r'LWr..-U/r,/'‘Js- 

We  define 

<1W 

and 

U9d) 

Thus  combining  Eqs.  (18a)  and  (19b),  we  find  that 

,20> 

Let  us  note  that  in  the  evaluation  of  ym  and  y^, 
the  T valley  can  be  neglected  just  as  in  the  case 
of  electron  density,  and  for  the  same  reasons. 
Thus  in  terms  of  the  integrals  of  Table  1, 

yJf)*ir,',P,(*,)U*0 

♦ r'/>5c[p1Ul)-al(*1)]  (21) 


rMr(f)  *r>  /a5c[(|  + *>,(*,)  - *,»,(*,)) 

+ |r,/,[sl  (*,)•*•  ■>,(<,) «- i<,p,(S|)(l'fyM)]. 

(22) 

Let  us  also  note  that  when  the  T valley  is  neglected , 


where 
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u(t)  •r'/,[v,u,)  + Vj (*,)]•*•  r,/,s0pl(z1).  (23) 


Summarizing,  Eqs.  (10),  (13),  and  (20)  are  a 
coupled  set  of  differential  equations  describing  the 
time  evolution  of  the  Fermi  energies,  t and  h and 
carrier  temperature  T.  For  convenience,  we 
summarize  our  working  equations  below. 


, dt 


<iT  v «/-(£) 


dt 


k 


dt 


8 t 


ao) 


dh 

dt 


...dt  .,.dr 

*(t)dt  dt  • 


(13) 


+ y-^0 


dr 

dt 


du 

+ 9u 

du 

dl 

DA 

tca  9t 

(20) 


As  the  light  is  absorbed  the  carrier  density  and 
temperature  pass  through  a region  in  which  neither 
the  low-temperature  approximation  nor  the  high- 
temperature  approximation  holds;  as  a result  we 
must  evaluate  these  expressions  numerically. 
First,  we  must  turn  our  attention  to  calculating 
the  absorption  and  relaxation  rates  appearing  on 
the  right-hand  sides  of  Eqs.  (10)  and  (20) 


B Dim!  absorption 


In  this  and  Secs,  fl  C and  II  D we  will  confine 
ourselves  to  presenting  only  the  results  in  order 
not  to  burden  the  discussion  with  too  much  algebra. 
Direct  and  free- carrier  absorptions  in  Ge  are 
dealt  with  extensively  in  the  literature As 
mentioned  earlier,  a number  of  simplifying  ap- 
proximations can  be  made  in  our  problem.  For 
this  reason  we  discuss  the  appropriate  electron- 
photon  and  electron-phonon  coupling  Hamiltonians 
in  the  text,  and  present  a detailed  summary  of  the 
calculation  in  Appendix  C.  The  appropriate  dia- 
grams describing  the  perturbation  terms  used  in 
the  calculation  are  shown  in  Figs.  8 and  9.  We 
confine  ourselves  to  the  lowest-order  perturba- 
tion terms  described  by  these  diagrams,  higher- 
order  terms  being  negligible. 

The  diagrams  in  Fig  8 describe  direct  absorp- 
tion and  emission  (of  a quantum  ffw0).  The  elec- 
tron-phonon coupling  is  described  by  the  well- 
known  Hamiltonian 


(24) 

where  |slt)  refers  to  Bloch  states  (s  Is  the  band 
index);  (,(fj)  is  the  unit  photon  polarization  vector; 


cr,c.« 


cr,;.a 


FIG.  8.  Direct  optical  transitions  Solid  lines  refer 
to  electrons  or  holes,  wavy  lines  to  photons. 


55  is  the  electron  momentum  operator;  cjf  and  6, 
are  the  electron  and  photon  annihilation  operators, 
respectively.  Electron  momenta  are  generally 
much  larger  than  the  photon  momentum;  thus 

(s'E'  |e,vr(x(3)  • p|«E)  “ £,(§)  • (s'k'  |j5|  sC).  (25a) 

For  direct  absorption,  generally  the  absolute 
square  of  the  interband  matrix  element  of  the 
momentum  operator  | (cic  j p j r-Jc>  | 2 is  needed.  This 
quantity  can  be  evaluated  with  the  help  of  the  so- 
called /-sum  rule.”  Near  the  T-valley  band  edge 
(BE) 

|<cC|jJ|pt>|,a,Jwi,£0(l/"'o+l/wi,)‘  m*E0/2mht  . 

(25b) 


The  laser  frequency  t!wn  is  such  that  the  direct 
absorptions  will  involve  the  F valley  and  the  heavy- 
hole  bands  only.  The  direct  absorptions  due  to  the 
transitions  between  the  light- hole  band  and  the 
heavy-hole  bands  are  neglibible  since  such  transi- 
tions are  restricted  to  regions  away  from  the  cen- 
ter of  the  Brillouin  zone  and  thus  restricted  by  the 
Fermi  statistical  factors,  and  since  the  density 
of  the  states  of  the  light-hole  band  is  small 
Using  (24)  and  (25)  with  the  help  of  the  diagrams 
shown  in  Fig.  8 and  the  approximations  indicated 
in  Appendix  C,  one  finds 
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9u 
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DA 
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■ N(t)aD(t), 


(26) 

(27) 


\ 


FIG.  9.  Froe-carrier  transitions.  Broken  lines  refer 
to  phonons 
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and 

^1  «-N(/)a0(/),  (28) 

91  ‘da 

where 

o0(/)  = a0[l-F(*,)-F(Oi.  (29) 

and 

a = Ae*S o(2m*mi/f"<S\-6a>)l,tw  (30) 

" 3 TjP  (m0+m,)  ) 

zAt  z,,  and  F are  defined  in  Table  I.  N(l)  is  the 
normalized  photon  density  in  the  interaction  region 
mentioned  earlier.  aD(t)  is  the  direct  absorption 
rate.  (v7^/r)a„  is  the  usual  absorption  coefficient 
(per  unit  length)  observed  in  Beer’s- law  region7 
(for  the  light  frequency  iu0). 

In  Fig.  10,  aD/a0  is  plotted  as  a function  of  elec- 
tron-hole temperature  for  fixed  electron-hole  den- 
sities. The  effective  masses  and  the  energy-gap 
parameters  are  the  same  for  all  the  curves  in 
the  figure. 

We  observe  that  as  the  density  decreases,  oD/a0 
approaches  1 for  all  T: 

lim  oB  = a0.  (31) 

» -o 

Also  <ip  has  the  same  behavior  for  all  densities 
as  T-«,  which  is  to  be  expected  since  as  T 
increases  the  probability  for  the  occupancy  of 
a given  conduction  state  decreases  to  essential- 
ly zero,  thus 

lim  oo  = a0.  (32) 

r • - 

As  7- 0,  aD  shows  interesting  structure.  We  see 
that  there  are  two  critical  densities  such  that 


yir,  10  Direct  ulmorption  coefficient  vh  clcclron- 
hnlc  temperature.  The  curvet  arc  generated  for  the 
gap  p.irrimcUTS  corresponding  loTp«300'K  The 
iMtcd  and  daahed  curve*  refer  to  poaallile  trajectories 
in  lids  a,  7 space  an  the  inline  pauses  thr  n-gh  the 
■ample. 


" <*o 

i«o 


lim  o0  * 
r-o 


° 


if  >1  <R,  , 
if  n =«, , 

If  nl<n<n7, 
if  n =n,, 
if 


(33) 


where  n,  and  n,  are  given  by 


n>N°=  (34a) 


and 


(34b) 


The  existence  of  n,  and  n2  can  easily  be  inferred 
from  Fig.  5.  As  7 — 0,  electrons  collect  at  the 
bottom  of  the  conduction- band  valleys  and  holes  at 
the  top  of  the  valence  bands.  When  n<n , and 
7-0,  Er<E’  and  E„>E'\  thus  direct  absorptions 
can  take  place  unimpeded  and  oD  - a0.  When 
n>nt  and  7-0,  EH<E'  and  the  valence-band 
states  which  can  absorb  J/u>0  quanta  are  empty. 
When  n>«j  and  7-0,  Er>E’  and  E„<E' , and 
thus  the  inverse  of  absorption,  i.e.,  stimulated 
emission  of  Jfuic  quanta  occurs,  and  aD  becomes 
-ar0.  The  two  densities  n,  and  n2  are  large.  For 
typical  values  of  the  effective  masses  and  energy 
gaps,  such  as  m0=»0.1,  m,=»0.34,  mec»0.22, 

Eg~  0.7  eV  and  £o~0.9  eV,  and  fiu)0  = 1.17  eV, 
we  find  «,A'0“4  x 10‘*  cm"’  and  >Jri\’0~2  x 1021  cm-5. 
In  the  experiment,  the  maximum  electron  densities 
achieved  remain  below  n3.  Finally,  we  observe 
that  for  n<n,,  the  minimum  of  aD  occurs  at  finite 
temperatures. 

As  we  shall  see  later,  the  direct  absorptions  re- 
main dominant  compared  to  the  free- carrier  ab- 
sorptions for  all  light  intensities  used  in  the  ex- 
periment, and  therefore,  the  qualitative  features 
of  the  experimental  results  can  easily  be  inferred 
from  a diagram  such  as  the  one  in22  Fig.  10.  e. 
a 16  for  Ge  23  Thus,  for  the  effective-mass  and 
energy- gap  values  quoted  above,  «0~9x  1013  sec"*. 
The  usual  (Beer’s  law)  linear  absorption  coeffici- 
ent is  obtained  from  at0  by  dividing  by  the  speed 
of  light  in  the  absorbing  medium,  i.e., 

afc„  = (Ve^/r )a0  = 1.2x  lO'cm*1 , 
in  agreement  with  the  well-known  value. 


C.  Fm  carmi  absorption 

The  electron-phonon  interaction  is  described  by 
the  Hamiltonian 
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X <sE'  I c- ; I s Oc’.p-.i' a»5«  * H.C., 

(35) 


where  p varies  over  (he  phonon  modes,  p is  the 
Ge  mass  density,  n„;,  is  the  phonon  frequency,  Qu 
is  the  electron- phonon  coupling  coefficient,  and 
a ^ is  the  phonon  annihilation  operator.  We  note 
that  for  the  Bloch  states 


<sf'|e'>-,|sC>»55(t'-  E - ?'),  (38a) 

(sE'|p|sE)  orfllJft’lE'  - E)  , (38b) 


where  s refers  to  the  conduction  or  the  valence 
bands.  Figure  9 Illustrates  typical  free-carrler 
absorption  diagrams.  There  are  altogether  16 
such  diagrams.  With  the  help  of  such  diagrams 
and  (35) — (36b) , the  rates  (9N/9I)  and  (9u/dt)  due 
to  free- carrier  absorption  can  be  calculated  as 
summarized  in  Appendix  C.  In  those  calculates, 
two  plausable  assumptions  are  made:  (a)  the  same 
electron-phonon  coupling  coefficients  are  taken  for 
the  L and  the  X valleys;  (b)  the  same  coupling  co- 
efficients are  used  for  the  electrons  and  the  holes. 
This  procedure  is  followed  for  the  purpose  of 
simplifying  the  calculations,  and  can  be  justified 
only  a posteriori,  alter  comparison  with  the  exper- 
imental results.  Let  ^ be  the  angle  between  q' 
and  the  jth  valley  direction.  The  electron- phonon 
coupling  coefficients  are  then  as  follows.  For  the 
longitudinal  optical  phonons''' 


<?lo  *Qac°s<t>,  ; 

(37a) 

for  transverse  optical  phonons 

9ro  x Qo  si"*/  ; 

(37b) 

for  longitudinal  acoustical  phonons 

9L*(q')'f'(Aa  + A,cos ’<*,)  ; 

(37c) 

for  transverse  acoustical  phonons 

QTA(4')  * simp,  . 

(37d) 

y0  Is  the  electron-optical-phonon  coupling  con- 
stant. Ttiis  coupling  constant  is  taken  to  be  de- 
pendent on  lattice  temperature:  <^0  - 6 x 10"*  at 
300  K and  Qn  - 2 x 10“  at  77  K.  A,  and  A.  are  the 
same  deformation  potential  constants  as  in  Herring 
and  Vogt.”  Also,  we  assume  a single  reststrahl 
frequency  for  the  optical  phonons,  and  make  a 
linear  approximation  for  the  acoustical  phonons”: 

^ u) . to  “Oq,  (38a) 

^u,t*  • (38b) 

Free-carrler  absorption  depends  upon  the  elec- 
tron population.  The  T- valley  electron  population 


la  negligible  compared  to  that  of  the  L and  X val- 
leys; therefore,  the  T valley  is  ignored  in  the  FCA 
calculation.  We  also  assume  in  Appendix  C that 
each  electron  after  a phonon- assisted  photon  ab- 
sorption remains  in  the  same  valley,  since  the  in- 
tervalley  scattering  rates  among  the  L and  X val- 
leys are  small.17 

With  these  approximations  and  others  indicated 
in  Appendix  C,  we  obtain  the  following  expressions 
for  the  absorption  rates: 


91  Iron 

— I •N(t)aTCK(i)  , 

91  I rex 

where 

arcA(f)  *{  0|[1  ♦ 2(u>0/n0)r 1 1 ♦ a,  Tl) 

* Uo(*») + 5(mc/m,)1t0(zl)| 

+ «,[ti  (*j)  ♦ 5(me/m,)*,,t;l  (*,)]  , 
a,  = eaQjm*/3ir,«<.p/f<fi0  , 
a,  = 2eJmjA,u>0/3ir,«.pc^84  , 
a,  = ^e2m1kA^(2lllJl^l>0)l,3/l%r7^mpcAH,  , 
A*  = Aj+  $AfA.+  $Aj  . 


(39) 

(40) 


(41a) 

(41b) 

(41c) 

(41d) 

(41e) 


For  the  definition  of  t functions  see  Table  I.  Here, 
Tt  * kBTL/ttu>0  Is  the  normalized  lattice  tempera- 
ture. For  Ge,  » 16,  p =*  5 .3  g/cm1,  A,  =■  -3.4 
eV,  A.«17eV,  Q0  »6  x 10“  erg/cm,  /ino=»1.88 
x 10“  eV,  and  cA  o<5.4  x 10'  cm, 'sec  at  300  ’K. 

Using  these  values  we  find  a,  5=10"  sec’*,  a, 
x 10u  sec'1,  and  a,  *=10*'  sec*1.  Although  as  seems 
to  be  quite  large  compared  to  a,  and  a,,  actually 
a,  occurs  in  as  a, Tt,  and  since  generally  rt 
-10“,  a,TL  isof  the  same  order  as  a,.  Note  that 


FIG.  11.  Free-carrler  absorption  coefficient  va 
electron- hole  temperature. 
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a,,  a2,  o,  are  quite  small  compared  to  a0-®x  10*’ 
sec-4 . 

Plots  of  arCA  as  a function  of  (tlw0/*,)r(»71  for 
fixed  n.\’„  are  given  in  Fig.  11  The  behavior  of 
orCA  can  easily  be  deduced  from  (41)  for  low-  and 
high-temperature  limits.  In  the  low-temperature 
limit  (for  fixed  n) 


For  the  values  of  the  density  that  occur  in  the  ex- 
periment 3r/2A,  < 1 and  3n/2Aa  <1,  and  therefore 

♦ 5(S;),t°(2,)  =*  Bn(J;+ S)  * 96"  <43> 

and  nr^orn.  At  large  temperatures,  i.e.,  as  J 


and 


. 5/a 


f 1/2 

(44a) 

f 1 

> hr  + 5 

— 1 ) Uy 

IA  \ 

* lOOOnf.  (44b) 

Therefore,  as  temperature  increases,  for  a given 
ii,  orCA  begins  to  vary  linearly  with  T . We  also 
observe  that  in  both  temperature  limits  aTCAacn. 

D Phononaiaisted  relaxation 

We  are  particularly  interested  in  electron-hole 
euergx  relaxation  via  phonon  emission  and  absorp- 
tion, as  it  relates  to  the  electron-hole  tempera- 
ture relaxation.  The  diagrams  in  Fig.  12  show 
phonon  emission  and  absorption.  As  discussed 
earlier,  electron  and  hole  transitions  are  assumed 
to  be  intraband  and  intravalley  in  phonon-assisted 
relaxation.  The  calculation  of  (®n/#f)*EL  from  the 
diagrams  of  Fig.  12  is  quite  similar  to  the  FCA 
calculation.  As  shown  in  Appendix  D,  we  find 


FIG.  12.  Phonon-aa- 

elated  rclaxat'on. 


FIG.  13.  Optlcal-phonon-aaststcd  relaxation  va 
electron-hole  temperature. 

(S)„L“[r°(0  + r',(01  ’ (<5) 

where 

r0«)«r,0'r(r-  rj 

x [ln(l  + e"*t)  + 5(»«c/m4)*  ln(l  + «"*>)] , (46a) 

rA(f)  = r"’T,(T  - r£)[x(ca)  ♦ 5(m,/mt)*«(<l)]  , (46b) 

r<0)  » 2QjmJui0/pir3N0)f 4 , (46c) 

r^'-ieA'mJ , (46d) 

and  k(x)  is  defined  in  Table  I.  ro(f)  is  the  relaxa- 
tion rate  due  to  the  optical  phonons,  rA(/),  due  to 
the  acoustical  phonons.  For  Ge,  r<0,»2  x 1010 
sec'1  and  rw,*5x  lO10  sec'1  at  300  CK.  Owing  to 
the  (linear)  dispersion  of  acoustical  phonons,  the 
temperature  dependence  of  ro(f)  and  rA(f)  is  dif- 
ferent. Plots  of  r0(l)  and  rA(f)  vs  (HujJkB)<r  (=T) 
are  given  in  Figs.  13  and  14.  For  large  tempera- 
tures (fixed  r) 


FIG.  14.  Acoustical- phonon- anaistrd  relaxation  va 
electron-hole  temperature. 
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llm  T0«) 


a,7T 


*9^,0'nrl,, , 


(47a) 


lim  T„(/)  - 


Aj~t 


*8T «•'n•r,/,  ; (4Tb) 

thus,  relaxation  via  acoustic- phonon  emission 
and  absorption  becomes  more  important  at  high 
temperatures.  In  this  limit,  ro  + r A“n.  Since 
lirnu  « it  as  f -«o,  the  rate  of  relaxation  ol  7 is  in- 
dependent of  the  density  n as  T — On  the  other 
hand,  when  T is  near  Tt  (but  »i  is  much  greater 
than  the  thermal  equilibrium  density), 

T0(f)-(3/2AJ*/1[l  + i(V'41)</3l 

xr'V’ir-  rt)  = 7r«°v/3(r-  rt) , (48a) 


rA(f)'i(3/2AJ)</1[l+(l/55/s)(4lMJ)</1| 


xr|,V/,(T-  rl)*6r<a'n,,'(r  - Tl)  ■ (48b) 


Thus  the  density  dependence  of  T0(f)  and  TA[t)  is 
different . Note  also  that  the  relaxation  of  7 is 
density  dependent  since  (3u/&T)„rx  itWl  in  this 
limit. 

The  specific  heat  of  the  lattice  is  approximately 
constant  for  Ge.”  Therefore 


duL  I 9u  \ <iT L 

dt  ="'97/«L  = rt  dt 


(49) 


where  cL*  CL/SQkB  and  CL  is  the  specific  heat  of 
the  Ge  lattice  in  erg/cmJ  ”K.”  Thus 


m—  [ro(0  + rM(0| . 

CL 


(50) 


to  the  tact  that  direct  absorption  populates  only 
the  T valley.  Therefore,  the  Fermi  energy  of  the 
T-valley  electrons  is  perturbed  relative  to  the  Fermi 
energy  of  the  L-X  valley  electrons  when  the  exci- 
tation pulse  is  on.  This  relative  perturbation  is 
rapidly  damped,  as  the  two  Fermi  energies  try  to 
rapidly  equalize  by  means  of  the  phonon- assisted 
intervalley  scatterings  between  the  T valley  and 
the  side  valleys.  It  is  Just  such  rapid  damping  that 
causes  the  substantial  broadening  of  the  plasmon 
resonance.  Note  that  Coulomb  collisions  (con- 
sidered in  Appendix  B)  are  ineffective  in  transfer- 
ring electrons  between  the  T valley  and  the  side 
valleys,  and  thus  would  not  damp  the  relative  per- 
turbation of  the  Fermi  energies  referred  to  above. 
Electrons  and  holes  can  of  course  relax  via  intra- 
valley and  intraband  phonon  emissions,  and  thus 
perturb  the  Fermi  energies.  But  these  are  slow 
processes  compared  to  the  r~L,X  scattering. 

Their  contribution  to  the  plasmon  resonance 
broadening  is  negligible. 

As  a result  of  the  broadening  of  the  plasmon  re- 
sonance, plasmons  have  short  lifetimes  and  plas- 
mon emission  will  essentially  be  spontaneous 
rather  than  stimulated.  An  electron  in  the  T valley 
will  recombine  with  a hole  near  the  top  of  the  val- 
ence bands  via  spontaneous  emission  of  a plasmon. 
Largely  long- wavelength  plasmons  will  be  emit- 
ted.” 

For  plasmons  which  have  long  wavelengths  and 
infinite  lifetime,  the  electron- plasmon  coupling 
(in  “random- phase  approximation”)  is  given  by  the 
Hamiltonian30 


»-f 


.t 


x<st|e,,*,|s'lc')i4^tc,.;.  + H.c. , 

(52a) 


The  energy-gap  parameters  Sc  and  S0  are  func- 
tions ol  7 L.  Therefore  relaxation,  in  addition  to 
retarding  the  growth  ol  7 , will  also  influence  n, 

«,  etc.,  through  the  variation  of  S0  and  S0  as  7t 
increases. 


where  A;  is  the  annihilation  operator  for  a plasmon 
of  momentum  5-  We  are  particularly  interested  in 
the  interband  transition  terms  (i.e.,  recombina- 
tion terms)  In  (52a).  Near  the  T-valley  band  edge, 
and  again  for  long- wavelength  plasmons  (3-0), 


E Plumon  >vmlrd  recombination 


(ct|e';-,|t-t)--[(3-Pc,)/m£0l61(C'-C  + 5),  (52b) 


The  combined  electron- hole  plasma  frequency  is 
given  by 


As  the  electron  and  hole  densities  build  up  u>,  in- 
creases in  magnitude,  iiecoming  comparable  to 
the  Ge  gap  frequencies.  The  plasmon  resonance 
in  substantially  broadened  when  the  excitation 
pulse  is  passing  through  the  sample.  This  is  due 


where  Pc**(rC|p|t,lt),  the  intraban  1 matrix  element 
of  the  momentum  operator  as  in  (25b).  From  W,, 
we  find  that  the  recombination  rate  for  a T-valley 
electron  with  momentum  k via  spontaneous  emis- 
sion of  a plasmon  (of  Infinite  lifetime)  is  given  by 


in£„ 


xo<fc(0-E,(e-3)-frw,). 


(53) 
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FIG.  IS.  Plssmoo-ss- 
ilitwj  recombination. 
Double  broken  line*  refer 
to  plasmon*. 


A,,  ■ H/2t0  represents  the  plasmon  resonance 
broadening  mentioned  earlier  [see  Eqs.  (56c)  and 
(56d)|.  We  set  (2/*)6„  equal  to  the  T-L.X 
scattering  rate  calculated  in  Appendix  A: 


However,  in  our  problem,  plasmons  hare  a 
short  lifetime,  and  therefore,  the  A function  in 
(53)  is  replaced  by  a Lorentzlan.  More  precisely, 
(53)  is  replaced  by11 


In>(«rt,  [£„(£)  - Ejf.  - ®)/N)  )’ 

<M) 


When  the  plasmon- assisted  recombinations  are 
appreciable,  electron  and  hole  densities  are 
sufficiently  large  that  uyr0>  1,  and  near  the  plasmon 
pole  (w~i<>,) 

€j(w) » I - w*/w*  *(2/ai,)(w-  <•»,),  (56a) 


where  *(q,  u>)»«,(5,  ui)  + «c,(q,  u>)  is  the  dielectric 
function.  Eq'^.tion  (54)  Is  obtained  from  the  self* 
energy  diagram  shown  In  Fig.  15.  We  will  approx- 
imate r „ by  the  rate  for  those  electrons  and  holes 
which  are  near  the  band  edge  at  T;  i.e.,  we  set 
*»0  and  E#0 -£.»£„♦  *V/2m».  We  will  also  as- 
sume a high-frequency  Drude  form  for  the  dielec- 
tric function: 


«,(w)~1/u»,t0, 


(56b) 


C(u>)  * 1 


U)(o>+  </t#)  ' 


(55)  From  these  approximations  and  (25b),  we  obtain 


2r*u).  f , , 6/2ra 

Integration  over  q gives 

r«  ’ '»♦£.-  iruyF'  • 


(56d) 

(57) 


Finally,  we  can  write 

(M) 

where 

«r’»',,,P»(*i). 

(59a) 

"r  m*l,: V»(*.> 

(59b) 

(see  Table  I). 


F Rile  equations  for  absorption  and  Iransmiaaion 

We  are  now  ready  to  derive  the  rate  equations  for 
describing  the  absorption  and  transmission  of  the 
excitation  pulse.  Let  us  consider  the  rate  equation 
for  N(l)  first.  It  can  be  written 


2£.-a(f)A /(/)♦*</)-/(/),  (60a) 

<*(/)- <*„(/)+ orrCA(f),  (60b) 

s(f)  is  the  normalized  photon  density  per  unit  vol- 
ume and  per  unit  time  that  enters  the  interaction 
region  (source)”  and  /(f)  is  the  normalized  photon 
density  per  unit  volume  and  time  that  leaves  the 
interaction  region.”  We  must  calculate  /(f)  in 
order  to  complete  the  equation  above;  however,  we 
would  first  like  to  make  a few  remarks  about  (60a). 

Let  the  excitation  pulse  enter  the  interaction 
region  at  f»0,  and  let  A/(/»0)» 0.  Then,  the  formal 
solution  of  (60a)  is 

N(l)  -J'dt'  [s(f')-/(/')J  exp(-  £ o(f')df').  (61) 
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Ut  ua  define 

#/)•  f'dftU').  (Mr) 

£(0-  f'dt'Ht').  (63b) 

S(0  represents  the  Integrated  photon  density  enter- 
ing the  Interaction  region  up  to  the  time  I;  L(t ) re- 
presents the  Integrated  photon  density  that  has  left 
the  region  up  to  the  time  <[S(0)  »0,£(0)  »0).  Thus 

M0«S(0-£(0 

-[jfV  ar(/')[s«') -£(/')] 

x exp^-  J a(f')rf/*j].  (62c) 

From  this  form  it  is  clear  that  the  last  term  on  the 
right-hand  side  represents  the  radiation  which  Is 
absorbed  up  to  time  t: 

MO  - jf  V a(/')[S(/')  - L(t')\  exp(-  f*  a(/')<tt'). 


(62d) 


Thus 

S{t)  = A(t)  * N(t)  * L(t) , 


(62e) 

which  Is  an  obvious  result  based  on  the  conserva- 
tion of  energy.  In  the  experiments  L«S,  and 
L«A*N~A.  Note  that  dA/dt « a(l)N(t). 

The  absorption  coefficient  a is  generally  within 
the  range  of  10M-10,,  sec"*.  Initially  we  have 
a(/-0)~10“  sec*1.  Clearly,  the  exponential  factor 
In  (61)  dominates  rapidly  [i.e.,  a(0M0»  Af(/)]t 
and  as  MO  Increases  from  its  initial  zero  value, 

It  will  quickly  reach  the  asymptotic  value  Implied 
by  (60a): 


MO  -[*(/)  -/(/))/<*(/) 


(63) 


In  a time  of  order  (1/a) -10*“  sec.  Since  the  pulse 
width  Is  much  larger  than  this  value,  the  Initial 


FIO.  IS.  Partitioning  of  optical  pulse  to  calculate 
trunamlealon 


rapid  transient  behavior  of  MO  can  be  neglected, 
and  MO  can  be  taken  as  given  by  (63). 

la  order  to  calculate  /(/),  let  us  consider  a thin 
Slice  of  the  pulse  between  I and  t*6t  (see  Fig.  16) 
that  has  entered  the  crystal.  The  number  of  pho- 
tons In  that  slice  is  given  by  s6t.  The  passage  time 
of  the  slice  is  /,  = (*.)*/a£/c,  where  £ is  the  width 
of  the  Interaction  region  (i.e.,  the  crystal).  f, 

»7«  10*“  sec  for  £*5.5  pm;  thus  f,  is  much 
smaller  than  the  pulse  width  T.  During  the  pass- 
age of  the  "slice  of  light”  through  the  crystal,  a(0 
will  not  change  much  since  /,«  T , and  a(0  can  be 
considered  to  a good  approximation  as  a function 
of  A(t),  and 


A(t*tJ-A(t)  t , dA 

m *mD  dt  ■ 


(64a) 


Since  dA/dt  = aN,  and  from  (63)  (where  we  have 
noted  that  f«  s)  and  from  (62e)  A- S~ st 


A(t*t.)-A(t) 

Mi) 


(64b) 


for  f,«  t-T.  The  quantity  /,<*  >1  and  multiple  re- 
flections can  be  neglected.  For  this  thin  slice,  the 
usual  absorption  equation 


(64c) 


can  be  Integrated  from  x »0  to  x »£  and  then  mul- 
tiplied by  (1  - r)  to  account  for  the  reflection  from 
the  back  surface  to  obtain  the  density  of  transmit- 
ted quanta, 

6L  - (1  - r)6S(x  -£)  - (1  - r)*6te~,»m,n , (64d) 

where  r is  the  reflection  coefficient,  thus 

l(f)»(t-r)s(f)e-V‘,>  (65) 

and 

.« 


1(0  - j dt'  (l  - »•)*(/')*-'»•“'». 


(66) 


The  transmission  coefficient  for  the  excitation 
pulse  la  given  by 

(Va -*-)»(>')*-•#•<'•>,  (67) 

4i  si  Jt 

where  SL  represents  the  normalized  density  of 
photons  in  the  laser  pulse  {Kw0N£l  is  the  ratio  of 
the  laser  pulse  energy  to  the  Interaction  volume.) 
Note  that  for  a rectangular  pulse  of  width  T 


s(0 


1“ 

( 0 


- r)St/T  for  Os  IsF, 
for  T<t, 


(68) 


where  (1  - r)  in  front  of  SL  accounts  for  the  reflec- 
tion from  the  front  surface  of  the  sample,  and  thus 
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The  reflection  coefficient  was  experimentally  mea- 
sured for  all  intensities  of  the  laser  beam  used  in 
the  experiment,  and  was  found  to  be  nearly  con- 
stant .**  It  is  given  by 

r = (V«I-l)V(V«:  + l)*-0.J6.  (70) 

In  order  to  collect  together  our  working  equa- 
tions, we  substitute  the  results  given  by  equations 
(26),  (SB),  and  (63)  in  (10)  and  obtain 

>.(f ) ^ ^ (s  - /)  - (nr  - »r>r..  (71a) 

Likewise  we  substitute  (27),  (40),  (45),  and  (63)  in 
(20)  in  order  to  obtain 

* ^ (S - Z)  ♦ ^ (s-0- [ro(0 ♦ r „(/)).  (71b) 

Solving  (71a)  and  (71b)  for  df./dt  and  dT/dt,  taking 
dh/dt  from  (13),  dTjdt  from  (50)  and  /(/)  from 
(65)  we  have  our  complete  set  of  working  equations: 


n-(1, fK-'w. 

-r.-rj, 

(72a) 

♦(5)(«-*-r  0-r4>. 

(72b) 

dh  d*  dr 

di’-y+di*‘r'dt  ’ 

(72c) 

(72d) 

I(f)-(l-r)s(f)«*'»*<*', 

(72e) 

where 

■{r*/4[p,  (*,)-<>,(*,)]} 

x{l0[vlu,)+  v,UJ)|[p,(£, )-<;,(*,) )[p,(*,)  - 0,(1,)] 

- 9[P,U.)  r[2p,U,)  - o,U.)  - »,(*,)]}• 

(72/ ) 

Given  s(f)  and  «(0),  *(0),  r( 0),  rt(0),  these  equa- 
tions can  be  solved  numerically  for  <((),  h(l),  T{t), 
rL(t),  and  1(f)  (or  L). 
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G.  Integration  of  dynamical  equations 

We  now  make  a few  qualitative  remarks  about 
(72a)- (72e)  and  present  their  numerical  solutions. 
At  the  initial  phase  of  the  excitation  pulse  absorp- 
tion, T increases  quite  rapidly  compared  to  n 
since  after  a direct  optical  transition  an  electron 
in  the  T valley  scatters  to  one  of  the  L-X  valleys. 
After  such  a scattering  event  an  energy  of  approx- 
imately 1-  Sc  goes  into  electronic  thermal  agita- 
tion. Therefore,  r approaches”  i(l  - So)(~1500'X) 
rapidly  and  c decreases  lo  a negative  value.  In 
this  initial  phase,  plasmon  recombination, 
phonon  relaxation  and  free-carrier  absorption 
terms  are  completely  neglibible;  thus,  we  may 
apply  a nondegenerate  (high  temperature)  approx- 
imation to  the  integrals  in  the  y coefficients  (see 
Table  II).  After  this  initial  state,  T may  decrease 
or  increase  depending  upon  the  comparative  rates 
of  the  phonon  relaxation,  plasmon  recombination, 
and  the  free-carrier  absorption. 

As  n increases,  initially  the  transmission  also 
increases.  Eventually  if  n is  sufficiently  large, 
plasma  recombinations  begin  to  compete  with  di- 
rect absorptions  and  n begins  to  saturate.  The 
saturation  region  is  reached  when 

(aa/a)s-(nr-;ir)rV  (73a) 

This  corresponds  to  dn/dt- 0.  In  this  region, 
those  electrons  which  have  just  made  direct  op- 
tical transitions  can  recombine  with  holes  via 
emission  of  plasmons  and  become  available  again 
for  direct  absorption.  Plasmons  rapidly  decay  and 
transfer  their  energy  to  electrons  and  holes.  Thus, 
in  the  saturation  region,  the  main  result  of  direct 
absorption  is  the  heating  of  electrons  and  holes. 

TABLE  II.  High-temperature  approximation  flow 
densities). 

o,U)-(V?/4V7Uk'‘' 

A (S'* 

■ - - <V7/2u,r,/,«-« 

y,  **  n/r 

vr—  <*/r)<i  5 -«,)<•  -u*,/r) 

Y«a-> 

ym~<*/rHS0*sr) 

y Wf  - W/TW  (1  • : s ♦ a ,1 >(Sc  ♦ 3*' (2  » ‘ a , >1  • <*a  ,/rius c ♦ 3f1 

n-3sVr 
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FIG.  17.  Temporal  behavior  of  Fermi  energies  during 
pan  a age  of  excitation  pulse 


In  this  region,  T then  varies  according  to 

5-£(s- r.-r.)  (73b) 

and  increases  rather  rapidly.  In  the  saturation 
region,  the  trajectories  of  the  system  remain 
close  to  the  constant  n curves  in  Fig.  10. 

The  following  several  figures  present  the  re- 
sults of  numerically  integrating  (72a)-(72e)  during 
the  period  the  excitation  pulse  is  present  in  the 
sample.  All  numerical  solutions  have  been  gen- 
erated assuming  a rectangular  optical  pulse  of 
width  7.  Figures  17  and  18  Illustrate  the  temporal 
evolution  of  the  Fermi  energies  of  the  electrons 
and  holes  and  of  the  carrier  temperature  as  an 
excitation  pulse  of  constant  energy  traverses  the 
sample.  We  have  presented  these  curves  for  sev- 
eral different  pulse  widths  to  emphasize  that  the 
pulse  width  has  a significant  influence  on  the  temp- 
oral evolution  of  these  quantities.  Note  that  for 
the  large  excitation  pulse  energy  used  in  these  fig- 
ures the  temperature  begins  near  1500  “K  and 


FIO.  IS.  Temporal  behavior  of  electron-hole  tempera- 
ture. 


FIG.  19.  Electron  density  vs  excitation  pulse  density 
for  different  pulse  widths  and  Initial  lattice  tempera- 
ture of  77  *K. 

rises  due  to  plasmon  recombinations  as  discussed 

above. 

Figures  19-24  show  the  final  carrier  number, 
carrier  temperature,  and  lattice  temperature  im- 
mediately after  the  passage  ol  the  excitation  pulse 
as  a function  of  excitation  pulse  energy  (or  lattice 
temperatures  of  77  and  297  K.  Again  the  depen- 
dence of  these  curves  on  pulse  width  is  emphasized. 
Notice  that,  due  to  plasmon  emission,  the  elec- 
tron number  begins  to  saturate  somewhat  and  the 
carrier  temperature  is  elevated  for  large  excita- 
tion pulse  energies  as  discussed  earlier.  Also 
note  'hat  this  model  does  not  provide  for  sufficient 
lattice  heating  to  account  for  the  surface  damage 


FIG.  20.  Electron  density  ve  excitation  pulse  density 
for  different  pulse  widths  and  Initial  lattice  temperature 
of  297  *K. 
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Incidan'  Rodtalioo  Oanaiiy  (cm'*l  Inodant  Rodiotioe  Ddmdy  (cm**) 


FIG.  21.  Electron-hole  temperature  vs  excitation  FIG.  23.  Lattice  temperature  v*  excitation  pulse 

pulse  density  for  different  pulse  widths  and  Initial  lat-  density  for  different  pulse  widths  and  Initial  lattice 

tlce  temperature  of  77  *K.  temperature  of  77  *K. 


we  observe  at  the  largest  pulse  energies.  In  Figs. 
25  and  26,  we  show  the  corresponding  single  pulse, 
or  excitation  pulse,  transmission.  These  curves 
are  presented  as  a function  of  pulse  width  for  lat- 
tice temperatures  of  77  and  300  “K. 

For  t >T,  we  set  s <0,f  >0.  The  system  begins 
to  evolve  according  to 

£■(*)«'•• r-’’  <«*> 

<™> 

rIT L dt  and  dh/dt  are  still  given  by  (72c)  and  (72d). 
In  (74a)  and  (74b)  we  have  omitted  the  plasmon-as- 
sisted  recombinations  since,  as  discussed  previ- 
ously, when  the  excitation  pulse  is  turned  off  the 
plasmon  process  is  likewise  rapidly  turned  off. 
Notice  that  in  the  experlmer.t  the  largest  probe 
pulse  delay  was  about  300  psec  and  that  direct 
and  indirect  optical  recombinations  are  still  neglig- 


ible for  such  delay  times.  Since  the  experimental 
probe  pulse  intensities  were  small  compared  to 
the  excitation  pulse  intensities,  the  system  evolves 
according  to  (74a)  and  (74b)  and  (72c)  and  (72d)  when 
the  probe  pulse  passes  through  the  interaction  region, 
rather  than  (72a)  and  (72b)  with  s = , l * lfnt, . 

For  a rectangular  pulse  of  width  T,  the  probe 
transmission  is  given  by 

L r+.  * / df  (1  - r)spt*.(f')e-'»«<«  (75) 

and  the  transmission  coefficient  by 

Tf**  fV,/,rff' (76) 

where  tt  is  the  delay  time  of  the  probe  with  respect 
to  the  excitation  pulse.  For  tt>l  , Tf0*  Increases 
at  first  due  to  electron- hole  relaxation.  As  elec- 
trons and  holes  relax  via  phonon  emissions,  they 
clog  the  optically  coupled  states  E'  and  E"  in  Fig. 


Incidtni  Rodiouofl  Dtnvty  Ictn'*) 


FIG.  22  Electron-hole  temperature  vs  excitation 
pillar  density  for  different  pulse  widths  and  Initial  lat- 
tice temperature  of  297  *K. 


tnodant  RodiolKm  Oanaily  (cm'*) 

FIG.  24.  Lattice  temperature  vs  excitation  pulse 
density  for  different  pulse  widths  and  Initial  lattlca 
temperature  of  297  *K. 
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FIG.  25.  Excitation  pulse  transmission  vs  excitation 
pulse  density  for  Initial  lattice  temperature  of  77  *K. 


5 and  reduce  Absorption.  Further  reiAXAtion  can 
eventually  free  these  states  (or  optical  transitions, 
and  probe  transmission  can  decrease.  These  ideas 
are  summarized  in  Fig.  7.  Figures  27  and  28  show 
probe  transmission  versus  delay  time  (or  different 
excitation  pulse  radiation  densities  and  (or  the 
two  initial  sample  temperatures  77  and  297  JK. 

For  these  solutions,  the  pulse  width  is  assumed  to 
be  5 psec.  As  the  excitation  pulse  intensity  de- 
creases, the  probe  transmission  curve  is  general- 
ly lowered.** 

Ill  COMPARISON  WITH  EXPERIMENTAL  DATA 
AND  CONCLUSIONS 

The  experimental  data  and  'he  theoretical  curves 
for  the  excitation  pulse  transmission  are  shown  in 
Fig  29.  The  theoretical  curves  are  for  a 5-psec 
pulse  width  and  a 3 x 10-,-cmJ  interaction  volume. 
This  volume  corresponds  to  a focused  spot  diam- 
eter that  is  in  accord  with  the  experimentally  es- 
timated value  of  about  250  pm.  The  agreement 


FIO.  26.  Excitation  pulse  transmission  va  excitation 
pula*  density  for  Initial  temperature  of  297 *K. 


FIO.  27.  Probe  transmission  vs  delay  time  for  differ- 
ent excitation  pulse  densities  and  initial  lattice  tempera- 
ture of  77  *K. 

between  data  and  theory  is  good  at  lower  excita- 
tion energies.  At  the  extremely  high  excitation 
energies  (-4  x 10M  quanta)  the  theoretical  enhanced 
transmission  does  not  saturate  as  does  the  experi- 
mental transmission.  However,  since  we  arc  un- 
certain whether  the  saturation  in  experimental 
transmission  is  due  to  bulk  effects  or  due  to  the 
onset  of  surface  damage,  we  shall  not  dwell  on 
this  point. 

Figure  30  shows  the  probe  transmission  data 
(minus  the  zero  delay  time  spike)  and  the  corre- 
sponding theoretical  curves.  The  theoretical 
curves  were  generated  for  2 x 10JO-cin‘'  incident 
excitation  quanta  density,  5-psec  pulse  width, 
and  lO^-cm5  interaction  volume.  The  agreement 
between  the  theory  and  the  experimental  data  ts 
reasonable  lor  both  77  and  297  "K. 

The  overall  good  agreement  between  the  theory 
and  the  experiment  permits  us  to  draw  a number 
of  conclusions.  An  important  assumption  in  our 
calculations  was  that  the  large  number  of  free 
carriers  created  by  the  excitation  pulses  be  de- 
scribed by  Fermi- like  distributions.  Apparently, 
this  assumption  works  quite  well  despite  the  in- 


FfG.  28.  Probe  transmission  vs  delay  time  for  dif- 
ferent excitation  pulse  densities  and  initial  lattice  tem- 
perature of  297 'K. 
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FIG.  29.  Theoretical  and  experimental  excitation 
pulse  transmission  vs  excitation  pulse  intensity. 

tensity  and  the  short  duration  o(  the  optical  pulses. 
Clearly,  in  order  to  observe  a large  deviation  from 
a Fermi- like  distribution  in  a thin  sample,  one 
must  go  to  even  shorter  pulse  widths  than  those 
achieved  in  these  experiments.  It  also  appears 
that  if  the  sample  is  only  a few  micrometers  thick, 
then  one  can  average  over  the  spatial  effects  and 
obtain  good  results. 

Many-body  (plasmon)  effects  are  extremely  im- 
portant in  intense  nonlinear  light  absorption  owing 
to  the  large  number  of  carriers  created  and  the 
short  lifetime  of  the  optically  coupled  states.  Plas- 
■non-assisted  recombinations  play  a role  in  the 
behavior  of  both  excitation  and  probe  transmission 
by  saturating  the  free  carrier  densities  while  in- 
creasing their  temperature.  Plasmon- assisted 
recombinations  should  have  observable  conse- 
quences in  similar  experiments  with  other  mater- 
ials. How  they  would  affect  the  excitation  and 
probe  transmission  would  depend  on  the  particul- 
ars of  the  energy-band  structure  of  each  semi- 
conductor. 

Let  us  emphasize  again  that  the  highest  plasma 


frequencies  achieved  in  the  experiments  are  actual- 
ly less  than  EJ 2tl  The  maximum  electron  den- 
sities occurring  in  the  vicinity  of  2 x 10”  cm'1, 
which  corresponds  to  (maxw,)/^  =»0.3.  There  is 
a directly  observable  quantity,  namely  the  refrac- 
tive index  for  the  optical  frequency  u>0,  which  is 
directly  related  to  the  ratio  up/w0.  For  the  elec- 
tron density  quoted  above,  the  maximum  change 
in  the  refractive  index  divided  by  the  initial  re- 
fractive index  is  given  by37 


max 


( ^r»t  lain 
\ref  index 


-0.05. 


This  is  in  good  agreement  with  the  direct  measure- 
ments of  the  refractive  index.4 

In  Sec.  II,  the  role  of  the  phonon-assisted  re- 
laxation was  discussed  in  detail.  We  have  seen 
that  the  reason  the  relaxation  extends  over  a time 
span  of  a hundred  picoseconds  is  due  to  the  fact 
that  the  phonon- assisted  relaxation  in  Ge  is  strong- 
ly dependent  on  the  difference  between  the  temper- 
atures of  the  free  carriers  and  the  lattice  and  due 
to  the  existence  of  conduction  band  valleys  in  Ge. 

In  different  materials  with  different  energy-band 
structures  and  electron-phonon  coupling  constants, 
the  behavior  of  the  probe  pulse  transmission  would 
be  (and  is)  drastically  different. 

As  mentioned  earlier,  the  lattice  temperature 
increase  is  generally  small.  Here,  "small”  can 
mean  as  much  as  10-30 'K  at  large  light  intensi- 
ties. This  may  be  sufficiently  large  to  be  experi- 
mentally measurable. 

An  important  assumption  in  the  calculation  of  the 
phonon- assisted  relaxation  rates  was  the  fact  that 
phonon  distributions  were  given  by  the  usual  Bose 
distributions.  This  assumption  apparently  works 
well.  However,  at  large  intensities  (near  1015 
laser  quanta),  large  numbers  of  phonons  are  gen- 
erated. It  is  quite  possible  that  a phonon  instability 
may  occur,  i.e.,  the  number  of  quanta  in  some 
phonon  modes  may  grow  to  larger  values  than 
those  given  by  equilibrium  values.  This  should  be 
investigated  in  the  future. 

An  interesting  feature  of  the  theoretical  model 
as  developed  is  the  behavior  of  the  transmission 
for  different  pulse  widths.  One  can  see  from 
Figs.  25  and  26  that  transmission  changes  signi- 
ficantly near  the  saturation  region  when  pulse 
width  is  increased.  Sir  e it  is  possible  to  split  an 
optical  pulse  into  separate  parts  and  then  rejoin 
them,  it  is  possible  to  perform  the  experiments 
with  variable  duration  pulses.  These  and  other 
experiments  are  currently  under  way. 
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appendix  a r - L.X  transitions  via 

PHONON  EMISSIONS 

Let  us  consider  a particular  side  valley,  say 
(lll)valley.  Usingthe  electron-phonon  coupling 
Hamiltonian  given  by  (35),  the  transition  rate  for 
an  electron  going  from  the  T valley  to  the  (111) 
valley  via  spontaneous  emission  of  optical  and 
acoustical  phonons  can  be  written 

* |Wu  W)|  *Kc,k  - q’  |e",**‘,|c ,k>|* 
X6(£e(k)-£c(it-q') (Ala) 

where  it  is  sufficiently  small  such  that  )c,k>  re- 
fers to  an  electron  in  the  T valley  and  -5'  is  suf- 
ficiently large  and  generally  pointing  In  [ill] 
direction  such  that  |c,k-q')  refers  to  an  electron 
In  the  [ill]  valley.  If  we  let  ktlll)  represent  the 
wave  vector  of  the  lowest  point  (in  energy)  of  the 
[ill]  valley,  thenq'  and q'  »k.  Since  is 

smaller  than  the  reciprocal -lattice  vector,  |(r,k 
-q'|e',J,'r|c,k)|8*  1.  Electron-phonon  coupling 
parameters  W„(q')  were  discussed  in  the  text. 

For  optical  phonons  VU)»(<i0cosif>  and 
»V0sln<J>,  where  <f>  is  the  angle  between  -q'  and 
k,lu).  Since  -q'  is  assumed  to  be  generally  in  the 
direction  of  kUu(,  Q w *y0  and  (jfo  * 0.  Thus  for 
optical  phonons  the  right-hand  side  of  (Ala)  be- 
comes 

£ M£c(k)-£c(k-q')-frn0).  (Alb) 

P*‘»  r 

For  acoustical  phonons  (A4  ♦ A,  cos’0) 

»q'(A4*  A.)  and  yTA»  -q'A,  cosO  sind>  ■ 0.  Since 
« **tiuj.WL*%<'(mi(A<+ A„)»yx.  Also,  cAq' 
"cA*Ciui“no.“  Therefore,  we  obtain  the  contri- 
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button  of  both  optical  and  acoustical  phonons 

*<*)r^r„.i 

£ Mi.dt-E.S-q')-*!!,).  (AlO 

P**0  f' 

In  our  approximate  band  model,  the  electron 
energies  are  given  by 

S,(S)-C0*»V/2*«0  (Aid) 

and 

£,(k  - q) » £c  ♦ [fT*(k  - q'  - 5,ull)*]/2m,.  (Ale) 


The  phonon  energy  in  the  5 function  can  be  neglec  - 
ted,  since  it  Is  qul*e  small  compared  to  electron 
energies,  letting  -q  - k,,,,,*  q’  , where  g* 

<Aq',k tuil-  ,h*  >ot*8r®l  o\«r  the  * function  is 
trivial,  and  yields 


"Wr-*nu» 


/ IV 


(A2) 


For  an  electron  just  raised  to  the  F v.  llev  by 
optical  transition  £„♦  h‘k'  Thus,  sum 

mlng  over  all  valleys,  we  obtain 


1 . tr  * , via 


(A3) 


For  mc»0.22»i,  £*»„»  1.17  eV , £c»v.1  eV,  p 
* 5.3  g/cm1,  ftli0»  0.02  eV,  v0”  8 * l0"1  er«  cm 
and  V4  * 2 * 1 0*’  erg/cm,  1 A,  =»  1 * 10M  sec'1  " 'the 
absorption  rate). 


APTtNDIXB  COULOMB  THERMAL!/ ATION 

In  this  appendix  we  will  calculate  the  lifetime 
of  an  occupied  state,  which  is  finite  due  to  Cou- 
lomb scatterings  among  electrons  and  holes  if 
the  distributions  are  not  Fermi-like.  For  the 
sake  of  simplicity,  we  will  neglect  the  long-wave- 
length screening  of  the  Coulomb  fields.  For  an 
electron  state  k in  the  ith  valley,  the  scattering 
rate  out  of  the  state  k is  given  by 


* (S/.,<*Xl  -/„(£'  ♦ q)]6(£el(k  -q)  - £.,<£)♦  £.,(£'  ♦ q)  - £„(*» 

♦ £ ( 1 - ',.(*' )]/.,(*  ♦ q)«<£t((k  -$)  - £.,(£)  - £,.(*  ♦ i)*  £„(V')>). 


/ varies  over  the  conduction-band  valleys,  i over  the  two  valence  bands.  ftl  and  ft , are  electron  and  hole 
distributions,  respectively;  Etl  and  £,,  are  electron  and  hole  energies  (£,„  is  positive). 
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The  scattering  rate  Into  the  state  k la  given  by 


x 


(Ell  + q)M£c<(*  - <?)  - £„(*)  + £„<k'  * q)  - £t,(k')> 

-El1  -/a.S' + q))/»,(*')M*.,(k  - 5)  - £.<(i)  - £,.(*'  - q)-  £».(*)))• 

v 


/« 


(B2) 


Of  course,  if  the  distributions  are  Fermi-like, 
with  a common  temperature  for  electrons  and 
holes,  then  r|*«  f^"’,  for  any  state  k.  We  will 
evaluate  T’s  for  two  kinds  of  rather  simple  dis- 
tributions which  are  energy  dependent  only.  To 
simplify  the  algebra,  we  will  take  the  effective 
masses  of  electrons  and  holes  to  be  the  same  m* 
(since  nij*  0.22m , * 0.34»i , this  is  not  too  dras- 

tic). Also,  the  total  electron  density  Is  always 
taken  to  be  equal  to  the  total  hole  density. 

(a)  l.Ct  the  distributions  be  Lorentzian  with 
narrow  widths,  such  that 


!1  if  ir*(/C  -|A)V2m*<£,(k) 
<*J(A+iA)*/2w\  (B6b> 

0 elsewhere; 

and  2.5A<  A < (1.5+  v^A  (i.e.,  the  width  A is 
large). 

Let  us  consider  an  occupied  electron  state  k in 
the  rth  valley.  The  contribution  of  the  conduction- 
band  electrons  to  the  PJ"*  is  given  by 

r;»‘<«>.(5eVi*A»,«iV(*,K.),  (B7a) 


/„(*)-  *>/{[£c/(k)  - £c  - £]>  + A*} 
-»*  A6(£c/(k)  - £c  - £), 


(B3a) 


where 


(B7b) 


/„<*)•  (5A)Vil£.v(k)-£]**(5A)*} 

» 5*  A6(£„(k)  - £),  (B3b) 

i.e.,  6 • A/£«  1.  ^ 

(i)  Let  us  consider  an  occupied  state  k:  £,,(*) 
«£+£0.  Then, 

rj*-0,  (B4a) 

rj,,-(l60e4m*/3A,ei,)(l-6)>/,/6.  (B4b) 

(li)  For  the  states  k'  such  that  £cl(k')<  £♦  £c 
- A or  £„(£')>£♦  £c+  A, 

rJf-feOcNn’/W’ciHl  -6),n/b,  (B5a) 


fk*K.\  „ A.  . 2A'.(t-<+Kn(fe,+  i£!) 

/(*.£.)- 4 !„(_)_  2T+  — 


8 kK\  4 K]  - 8 k*K. 


(B7c) 


Let  *«A.+  a,  A s 3-0, 


/(*,  A.) 31  4 ln(4A'./6)+  J-  (A./6)s  + 0(l/6J)  - «.  (B8) 

For  fe»A+  I A, 

16  if  A*  2.5A,  (B9a) 

/ <*  26  if  A • 2.9A.  (B9b) 


rf'-O.  (B5b) 

Setting  ni  * "»  0.22m , «„  • 16,  we  find  80e4m  ‘/M’el 
» 101'  sec’1.  Thus  rj“,,ri?  "(light  absorption 
rate). 

(b)  Let  the  distributions  be  of  the  following 
form: 


Thus,  x lO11  sec'1.  The 

contribution  of  the  holes  is  about  the  same;  thus 
rj"^'  6 x ioM  sec*1  which  is  greater  than  the  ab- 
sorption rate.  Finally,  we  note  that  when  we  in- 
clude further  screening  effects  this  rate  drops  to 
around  10“  sec-1. 


!1  If  /T*(A  - 1 A)V2w  * < £,(£)  - E0 
<»>(A+lA)*/2m\  (B6a) 

0 elsewhere; 


APPENDIX  C RADIATION  ABSORPTION  RATES 

(a)  Direct  absorptions.  From  the  diagrams 
shown  in  Fig.  8,  we  readily  find 
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•*ye 


- -2 


800 

KLl<cr.t*q|e‘*-»«,(q)-p|i;,it>|1  [l -/*(k)J[l -/cr(k*5)]M£er(k  + 5)+ £«(k) -»w.) 

1)SKt'.kk'<,-,l(q)-pkr,k+^|*/er(k+i)/),(k)5<£cr(k+5) 

♦ £*(k)-Ku>,). 

(Cl) 

Here,  i4(q)  denotes  photon  polarization.  The  meaning  of  the  other  symbols  are  as  defined  in  the  text.  The 
sum  Ej  also  includes  spin  summation.  The  extra  factor  of  2 in  front  of  the  terms  is  due  to  the  fact  that 
we  have  two  degenerate  valence  bands.  The  photon  momenta  are  completely  negligible  compared  to  the 
electron  and  hole  momenta.  Therefore,  q's  in  the  interband  matrix  element  of  p,  in /' s and  in  the  fi  func- 
tions can  be  omitted.  Also,  the  radiation  is  primarily  at  the  circular  frequency  u>0  and  the  corresponding 
Nirta>  » 1.  Therefore,  the  spontaneous  emission  term  can  be  neglected.  To  obtain 

M[|  . 1 I 


it 


'o  1) 


it 


we  perform  Ekj  on  (Cl),  and  find 

| da  * -(j^— )^(0  / l<cr, S|p|  1 -/cr(k)  -/*(k)]a<£er(k)+  £»  -*u.0). 

The  Interband  matrix  element  of  p near  the  T-valley  band  edge  is  given  by 

|<cr,k|p|t),k>|’a,m,£0(m0+  mt)/2m^mk. 

in  our  approximate  band  structure,  the  electron  and  hole  energies  are  given  by 

£cr(k)»  £0*Htk’/2m0,  £«(£)■ 

The  integral  over  the  5 function  is  trivial,  and  yields 

9JL I . su)(  4<*f°  ) (2-n‘C * 

»t  |DA  ' \3tjr*u0/\  H\ni0+n;k)  ) 

* (l  -jw«p[*  - *•>♦*.)]}' ‘ - ji  ♦ (Eo*  - 


(C2) 

(C3) 

(C4) 


(Ku>0  -£.)_£ 


>-)!!  ")■ 


(C5) 


r 


^jr|  - -\(f)ajl  -£(z4)  -£(z,)] 

- -S(t)aD(t).  (28) 

To  find 

’"I  1 f 2d»k  9A,r(k)  I »/,r(k)  I 

w|»»  So  j (SP  at  I da ’ 9t  Id* 

can  be  calculated  in  completely  analogous  manner. 
It  is  given  by  the  negative  of  the  righthand  side  of 
(Cl)  If  Ef  is  replaced  by  E,t.  Thus,  in/it  |D»  is 
given  by  the  negative  of  the  right-hand  side  of 
(C2),  and 

§|m-NW-J1  -H*,)-Hx,)) 

•At(/)«a(f).  (C«) 


To  find  in/it\DA  we  need  ifH(H)/it  | DA: 


(16a) 

Because  of  the  two  degenerate  valence  bands 

WSI  . (C7) 

* DA  “ • DA 

Thu*, 


s?L*  / wlr-r<e>*  *■<»>  “Sr  L 
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Because  of  the  8 function  6(£cr(it)  + £„(£)-  8 j0)  in 
8'tr(k)/8f  |DA,  £cr(k)  ♦ E*(k)  = *<*>„!  »nd 


1 (2cP k 8/.r(k) 

Vl  DA^y0J(2i7  »/ 


I » N(t  )aB(t) . 
Ida 

(27a) 


(b)  Free-carrier  absorptions.  Let  us  consider 

J 


the  conduction  band  electrons  first.  Since  the  free- 
carrier  absorptions  strongly  depend  on  density, 
and  since  the  density  of  electrons  in  the  T valley 
is  negligible  compared  to  the  L-X  valleys,  we  Will 
consider  the  L-X  valleys  only,  and  ignore  the  T 
valley. 

With  the  help  of  the  diagrams,  like  the  ones 
shown  in  Fig.  9,  we  find 


aA’^1 

81 


2 n'1 

re*  -tt  ' * 


X/C<#)U  -/./(k  + q-  q')]MEc/k)  + rtwt-  £e/(£*q-  o’)  - ***.,•> 

..  | q'  I e*1*'**  I c? , it  + q)(c;  ,k  + ql  e***fCt(q)  *pl  cj , £) 

£ei(i«)  + '!<*>,-  + q) 

Qu(<r)(cj.£-q'l<,~<<''>lcj'*>fo'k  + q-q'lcl<'>^(q)'Plcj-k~q^  1 

£c,(k  + 5-T)-£el(k-T)-*w, 


( n \ 

( e ) 

\2pn„J 

\mc  / 

x/e/(k)[l  -/e/E  + q + q')]6(£Cj(k)  + Hu>t  - ECJ( k + q ♦ q')  ♦ .) 

x | Vu(q,Kc;,it-*-q'*-q'|gl<'‘>lc;>^*q><C;,k'tqlf<<'f4>(q)~Plcj.K> 

' EeJ(k)+  flu>t  - Ee/(k  + q) 

+ Qjq,Kcj.ic^qMel<-flcj,k)(c;,i<^qtq,le,i-fC(q)-p'cj.it*q,>  | 

q')  -*<•». 

-/*<(£-  q-  q')]8(Ecy(C)  - »u>,-  wn„,,-  Ee<(it-  q-q')) 

I 9.■(q,K^j,k-q-q,'g•‘t••^cj,k-q)<cf,g-qi<,•<<•fC(q)•glc>■g-q,) 
Kl  Et/(S)-E,#(E-®-fc-t 

i«>Jq,Kcj.g-q,l«’'‘<"f|0.gKcj.g-q-q,|<,',i'f^<q)'Plc^g-^,)  I 
Ew(K-q-q')«-*<*>f-  £.>(*- q'> 

x /,/£)(  1 -/«<(^+q'  - q)l®(£Ci^)+ Eci(E*q'-4)) 

, I V- (qO(c>,k.q>-$l^"^c;,k-q)(c;,i;-ql^-^(q)-glcj,g) 

I E#i(E)-E„(S-5)-*u>t 

, V . (q'Xc)  ■ * ♦ q' 1 e‘« ' ,f I cj , E)(cj , E ♦ q'  - ql  f'* ’ r(k(j)  • ff  I cj , C ♦ qQ  I* 
£,,(&  + q'  - q)  - £C,(E  ♦ q')  ♦ 

• mQ\ 


As  before,  j refers  to  the  valleys,  p to  the  phonon  modes.  n„,.  and  ijBf,  are  phonon  frequencies  and  phonon 
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dlstr  I Indians,  resisTtlvely.  (q')  arc  the  phonon- electron  coupling  parameters  discussed  In  the  text.  We 
will  annum*  that  there  are  no  phonon  Instabilities,  and  that 


(C») 


As  before,  the  light  momenta  q '<  k.q’  and  can  be  neglected  in  /'s,  A functions,  and  the  matrix  elements. 

In  the  >th  valley,  for  the  Bloch  states,® 

<e*,£*q|e,«'*tk(q)-p|c;,£>“(c*p£U1(q)-J  |cj,f>  »*(w /me) f4(q)  • K . (CIO) 

Since  q'  is  small  compared  to  the  reciprocal- lattice  rector® 

<c7,Ai<;'|e*,,-,|c;,A>*l.  (Cll) 

Phonon  energies  h n„  ,<<Ec.  h'jL>t  and  can  be  neglected  in  the  A functions.  We  will  also  neglect  spontane- 
ous photon  emission  terms,  as  done  before  in  direct  absorptions.  With  these  approximations,  and  after 
a redefinition  of  integration  variables,  we  obtain 


8 

8/ 


I FCa-EL.  V / A * I*  n**'  ' '*“»«•  ' 


X l/<i(k)-/«/iJ*q')]A(fc/(5)  + Ao;,-.£CJ(t  + q')).  (C12) 


We  again  average  over  light  polarization,  then  perform  (l/iV0)  Eu,  and  find 


8tf 

at 


- "(ssefe)  ??$  (*•*«£) 


(Cl  3) 


In  our  approximate  band  structure, 


Substituting  this  in  (C13)  and  summing  over  the  valleys,  changing  the  integration  variable  k,  we  find 
8.V  | „/  2(VV»  [ /VW/2mc-£,'|  „(Ea+  frV/2«i„  + Aw.- E, \] 

8<  IrcA-tt-1"'  V3mJ«.w*p ) J (2?T  l V k,T  / \ k.T  /J 

where  F(x)*  (1  ♦ r')*‘.  For  optical  phonons, 

Owing  to  the  A function,  integration  over  q*  is  quite  simple: 

A ^(kiq')^  2„m’  ir/^*2  . \l/*  r /*•**  . \,/a  f#V\x,r\*\ 

J -(ssr)  J| 


(CIS) 


For  acoustical  phonons,  n„f.  * c Aq'  and 
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One  of  the  Integrations  Is  the  same  as  In  (Cl 3).  The  other  gives 

**  wft  Z.  • *"•  - ’’’fa,?'1') 

Substituting  (C 15)  and  (C16)  in  (C14),  performing  the  remaining  angular  integrations  (of  k),  then  letting 
y -KJkr‘/2K)0mc,  we  obtain 

ir  ('♦  W)  ‘(wqS*-) Jf ' * w")‘ 

* f~dy{Hy  + i)l,i  + {y)'''T- !(>♦  D*/> -(>)*/aP}  [r  (y  * * ■-/*-)]  j . (C17) 

Thus, 

Sr  |,c*.«L="vta,x5(^)  (l+2^  r0to(£,)+a,rtX5(^)  toUl)+°,x5fe)  £lU,)] ' <CI8) 

a-V/a/|TCA  HOLES  can  be  obtained  from  (C18)  by  simply  omitting  (m^/nt,)  factors,  replacing  z,  - (Sc-  t)/T 
by  z1  = h/T , and  dividing  by  5.  Adding  electron  and  hole  contributions  we  obtain  8N/al|,CA  8>ven  by  (3*)> 
8u/af  |rCA  can  be  obtained  by  a simple  energy-conservation  argument: 


du 


1__  8£|  -1  3(*y/0AQ 


I FC A ^0^0  I fCA  ^W(A( 


a/ 


;arCA(OMl)- 


(Cl  9) 


re* 


APPENDIX  D PHONON-ASSISTED  RELAXATION 

Phonon- assisted  relaxation  may  be  represented  by  the  diagrams  shown  in  Fig.  12.  Let  us  consider  the 
L-X  valley  electrons  first.  Using  the  electron-phonon  coupling  Hamiltonian  given  by  (35),  we  find 


:>/cA> 
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Bit. 


2ir  IQ.(q,)l*>l 
n VV  2pn„. 

* [(•>„, . ♦ 1 H | (CJ , c I ■ ’ ■ ■ ! » I a , : E ♦ q'>  | yc,(C  ♦ q' ) [ 1 - /c,(k)  )B (FC,(E  ♦ q' ) - Fe/C)  - , .) 

- |W.C-q'|p-','-,|^.K)|Vci(S)U-/^-q'))5(£«/(C)-fc/(K-q')-*n...)} 

♦ V.  J | tt , 1 |e*« ' •' » | c) ,| k - qO  | •/.,£  - q' )|  1 - /C,(K) )6(£c>(C)  - F c< (E  - q' ) - >10. ,.) 

-|W,t*q'U',’"I^.K>,/e/K)U-/«<(C+q'))fl(^A*q’)-Ee/C)-fin.,.)}].  (Di) 


Using  (CU),  we  find 
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- l/.y«>  -/.,<* - ?)]6(£cy(E)  - f„(C  - ?)  - *<W>} 
♦/ey(£  * $'>(1  - fcfa  MB'fi- ♦ - *«<*>  - *W> 

-/,y(£)[i  -/.,<£)  - qMlME.itf)  - - q')  - 

We  multiply  (D2)  by  Eey(k)  and  sum  over  the  electron  states  and  the  valleys: 


*_YX  -/.!«>]  6lEey(£+<T)-  *.,<£>  - «W> 

p4V0/f<a;0  y(  M|»  *•*»*• 


r/  /C\  s /C 


♦£,/(*)/«,(*  ♦*')[»  -/.,(£))»(£«*  ♦4r>-  E./c>-  *n».‘> 
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(D2) 


(D3) 


Using  the  Identities 

/,/(k*q')[l  - /„(£)  j6(EeJ(k  ♦ q')  - Eey(£)  - AO, ,.) 

. + i)-l][,  _ (i  * e»'«e/,t>*«P,)-,)6(Ee/(£  + q')  - Ec/(£)  - 

»(e*,*0*«'>-  1)'*[(«*,,«<,I,'*,,+  l)"-(e*l*«7<l'*,‘0o-«r’4.1)-,]fi(£<iy(k  + q')-£<:y(k)-»n^) 
. («*«•"..*•  _ l)*‘l/cy(£)  -/cy(£ ♦ q')6(Ecy(k  ♦ q')  - E.y(£)  - >1^)  , 

we  find 


SL’kfcr)  ? 


(DA) 


last  ,j  ,» 

x <E  ,/*)[/,/*  > <7 ')  -/cy(*)  ] «(E  cy(*  + 9 ')  - E ty(fe)  - «1  ) 

- £«y(*)[/«y(*)  ~/.y<*  - «')]  «<£„(*>  - EcM  - 90  “ *0.^  » • 

Since  phonon  energies  are  small, 

!}.,  - (e**nM'  - 1)-*  -*.rt/»in„.  - k'T/hn.,  . 

Letting  k-k  + 5'  *n  tbe  second  term,  we  obtain 

7?l  ../(ssk)  2: 5 

x «<£„(£♦  $0-E,y(£)- 

Owing  to  the  S function, 

RSL  0 0 jg  M< 

Redefining  k in  each  valley,  we  obtain 

•SrI  . imiLz>sni  y [di  fay 
•f  l„L  p/V0Rw0(2«)’  7 q 

x [/tE+E0  + »n^)-./(E+E6)|o(* *-»yco?g  ♦ , 


(D5) 


(D6) 


(D7) 


(D8) 


(D9) 
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where  cost*  - k‘t)'  and 

h *.V/2m,  , f(F.  ).(«•'**'  •♦I)*1  , 
Mime  phonon  energies  are  small, 


nt+Ec  + h -/(£  (a-(^~ ■) 

can  be  omitted  from  the  6 (unction.  Thus, 

? / '*/  ^ •«£)• 


Substituting  the  properly  ahgular  averaged  electron- phonon  coupling  coefficients  we  find 

»f  hil  AVoW0(2»F^  y V q 'Vo  9 «A  \ >«c  2mf  ) 

The  remaining  integrals  are  trivial,  and  we  find  (in  terms  of  the  normalized  quantities) 

dn. 


* -r,0,5(— 'V  r(T - r, ) in(i  te-'i)  - r<*’5(— )'r2(r  - r ,)«(«,) . 
»'  'bil  V'"./  V"1./ 


(DIO) 


(Dll) 


(D12) 


(D13) 


(D14) 


9u/H,»u  can  be  obtained  from  (DM)  by  replacing  rlbyr1,  dividing  by  5,  and  omitting  (mc/mt)  factors. 


APPI  M>IX  t tXPl  RIMENTAL  METHOD  AND  DETAILS 

(a)  Introduction.  This  appendix  contains  an  elab- 
oration of  the  experimental  method.  This  presen- 
tation of  the  experimental  apparatus  is  given  in  de- 
tail because  picosecond  pulse  techniques  are  not 

a standard  tool  of  solid-state  physicists,  and  some 
aspects  of  sample  preparation  may  be  of  interest 
to  the  "nonsolid  staters."  The  first  experimental 
section  contains  a description  of  the  mode- locked 
Nd glass  laser  system,  electro-optical  switch- 
out  apparatus,  and  the  pulse  splitter  and  delay  op- 
tics. The  second  section  presents  properties  of 
the  germanium  sample,  including  purity,  surface 
preparation,  and  thickness  measurement.  The 
third  section  discusses  the  detection  system,  and 
the  fourth  discusses  possible  sources  of  experi- 
mental error.  Figure  31  shows  the  entire  apparat- 
us in  block  diagram  form. 

(b)  I’ul'.e  prc/nirtihon  The  production  of  very 
short  optical  pulses  by  mode-locking  requires  a 
laser  system  with  an  active  medium  that  provides 
a large  gain  bandwidth  for  light  amplification. 

Nd“  glass,  with  a fluorescent  bandwidth  of  2.7 

> in1'  Hz,  has  the  potential  of  producing  pulses  as 
short  as  0.3  psec.  In  practice,  pulses  between  5 
and  15  psec  duration  are  obtained.  This  deviation 
from  the  theoretically  calculated  pulse  width  is 
caused  by  the  broadening  of  short  pulses  in  a dis- 
persive and  nonlinear  medium  and  by  spectral 
narrowing  due  to  preferential  amplification  of  fre- 
quency components  close  to  the  gain  maximum.  In 
fact,  pulse  widths  would  be  wider  than  those  ob- 


served if  It  were  not  for  the  pulse  shortening  effects 
of  the  nonlinear  saturable  absorber  that  is  used  to 
passively  mode  lock  the  laser.  For  a discussion  of 
the  basic  principles  of  mode-locking  and  the  pro- 
duction of  ultrashort  laser  pulses,  the  reader  is 
referred  to  von  der  Linde.1 

The  mode- locked  output  from  the  Nd. glass  laser 
system  employed  in  these  experiments  consists  of 
a series  of  approximately  100  pulses  separated  one 
from  the  next  by  approximately  8 nsec.  Each  pulse 
is  typically  5-10  psec  in  duration  and  has  an  ener- 
gy of  0.1  mj  at  a wavelength  of  1.06  pm.  The  laser 
produces  single  well  mode-locked  pulse  trains  on 
more  than  80%  of  its  firings.  A brief  description 
of  the  laser  hardware  follows. 

The  active  lasing  medium,  a Brewster-cut 
Owens-Illinois  ED-2S  Ndiglass  rod  approximately 
19  cm  long  and  6.4  mm  in  diameter,  is  housed  in 
a Korad  K1  laser  head.  The  glass  rod  is  pumped 
by  a helical  xenon  flash  lamp  with  an  arc  length  of 
26  in.  A high-voltage,  symmetrically  shaped  pulse 
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FfG.  31.  Block  diagram  of  the  experiment. 
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from  a 10-kV  power  supply  Is  used  to  Ionize  the 
xenon  gas.  The  high-voltage  pulse  la  0.5  msec  In 
duration,  and  It  typically  delivers  an  energy  of 
000  J at  4.5  kV  Both  flash  lamp  and  glass  rod  are 
cooled  by  circulating  water  maintained  at  room 
temperature.  The  power  supply  consists  of  a 60- 
pF  capacitor  bank  plus  pulse  shaping  electronics. 
The  flash  lamp  and  power  supply  are  fired  approxi- 
mately once  per  minute. 

The  laser  oscillator  cavity  Is  120  cm  long  with  a 
spherical  mirror  of  radluB  2.5  m with  a reflectivity 
of  99. 7%  at  1.06  pm  at  one  end  and  a flat  mirror 
with  a reflectivity  of  27%  at  the  output  end.  The  la- 
ser Is  mode- locked  by  placing  a 1- mm- thick  cell 
containing  at  10:1  solution  of  dichloroethane  and 
Eastman  Q- switch  solution  A9860  In  contact  with 
the  output  mirror.  The  output  mirror  forms  one 
face  of  the  dye  cell.  Thus  it  must  be  resistant  to 
the  solvent  dichloroethane  as  well  as  withstand  the 
hlgh-power  levels  of  the  optical  pulses.  The  mode- 
locking dye  is  replenished  after  each  laser  firing 
by  siphoning  fresh  dye  solution  from  a rese  vior 
with  a hypodermic  needle. 

The  single  most  critical  component  of  the  entire 
laser  system  is  the  glass  laser  rod.  Self-focusing 
of  the  pulse  propagating  in  the  glass  rod,  due  to 
the  nonlinear  refractive  Index  n,,  Is  a major  prob- 
lem In  achieving  reliable  mode- locking  in  this 
laser.  Using  Owens-Illinois  type  ED-2S  laser 
glass,  we  obtained  clean  pulse  trains  on  approxi- 
mately 80%  of  the  firings. 

In  order  to  perform  the  two  experiments,  asliyte 
pulse  must  be  selected  from  the  train  of  picosec- 
ond pulses.  Figure  32  shows  the  electro- optical 
shutter  used  to  perform  this  function.  The  pulse 
train  passes  through  a Pockets  cell  placed  between 
two  Gian  prisms  which  act  as  nonabsorbing  crossed 
polarizers.  If  no  voltage  is  applied  to  the  Pockets 
cell,  the  light  passes  through  with  no  change  In 
polarization  and  is  rejected  by  the  second  crossed 
polarizer.  When  the  half  wave  voltage  Is  applied 
U>  the  Pockels  cell,  the  light  polarization  is  rotated 
by  90*  and  thus  the  light  is  transmitted  by  the  second 
polartser.  Voltage  is  held  off  the  Pockels  cell  until 


FIG.  3d.  Pulse  splitter  and  delsy  optics. 


the  first  part  of  the  pulse  train,  which  Is  reflected 
by  the  second  Gian  prism , breaks  down  a spark  gap. 
Upon  breakdown,  a voltage  pulse  from  a charged 
transmission  line  Is  applied  to  the  Pockels  cell;  the 
length  of  the  charged  line  is  chosen  such  that  the 
duration  of  the  voltage  pulse  Is  equal  to  the  tempor- 
al pulse  separation  In  the  train  In  order  that  one 
and  only  one  pulse  is  switched  out  of  the  train. 

In  order  to  perform  the  two  pulse  or  probe  ex- 
periments, the  single  pulse  selected  from  the 
train  must  be  split  In  two  and  a variable  delay  must 
be  Introduced  between  these  two  pulses.  Figure 
33  Is  a detail  of  the  pulse  splitter  and  delay  optics. 
Although  a weak  probe  pulse  is  generally  desired 
for  the  probe  experiment,  the  pulse  is  split  by  a 
60-80  beam  splitter  for  ease  in  alignment  of  the 
taro  pulses  coincident  on  the  sample  surface.  The 
variable  optical  delay  is  Introduced  by  means  of  a 
translatable  Porro  prism.  The  plano-convex  lens, 
focal  length/” 20  cm,  serves  to  Increase  the  optic- 
al energy  density  of  the  pulse  and  to  angularly  sep- 
arate the  two  pulses  for  detection.  The  lens  is 
placed  after  the  pulse  splitter  and  delay  optics  in 
order  that  the  probe  pulse  spot  size  in  the  plane  of 
the  sample  remains  independent  of  the  prism  posi- 
tion. Alignment  of  the  optical  components  and 
translation  of  the  prism  are  critical  in  order  that 
the  two  (xilses  are  spatially  coincident  on  the  plane 
of  the  sample  and  remain  that  way  during  the  ex- 
periment. Alignment  of  the  two  pulses  Is  confirmed 
by  Inspection  of  burn  patterns  on  exposed  Polaroid 
film  placed  In  the  plane  of  the  sample , because  the 
red  beams  used  for  alignment  of  the  laser  are  not 
coincident  with  the  Infrared  laser  output  after  pas- 
sing through  all  the  switch- out  and  pulse  splitter 
and  delay  optics. 

(c)  Crystal  sample.  The  semiconductor  material 
used  In  these 'experiments  is  intrinsic  germanium. 
The  minimum  resistivity  of  the  high-purity  single 
crystal  sample  is  40.0  fl  cm  at  room  temperature. 
The  sample  was  prepared  from  a crystal  wafer  1 
In.  In  diameter.  First  one  surface  was  polished 
flat  and  smooth  and  then  etched  with  Syton.  The 
purpose  of  etching  la  to  remove  the  damaged  (from 


FIO.  32.  Electro- optical  shutter. 
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mechanical  polishing)  surface  layer  and  clean  the 
surface  of  impurities  in  order  that  the  optical  mea- 
surements reflect  the  bulk  properties  of  the  solid. 
As  is  well  known,  etching  the  sample  surface  has 
an  important  effect  on  the  observation  of  recombin- 
ation radiation  of  germanium  excited  by  a thermal 
f rce.  After  the  first  surface  was  prepared,  the 
sample  was  bonded  to  a KZF-2  glass  substrate, 
which  has  approximately  the  same  coefficient  of 
thermal  expansion  as  germanium  and  is  transpar- 
ent at  the  laser  wavelength.  Then  the  sample  was 
ground  down  to  the  desired  thickness  and  the  sec- 
ond surface  was  also  polished  and  etched.  Sample 
thickness  is  determined  to  be  5.2  pm  by  the  inter- 
ference fringe  spacing  of  the  reflection  spectra 
measured  in  the  infrared  region  of  3-15  iim. 
Crystallinity  of  the  sample  is  verified  by  observing 
the  x-ray  diffraction  pattern  of  the  sample.  During 
the  measurements  the  sample  on  the  substrate  is 
mounted  in  an  evacuated  Dewar  with  quartz  win- 
dows in  order  that  the  sample  can  be  cooled.  The 
sample  temperature  is  measured  by  a thermo- 
couple in  contact  with  the  crystal  surface. 

(d)  Detection  system.  Two  different  types  of 
detectors  are  used  in  these  experiments.  The  en- 
tire pulse  train  (minus  the  switched  out  pulse)  is 
monitored  in  order  to  verify  that  the  pulse  train 
is  well  mode-locked.  Todothls,  the  laser  output 
is  displayed  on  an  ITT  model  F4000  biplanar  vac- 
uum photodiode  and  a Tektronix  model  519  oscillo- 
scope; the  combined  response  time  is  approxi- 
mately 1 nsec. 

The  primary  measurement  in  these  experiments 
is  the  measurement  of  energy  in  a single  (incident, 
reflected,  or  transmitted)  picosecond  pulse.  This 
is  done  by  means  of  United  Detector  Technology 
model  6D  silicon  p-i-n  photodiode  operated  in  the 
pholoconductive  mode.  When  the  laser  pulse  is 
incident  on  the  large  area  photodiode,  electron- 
hole  pairs  created  in  the  intrinsic  region  junction 
arc  swept  out  by  the  reverse  bias  and  charge  an 
integrating  capacitor.  Thevoltage  on  this  capacitor 
is  proportional  to  the  energy  in  the  pulse  and  is 
read  by  means  of  a peak- detector- and-hold  cir- 
cuit and  displayed  on  a digital  voltmeter.  The  op- 


tical collection  system  is  designed  such  that  the 
detectors  intercept  radiation  only  from  a small  sol- 
id angle  completely  filled  by  the  sample.  Electric- 
al shielding  reduces  pickup  of  electrical  noise 
from  the  flash  discharge.  Detection  of  recombina- 
tion radiation  from  germanium  is  eliminated  by 
the  sharp  cutoff  in  spectral  sensitivity  of  the  sili- 
con photodiodes  at  1.1  pm. 

The  pulse  energy  detectors  were  calibrated  both 
absolutely  and  relatively.  Absolute  calibration  was 
performed  against  a calibrated  Quantronix  model 
500  energy  receiver  and  model  504  energy/power 
meter.  Relative  sensitivities  of  the  detector  units 
were  measured  with  the  units  placed  in  their  re- 
spective experimental  configurations. 

(e)  Experimental  error.  We  believe  the  main 
sources  of  error  for  both  excitation  and  excitation- 
probe  measurements  to  be  the  variation  in  pulse 
width  and  in  energy  from  one  laser  firing  to  the 
next.  The  laser  pulse  train  envelope  is  somewhat 
irreproducible,  and  the  pulse  width  and  shape  vary 
with  the  location  of  the  pulse  in  the  tram  of  pulses, 
as  described  by  von  der  Linde.*  This  irreproduci- 
ble and  random  nature  of  the  pulse  evolution  within 
the  laser  cavity  precludes  the  continued  selection 
of  identical  pulses  by  the  laser- triggered  spark 
gap.  The  pulse  width  (unmeasured)  is  believed  to 
range  from  5 to  10  psec  in  duration.  The  depend- 
ence of  the  excitation-probe  data  on  the  level  of  ex- 
citation is  clearly  depicted  in  Fig.  4;  the  theoreti- 
cal dependence  of  the  single  pulse  (excitation  trans- 
mission data  on  pulse  width  are  illustrated  in  Fig. 
25). 

In  addition  to  the  uncertainty  in  pulse  energy  and 
width  from  data  point  to  data  point,  the  transverse 
mode  structure  of  the  laser  is  also  uncertain  (un- 
controlled). Deviations  of  the  transverse  mode 
structure  from  the  TEMW  mode  will  lead  to  "hot" 
spots  on  the  surface  of  the  semiconductor  sample 
when  the  beam  is  focused.  Variations  in  the  posi- 
tions of  these  "hot"  spots  as  the  mode  structure 
changes  from  one  firing  to  the  next  will  result  in 
variations  in  the  degree  of  overlap  between  the  ex- 
citation and  probe  beams  and  consequently  will  re- 
sult in  scatter  in  the  probe  transmission. 
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"The  authors  would  like  to  comment  at  this  point  that 
the  absolute  values  of  the  probe  transmissions  shown 
In  these  figures  are  higher  than  those  experimentally 
observed.  This  Is  a result  of  the  approximations  made 
In  describing  the  Ge  band  structure.  Further  study 
ha*  shown  that  the  absolute  transmission  may  be  ad- 
justed by  romovlng  some  of  these  approximations. 
These  points  are  discussed  In  a following  paper. 

nThe  refractive  Index  for  the  optical  frequency  u0  Is 
approximately  given  by  (t.(l  - coj/u>2)] */ *- 
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ABSTRACT 

Over  the  past  half-decade  measurements  of  the  nonlinear, 
nonequilibrium  optical  properties  of  the  germanium  solid-state 
plasma  have  been  carried  out  in  several  laboratories.  These 
measurements  involve,  for  example,  the  ultrafast  relaxation 
of  optically-excited  electron-hole  distributions 

in  semiconductors  and  the  photoluminescence  spectrum  of  Ge  at 
high  excitation  intensities.  Here,  we  review  our  present 
theoretical  understanding  of  these  experiments  and  we  discuss 
the  theoretical  limitations.  We  also  extend  our  previous  model, 
through  a simple  calculation,  to  include  carrier  diffusion. 
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I.  INTRODUCTION 


I 

I 

The  development  of  lasers  that  can  deliver  high  intensity 
optical  pulses  of  extremely  short  duration  has  renewed  and 
stimulated  the  interest  in  the  measurement  of  the  optical  prop- 
erties of  semiconductors.  Such  measurements  give  information 
about  the  ultrafast  dynamics  of  photogenerated  electron-hole 
plasmas,  as  well  as  the  band  structure  of  the  semiconductors. 

In  this  paper,  we  review  some  of  our  experiments  in  which 

picosecond  optical  pulses  were  used,  their  results,  and  our 
theoretical  work  concerning  them.  We  also  present  the 
results  of  a simple  calculation  that  illustrates  the  effect 
of  including  diffusion  in  our  previous  model.  Other  major 
work  in  this  area,  performed  at  Bell  Laboratories,  will  be 
reviewed  by  D.  H.  Auston  in  an  invited  paper  elsewhere  in 
this  issue. 

In  this  paper,  we  will  separate  our  discussion 
into  two  parts.  The  first  part  (Section  II)  describes  the 
photogeneration  of  the  carrier  plasma  by  an  intense  ultrashort 
excitation  pulse  and  the  evolution  of  that  plasma  during  the 
time  this  excitation  pulse  is  incident  on  the  sample.  The 
second  part  (Section  III)  describes  the  temporal  evolution  of 
that  photoqenerated  plasma  after  the  excitation  pulse  has 
passed,  as  monitored  by  a weak  probe  pulse.  In  addition,  we 
also  present,  in  Section  III,  a simple  analysis  of  the  possible 

j effects  of  carrier  diffusion  on  these  optical  measurements. 

I 
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II.  OPTICAL  GENERATION  OF  HIGH  DENSITY,  HOT  ELECTRON  HOLE 
PLASMAS 

In  this  section  we  review  experiments  in  which  the  electron- 
hole  plasma  is  generated  by  the  absorption  of  a single,  ultra- 
short  optical  pulse  in  a thin  semiconductor  wafer.  The  trans- 
mission of  this  single  pulse  through  the  sample  provides  infor- 
mation concerning  the  generation  and  evolution  of  the  plasma 
during  the  period  the  pulse  is  present  in  the  semiconductor. 

The  nonlinear  transmissions  of  thin  germanium  wafers  have 

1 2 

been  measured  by  Kennedy  et  al.  and  Smirl  et  al.  using  ultra- 

short  (5-10  psec)  pulses  at  1.06pm  from  a mode-locked  Nd:glass 

laser.  Figure  1 shows  the  transmission  of  a 5.2pm-thick  Ge 

sample  as  a function  of  incident  pulse  energy  for  two  sample 

2 

temperatures,  77K  and  300K.  These  curves  display  two  regions. 
There  is  a region  in  which  the  transmission  is  nearly  constant 
and  there  is  a nonlinear  absorption  region  in  which  the  trans- 
mission rises  slowly  and  appears  to  peak  where  the  crystal  is 
damaged  by  the  optical  pulse.  Similar  measurements  have  been 
performed  on  thin  samples  of  Hg^_xCdxTe  by  Matter  et  al.3,  as 
shown  in  Fig.  2.  We  shall  concentrate  on  the  Ge  results,  since 
the  theoretical  model  for  Ge  is  more  nearly  complete. 

Recently,  Elci  et  al.*  have  presented  a theoretical  model 
that  accounts  for  the  transmission  of  single  optical  pulses 
through  thin  Ge  samples  in  terms  of  the  energy  band  structure 
and  other  well-known  parameters.  The  band  structure  of  Ge  is 
shown  in  Fig.  3.  The  data  presented  in  Fig.  1 can  be  accounted 
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for  in  terms  of  this  model  in  the  following  way.  When  the 

optical  pulse  enters  the  Ge  crystal,  it  is  absorbed  by  direct 

transitions.  Thus,  it  creates  a large  number  of  electrons 

(holes)  in  the  central  valley  of  the  conduction  (valence) 

band.  The  electrons  located  in  the  central  valley  of  the  con- 

-14 

duction  band  are  rapidly  (^lO  sec)  scattered  to  the  conduc- 
tion band  side  valleys  by  long  wave  vector  phonons.  Since 
these  electrons  are  emptied  from  the  central  valley  to  side 
valleys  at  a rate  that  is  comparable  to  the  direct  absorption 
rate,  any  decrease  in  the  number  of  states  available  in  the 
central  valley  for  direct  absorption  is  ultimately  determined 
by  a buildup  of  the  populations  in  the  side  valleys.  Carrier- 
carrier  scattering  events,  which  also  occur  at  a rate  comparable 
to  the  direct  absorption  rate,  ensure  that  the  electron  and 
hole  distributions  will  be  Fermi-like.  They  also  ensure  that 
the  Fermi  distribution  for  holes  and  the  Fermi  distribution 
for  electrons  will  reach  a common  temperature  that  is  in  general 
greater  than  the  lattice  temperature.  Since  the  photon  energy 
ftu>o  is  greater  than  either  the  direct  band  gap  Eq  or  the  indirect 
gap  Eg,  a direct  absorption  event  followed  by  phonon-assisted 
scattering  of  the  electron  from  the  central  to  side  conduction 
band  valleys  results  in  an  excess  energy  fioo  -E_  being  given 
to  thermal  agitation.  This  excess  energy  results  in  a carrier 
temperature  due  to  direct  absorption  that  is  much  higher  than 
the  lattice  temperature.  Thus,  direct  absorption  results  in 
the  creation  of  a large  number  of  electrons  (holes)  in  the  con- 
duction (valence)  band  with  a high  distribution  temperature. 
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Tf  the  above  mentioned  processes  were  the  only  ones  included 
in  the  model,  the  transmission  of  the  thin  Ge  wafer  would 
begin  at  its  Beer's  law  value  and  rise  as  a function  of  in- 
cident pulse  energy  because  of  the  saturation  of  the  optically 
coupled  states  as  a result  of  direct  absorption.  The  observed 
rise  in  transmission  as  a function  of  incident  pulse  energy 
is,  however,  much  slower  than  the  band  filling  due  to  direct 
absorption  would  predict.  The  observed  slower  rise  is  caused 
by  many  body  effects  in  the  electron-hole  plasma,  as  will  be 
discussed  below. 

As  a result  of  direct  absorption,  the  conduction  (valence) 
band  of  the  sample  contains  a large  number  of  electrons  (holes) 
with  a high  distribution  temperature.  Electrons  (holes)  loca- 
ted hig^in  the  conduction  (valence)  band  can  relax  within  that 
band  by  phonon  emission.  The  effect  of  this  relaxation  is  to 
reduce  the  carrier  temperature  and  increase  the  lattice  tem- 
perature. As  the  distribution  cools  by  means  of  these  intra- 
band phonon-assisted  transitions,  carriers  fill  the  states 
needed  for  absorption.  Thus,  phonon-assisted  relaxation  tends 
to  increase,  or  further  enhance,  the  sample  transmission. 
Phonon-assisted  relaxation  occurs  on  a time  scale  comparable 
to  the  optical  pulse  width. 

Many  body  effects  in  the  electron-hole  plasma,  which  be- 
come significant  as  the  density  of  the  carriers  created  by 
the  optical  pulse  increases,  cause  the  experimentally-measured 
transmission  as  a function  of  incident  optical  pulse  energy 
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to  rise  more  slowly  than  predicted  by  band  filling  caused  by 
direct  absorption  alone.  In  particular,  an  electron  or  a 
hole  moves  in  the  screened  coulomb  field  of  other  electrons 
and  holes.  The  quantized  long  range  part 

of  the  screened  coulomb  field  is  a plasmon.  As  the  carrier 
density  increases,  because  of  direct  absorption,  the  plasma 
frequency  of  the  carriers  increases.  If  the  carrier  density 
is  sufficiently  high,  an  electron  in  the  central  conduction 
band  valley  can  recombine  with  a hole  near  the  top  of  the 
valence  bands  while  emitting  a plasmon.  Normally,  plasmon- 
assisted  recombination  (PAR)  can  occur  only  if  the  plasma 
frequency  u>  is  larger  than  the  direct  band  gap  frequency 
E0/ft.  However,  in  our  case,  owing  to  the  strong  perturbation 
of  the  solid  by  the  optical  pulse,  the  plasmon  resonance  is 
broadened  considerably.  Therefore,  PAR  can  occur  at  plasma 
frequencies  below  E^/fi.  These  collective  oscillations  have 
extremely  short  lifetimes  and  rapidly  decay,  transferring 
their  energy  to  electrons  and  holes,  thus  increasing  the  carrier 
temperature.  Therefore,  as  the  carrier  number  increases,  the 
PAR  rate  increases , which  retards  the  growth  of  the  carrier 
number  and  raises  the  carrier  temperature.  Thus,  further  in- 
crease in  the  sample  transmission  is  slowed  as  the  electrons 
(holes)  are  heated  and  removed  from  the  optically  coupled 
states.  Free  carrier  absorption,  although  less  effective, 
also  serves  to  increase  the  distribution  temperature  while 
the  pulse  is  incident  on  the  sample. 
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Thus,  we  see  that,  physically,  the  single  pulse  trans- 
mission of  Ge  as  a function  of  incident  pulse  energy  (see  Fig. 

1)  can  be  accounted  for  by  direct  interband  transitions  fol- 
lowed by  heating  of  the  carrier  distribution.  The  solid  lines 
in  Fig.  1 represent  the  calculated  transmissions  using  the  model 

4 

of  Elci  et  al.  The  agreement  between  theory  and  experiment 
is  quite  good. 

Plasmons  play  an  important  role  in  determining  the  trans- 
mission properties  of  thin  Ge  wafers  on  picosecond  time  scales. 
Luckily,  there  is  a property  of  PAR  that  makes  it  possible  to 
test  experimentally.  PAR  is  sensitive  to  tiu)p/Eo.  If  the  direct 
band  gap  Eq  could  be  increased  slightly  while  keeping  the  car- 
rier density  relatively  constant,  the  PAR  rate  should  decrease 
appreciably.  This  would  result  in  a decreased  plasma  temper- 
ature for  a given  pulse  energy,  and  one  would  expect  the  sat- 
urated transmission  to  increase  owing  to  the  carriers  being 
located  lower  in  the  band.  It  is  possible  to  increase  Eq  by 
subjecting  Ge  to  pressure.  Van  Driel  et  al.^  have  measured 
the  nonlinear  transmission  of  picosecond  optical  pulses  (1.06un\) 
through  a thin  Ge  sample  as  a function  of  hydrostatic  pressure 
up  to  24  kbars,  as  shown  in  Fig.  4.  At  high  intensities,  the 
degree  of  enhanced  transmission  is  observed  to  increase  then 
decrease  as  the  band  gap  is  tuned  with  pressure.  In  terms  of 
our  model,  the  increase  and  then  decrease  in  enhanced  transmission 
with  increasing  pressure  is  caused  by  a decrease  in  the  PAR  rate, 
plasma  temperature,  and  carrier  density.  The  initial  increase 
in  the  degree  of  enhanced  transmission  as  the  pressure  is  in- 
creased from  1 bar  is  primarily  a result  of  a decrease  in  the 
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PAR  rate  caused  by  an  increase  in  the  direct  gap  Eq,  which 

allows  for  a slightly  more  dense  but  cooler  plasma.  This 

results  in  the  optically  coupled  electron-hole  states  being 
completely 

more^ occupied  than  at  the  lower  pressure  and,  therefore, 

causes  a higher  transmission.  This  trend  continues  up  to 

a pressure  of  approximately  10  kbar.  Above  this  pressure, 

other  effects  begin  to  dominate.  For  example,  the  direct 

absorption  coefficient  in  the  Beer's  law  region  will  decrease 

with  increasing  pressure,  reflecting  a decrease  in  the  joint 

density  of  states,  for  flcu  = 1.17  eV,  as  the  valence  and  con- 

o 

duction  bands  move  apart.  This  decrease  in  absorption  coef- 
ficient results  in  a decrease  in  the  carrier  density  and, 
thus,  the  plasma  frequency.  A more  subtle  effect  in  the  re- 
duction of  the  efficiency  of  the  PAR  process  occurs  because 
the  indirect  energy  gaps  do  not  increase  as  fast  as  the  direct 
gap  with  increasing  pressure.  Therefore,  the  ratio  of  the 
central  to  side  valley  electron  populations  decreases  with 
increasing  pressure.  Since  the  efficiency  of  the  PAR  process 
depends  on  the  number  of  carriers  in  the  central  valley  of  the 
conduction  band,  the  number  of  PAR  events  decreases  as  well. 

Because  of  the  role  of  the  conduction  band  side  valleys 
in  the  present  interpretation  of  the  single  pulse  transmission 

data,  van  Drielu  have  measured  the  photoluminesence  spectra  of 
et  al . 

r,e  as  a function  of  excitation  level  at  1.06pm.  Figure  5 shows 
the  results  such  experiments  performed  at  sample  temperatures 
of  295  and  174K.  Four  different  excitation  irradiances  (50, 
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500 , 10  and  10  W/cm  ) were  uBed  at  each  temperature. 

The  data  at  both  temperatures  show  a peak  at  1.76um  cor- 
responding to  the  indirect  gap  recombination.  In  addition, 
there  is  a second  peak  at  1.5um,  corresponding  to  direct  gap 
transitions.  The  direct  gap  photoluminescence  intensity  in- 
creases with  increasing  excitation  level.  The  main  features 
of  the  data  can  be  understood  by  a simple  argument.  Recall 
that  electron  scattering  from  the  central  to  side  conduction 
band  valleys  by  long  wave  vector  phonons  is  extremely  rapid. 
Thus,  any  buildup  of  the  central  valley  electron  population 
must  be  accompanied  by  a buildup  of  the  populations  in  the 
ten  side  valleys.  At  low  excitation  levels,  the  number  of 
electrons  created  is  so  small  that  an  insignificant  number 
exists  in  the  central  T valley;  nearly  all  recombination  is 
indirect.  At  higher  excitation  levels,  the  number  of  electrons 
in  the  central  valley  increases.  Although  this  number  is  al- 
ways much  less  than  the  number  in  the  side  valleys,  the  re- 
combination rate  across  the  direct  gap  far  exceeds  that  across 
the  indirect  gap  because  the  latter  requires  phonon  partici- 
pation. It  is  therefore  possible  with  increasing  carrier 
density  to  observe  more  direct  than  indirect  gap  luminescence. 
The  theory,  shown  by  solid  lines  in  Fig.  5 , is  seen  to  be  in 
good  agreement  with  experiment. 

We  now  make  a few  additional  comments  concerning  certain 
assumptions  made  in  this  model.  An  important  assumption  made 

4 

in  our  calculations  is  that  of  local  thermodynamic  equilibrium. 


Although  we  have  ured  the  rapid  rate  for  carrier-carrier  col- 
lisions as  a justification  for  taking  the  carrier  distribution 
functions  to  be  Fermi-Dirac , it  is  none-the-less  an  assumption. 

The  implications  of  this  assumption  are  currently  being  in- 
vestigated. If  the  distributions  are  assumed  to  be  Fermi- 
Dirac,  the  calculations  are  greatly  simplified.  The  plasma 
can  then  be  described  by  a local  temperature  and  Fermi  energies. 
To  this  date,  good  results  have  been  obtained  with  this  assump- 
tion, as  the  experiments  described  above  illustrate. 

A more  drastic  assumption  is  to  neglect  the  spatial  vari- 
ation of  the  parameters  (temperature  and  Fermi  energies)  that 
characterize  the  electron-hole  plasma  throughout  the  inter- 
action volume  of  the  thin  (approximately  5ym-thick)  Ge  wafer. 
Thus,  in  the  interaction  volume,  electrons  and  holes  are  des- 
cribed by  Fermi-Dirac  distributions  with  appropriate  Fermi 
energies  for  electrons  and  holes  and  a single  temperature  for 

both,  all  of  which  are  independent  of  position.  Therefore,  the 
the  the 

effects  of^longitudinal  diffusion  of^carriers  and  of  the  optical 
pulse  propagation  are  neglected.  This  is  a severe  limitation 
of  the  model  and  can  be  relaxed  by  allowing  the  plasma  parameters 
to  be  dependent  on  both  space  and  time.  We  include,  in  a sim- 
plistic manner,  these  spatial  effects  in  the  model  for  the  first 
time  in  the  next  section  of  this  paper.  Other  approximations 
and  simplifications  to  the  model  are  discussed  by  Latham  et  al.7 
in  another  paper  in  this  issue. 


III.  EVOLUTION  OF  THE  PHOTOGENERATED  PLASMA 
The  electron-hole  plasma  created  by  an  intense  picosecond 
optical  pulse  evolves  rapidly  after  the  passage  of  that  pulse, 
owing  to  the  rapid  relaxation  processes  that  take  place  in 
semiconductors.  Substantial  information  concerning  this  tempo- 
ral evolution  can  be  obtained  by  measuring  the  transmission 
of  a second,  weak  picosecond  pulse  that  is  delayed  in  time 
with  respect  to  the  excitation  pulse.  In  this  section,  we 
briefly  review  the  experiments  that  monitor  the  temporal 
evolution  of  the  hot  plasma  after  the  passage  of  the  excitation 
pulse  and  our  present  understanding  of  these  experiments.  In 
addition,  through  a simple  calculation,  we  consider  the  effect 
of  diffusion  on  the  present  model. 

The  relative  transmission  of  a probe  pulse 

as  a function  of  time  delay  after  an  excitation  pulse  for  Ge 
sample  temperatures  of  77K  and  297K  is  shown  in  Fig.  6 (from 
Ref.  2).  In  these  experiment^  the  probe  pulse  had  5%  of  the 
energy  of  the  excitation  pulse.  The  excitation  pulse  was  in 
the  region  of  nonlinear  transmission  and  the  probe  was  in  the 
region  of  linear  transmission  (see  Fig.  1 ).  Note  that  the 
probe  transmission  rises  rapidly,  peaks,  and  then  turns  over 
for  delays  of  the  order  of  100  psec.  The  solid  lines  in  Fig.  6 

are  the  probe  transmissions  calculated  from  the  model  of  Elci 

4 -assisted  carrier 

et  al.  , in  which  only  the  intraband  phonon. relaxation  is  taken 

A 

into  account;  diffusion  and  other  possible  interactions  are 
neglected.  In  view  of  these  approximations,  the  agreement 
achieved  is  good. 
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The  probe  transmission  can  be  understood  in  terms  of  the 

4 

model  of  Elci  et  al.  in  the  following  way.  After  the  passage 

of  the  excitation  pulse,  the  interaction  region  of  the  sample 

contains  a large  number  of  carriers  with  a high  distribution 

temperature.  The  electrons  (holes)  are  located  high  (low) 

in  the  bands  because  of  the  high  distribution  temperature, 

leaving  the  states  that  are  optically  coupled  available  for 

the 

absorption.  Thus,. probe  transmission  is  small.  Later,  as 
the  distribution  cools  by  means  of  intraband  phonon-assisted 
transitions  and  carriers  fill  the  states  needed  for  absorption, 
the  probe  transmission  increases.  Note  that  the  PAR  process 
is  essentially  turned  off  when  the  excitation  pulse  has  passed 
through  the  sample. 

One  might  plausibly  attribute  the  decay  of  the  probe 

8 9 

transmission  for  long  delays  (Fig.  6 ) to  carrier  diffusion.  ' 
Because  of  the  large  focused  spot  size  for  the  optical  pulse 
(250pm  in  diameter) , diffusion  transverse  to  the  direction 
of  light  propagation  is  entirely  neglible.  However,  diffusion 
of  carriers  from  a region  near  the  sample  surface  into  the  bulk 
crystal  in  the  direction  of  light  propagation  can  be  signifi- 
cant. We  term  the  diffusion  in  the  direction  of  light  prop- 
agation longitudinal.  Reflectivity  experiments  that  measure 
the  longitudinal  diffusion  of  an  optically-created  surface 

carrier  density  into  the  crystal  bulk  have  been  carried  out 

0 

by  Auston  and  Shank.  However,  note  that  longitudinal  carrier 
diffusion  will  not  necessarily  affect  surface  and  bulk  measure- 
ments in  the  same  manner.  When  a 


probe  pulse 
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transmitted 

is  through  a thin  sample,  it  sees  the  same  total  number 

A 

of  carriers  (neglecting  Auger  and  other  recombination)  for 
all  relevant  delays;  however,  the  carriers  will  be  spatially 
rearranged  because  of  diffusion.  To  investigate  how  the 
carrier  diffusion  affects  the  bulk  transmission  measurements, 
we  perform  a simple  calculation  based  on  the  calculations  of 
Elci  et  al.* 

To  study  the  effect  of  carrier  diffusion  alone , we  assume 
that  the  plasma  is  in  thermal  equilibrium  with  the  lattice. 

This  is  the  case  for  delay  times  greater  than  100  psec,  since 
on  these  time  scales  phonon-assisted  relaxation  is  essentially 
complete.  This  plasma  , although  inhomogeneous  in  spatial 

distribution,  has  a uniform  temperature  T = T^  (lattice  tem- 
perature) . We  will  assume  local  charge  neutrality  and  ambi- 
polar  diffusion.  Thus,  electron  and  hole  densities  are  equal 
everywhere.  Electrons  and  holes  have  different  Fermi  energies 
EF  and  E^ , respectively,  which  are  space  and  time  dependent, 
although  Ep  and  are  related  since  electron  and  hole  densities 
are  equal.  Finally,  we  will  assume  that  the  electron  (or  hole) 
density  is  space  and  time  dependent  only  via  E,,  (or  E„)  . Since 
diffusion  is  longitudinal,  we  have  a one-dimensional  diffusion 
problem,  where  Ep  = Ep(x,t)  and  EH  = F,H(x,t).  The  electron 
density  and  hole  density  are  given  by^ 
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Definitions  of  the  various  symbols  and  functions  used  in  this 
and  the  subsequent  formulae  are  given  in  Tables  I and  II.  Also 
the  carrier  density,  n(x,t),  must  satisfy  the  diffusion  equa- 
tion 


9n 

St 


(2) 


We  will  not  attempt  to  calculate  the  ambipolar  diffusion  con- 
stant D.  Instead,  we  will  use  the  value  measured  by  Shank  and 

g 

Auston.  If  an  initial  density  n(x,o)  is  given,  (1)  and  (2) 
determine  the  subsequent  E_(x,t)  (or  E„(x,t) ) . 

Since  the  optical  pulses  used  for  probing  are  weak,  we  can 
neglect  the  perturbation  of  the  distributions  induced  by  the 
probe  in  the  calculation  of  its  transmission.  Also,  since  the 
probe  pulse  duration  is  short  compared  to  the  diffusion  time 
scale,  we  can  neglect  the  effect  of  the  finite  probe  duration 
on  the  probe  transmission  as  well.  Within  these  approximations, 
the  probe  transmission  is  given  by 

2 -/J*  dxa  (x,t) 

Tpr  (t)  - (1  - R)  e ° , (3) 

where  n(x,t)  is  the  absorption  coefficient,  and  a(x,t)  depends 
on  x and  t via  Ep  and  EH«  The  absorption  coefficient  at  as  a 
function  of  Ep  and  E^  is  calculated  in  Ref.  4.  It  is  composed 
of  two  parts: 


a(x,t) 


aD(x,t)  + apCA(x,t) 


9 


(4) 
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where  r»D  is  due  to  the  direct  electronic  transitions  from  the 
valence  bands  to  the  central  conduction  band  valley  of  Ge 
and  a is  due  to  the  free  carrier  absorption.  Expressions 

r V./» 

4 

for  aD  and  «FCA  are  given  by  Elci  et  al.  as 


aD(x,t)  = aQ[  1 - - F (y  2 ) ] 


aFCA(x't)  = al[Co(lJ3)  + 2 


% 


C0(^4)] 


+ a2f  ( y 3 ) + 2 


“h 


5/2 


(U4 ) ] . 


(5) 


(6) 


In  order  to  calculate  the  carrier  distribution  n(x,t)  at 
some  later  time,  we  need  the  initial  distribution  n (x , t(jelay=0^  ' 
the  spatial  distribution  of  carriers  as  they  are  generated  by 
the  excitation  pulse.  We  are  currently  working  on  a model 
that  will  give  n ^elay”0*  * However*  for  the  present,  to 
avoid  considering  the  temperature  dependence  of  the  electron- 
hole  plasma,  we  take  the  time  origin  to  be  such  that  temperature 
relaxation  is  complete.  Thus,  n(x,0)  = n (x , tde^ay“100psec)  is 
used  as  the  initial  condition  for  solving  the  diffusion  equation. 
Of  course,  n (x , tde^ay»100psec)  is  itself  not  known;  we  must 
assume  an  initial  distribution  and  investigate  the  effect  of 
its  temporal  behavior  on  the  sample  transmission.  Again,  we 
emphasize  that  we  do  not  intend  this  to  be  an  exact  calculation, 
but  rather  a simple  investigation  of  the  possible  effects  of 
diffusion  for  long  time  delays. 
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To  understand  the  effects  of  diffusion  on  the  probe  trans- 
mission, we  perform  the  following  calculation.  Given  the  aver- 
age density  of  carriers  created  by  the  excitation  pulse  in  the 

4 

interaction  region,  n,  (see  Elci  et  al.  ) and  assuming  a spat- 
ially dependent  initial  carrier  density  n(x,0),  the  diffusion 
equation,  Eq.  (2) , can  be  solved  for  the  subsequent  carrier 
distribution  n(x,t).  Results  of  a typical  calculation  are 
shown  in  Fig.  7.,  where  an  initial  distribution  of  the  form 


n (x,  0) 


no  e 


X//W  is  assumed.  Since  the  electron  temperature 


is  equal  to  the  lattice  temperature,  the  Fermi  energies  Ep(x,t) 


and  E„(x,t)  can  be  calculated  from  Eq.  (1) 

H 


Fermi  energies,  (*D  and 


Next , knowing  the 
can  be  calculated  from  Eq.  (5)  and 


Eq.  (6) , and  the  transmission  can  be  calculated  from  Eq.  (3) . 

absorption 

The  spatially  dependent  direct  and  free  carrier^coef f icients , 


calculated  from  the  carrier  densities  of  Fig.  7 , are  shown 

in  Fig.  8 and  Fig.  9 , respectively.  Note  that  the  direct 

absorption  coefficient  is  always  much  larger  than  the  free 

carrier  absorption  coefficient.  Finally,  Fig.  10  shows  the 

resulting  probe  transmission  for  times  greater  than  100  psec 

and  for  various  initial  density  distributions.  The  average 

19  3 

electron  density  for  all  curves  in  Fig.  7-10  is  2.5x10  /cm  . 

Figure  10  illustrates  that  for  low  average  carrier  den- 

the 

sities,  diffusion  will  tend  to  decrease  probe  transmission,  as 
might  be  expected.  However,  when  n is  large  enough,  diffusion 
can  actually  cause  an  increase  in  probe  transmission,  as  shown 

in  Fig.  11.  While  this  may  seem  surprising  at  first,  it  be- 
comes more  reasonable  if  the  following  physical  argument  is 
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considered.  Consider  a simple  schematic  of  the  diffusion 
process  such  as  the  one  shown  in  Fig.  12.  Since  n stays 
constant,  the  probe  pulse  sees  the  same  number  of  carriers 
independent  of  time  delay.  However,  the  states  that  are 
optically  coupled  are  restricted  to  narrow  regions  in  energy 
in  the  band  structure.  So,  if  we  consider  only  direct  ab- 
sorption of  the  probe,  not  all  carriers  are  effective  in 
filling  the  optically  coupled  states  and  preventing  absorp- 
tion. Thus,  it  is  not  right  to  say  that  aD*  for  the  probe 

L 

is  proportional  to  f n(x,t)dx,  the  total  number  of  carriers, 
rather. 


neff(x,t)  dx  = Neff 


(7) 


where  N ,,  denotes  the  total  number  of  carriers  effective  in 
err 

clogging  the  optically  coupled  states  needed  for  absorption. 

Neff'  tota*  num*>er  of  carriers  effective  in  preventing 

light  absorption,  is  altered  by  diffusion  as  shown  schematically 
in  Fig.  12.  If  the  number  of  carriers  is  large,  when  they  mi- 
grate from  the  front  surface,  they  can  fill  the  states  needed 
for  absorption  away  from  the  front  surface  without  depleting 

the  optically  coupled  states  near  the  surface;  N increases 

ef  r 

and  absorption  decreases.  If  n is  small,  diffusion  will  de- 
crease neff  in  the  front  region  without  increasing  n f sig- 
nificantly in  the  back,  thus  the  transmission  decreases.  Thus, 
depending  on  the  total  carrier  number,  longitudinal  diffusion 


can  cauae  the  probe  transmission  to  rise  or  fall.  Experi- 
mentally, we  have  not  observed  a rise  for  long  probe  delays 
at  present  densities. 

Diffusion  should  also  affect  the  probe  transmission  for 
shorter  delay  times  as  well.  That  is,  carriers  diffuse  while 
they  lose  energy  to  the  lattice.  However,  when  the  carriers 
are  hot,  D is  smaller  than  its  thermal  equilibrium  value  for 
the  same  density.  In  a precise  calculation,  one  would  have 
to  allow  the  diffusion  coefficient  D to  depend  on  time  through 
the  carrier  temperature. 

Excitation-probe  measurements  have  been  performed  on 
Hg1_j(CdxTe3  as  well  as  Ge.  The  results  of  these  measurements 
are  presented  in  Fig.  13  for  comparison  with  the  Ge  results 
and  for  completeness. 


IV.  SUMMARY 


In  this  paper,  we  have  reviewed  recent  experiments  that 

measure  the  nonlinear,  nonequilibrium  properties  of  germanium 

and  mercury  cadmium  tellur ide.  In  addition,  we  have  described 
germanium 

how  the  experiments  can  be  accounted  for  by  a recently  pro- 

A 

posed  model  that  describes  the  generation  and  temporal  evolu- 
tion of  hot  electron-hole  plasmas  produced  in  semiconductors 
by  the  absorption  of  intense,  ultrashort  optical  pulses.  In 
these  calculations,  however,  several  assumptions,  omissions, 
and  simplifications  are  made.  The  major  assumptions  are  the 
following.  The  carrier  distributions  are  assumed  to  be  Fermi- 
Dirac  for  all  time  scales  of  interest.  All  carrier  densities 
and  Fermi  energies  are  allowed  to  depend  on  time  but  not  posi- 
tion. The  overall  good  agreement  between  the  theory  and  pre- 
sent experiments  indicate  that  these  assumptions  work  quite 
well  despite  the  intensities  and  short  durations  of  the  optical 
pulses.  However,  more  exact  calculations  and  definitive  experi- 
ments are  in  progress  in  an  attempt  to  investigate  the  impli- 
cations of  these  approximations.  These  calculations  will  also 
include  processes  omitted  from  the  original  model,  such  as 
Auger  recombination  and  diffusion.  We  believe  it  is  clear  that 
future  investigations  with  picosecond  optical  pulses  will  yield 
information  concerning  the  individual  characteristics  of  semi- 
conductors and  the  dynamics  of  the  fundamental  processes 
that  take  place  in  semiconductors. 
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II.  Definition  of  Integrals  and  Constants 
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TABU;  1 


w light  frequency 

c electronic  charge 

ca  high  frequency  dielectric  constant 

c speed  of  light 

c.  sound  velocity 

A 

m electronic  mass 

m electronic  effective  mass  in  the  central 

valley  of  the  conduction  band 

m^  density  of  states  effective  mass  in  L-valleys 

m^  effective  hole  mass 

T temperature 

E direct  gap 

o 

Eg  indirect  gap  / 

\i’  ^d  electron-acoustic  phonon  coupling  coefficients 
Qo  electron-optical  phonon  coupling  coefficient 

p mass  density 

fiQ  optical  phonon  frequency 

R reflection  coefficient 

l sample  thickness 

D diffusion  constant 


•fiui  = 1.17  cV 

16 

5.4  » 10^  cm-sec  1 

0.1  m 

0.22  m 
0.34  m 
300°K 

0.889  eV  (at  300°K) 
0.664  eV  (at  300°K) 
17  eV,  -3.4  eV 

6 * 10  erg  cm 

c t -3 

b.o  gm-cm 

hfl  = 0.02  eV 
o 

0.34 
5.5  pm 

2 -1 
230  cm  -sec 
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PIGURE  CAPTIONS 

1.  Single  pulse  transmission  of  a 5.2ym-thick  single  crystal 
germanium  sample  versus  incident  energy  in  units  of  quanta 
at  1.06ijm  for  temperatures  of  297K  and  77K. 

2.  Single  pulse  transmission  of  a Sym-thick  Hg.  ,QCdn  ,.Te 
sample  versus  incident  energy  in  units  of  quanta  at  1.06ym 
for  temperatures  of  297K  and  105K. 

3.  Ge  energy-band  (eV)  structure  at  300K. 

4.  Transmission  of  a 5.2ym-thick  germanium  sample  as  a function 
of  incident  energy  and  pressure,  for  a sample  temperature 

of  295K. 

5.  Emission  spectra  of  germanium  for  different  excitation 
irradiances  IQ  at  temperatures  of  (i)  295K  and  (ii)  174K. 

6.  Probe  pulse  transmission  of  a 5.2ym-thick  germanium  sample 
versus  time  delay  for  sample  temperatures  of  297K  and  77K. 
Relative  transmission  units  are  given  by  the  ratio  of  the 
probe  transmission  to  excitation  transmission  for  each 
delay,  normalized  such  that  the  peak  of  the  77K  curve  is 
unity. 

7.  The  variation  due  to  diffusion  of  the  carrier  density,  n(x,t) 
as  a function  of  longitudinal  position  and  time  in  a 5.5ym- 
thick  germanium  sample.  The  average  density  is  determined 
by  the  total  carrier  number/interaction  volume. 

8.  The  variation  due  to  diffusion  of  the  direct  absorption 
coefficient  as  a function  of  longitudinal  position  and  time 
in  a 5. Sum-thick  germanium  sample. 
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9.  The  variation  due  to  diffusion  of  the  free  carrier  absorp- 
tion coefficient  as  a function  of  longitudinal  position 
and  time  in  a 5.  Spin-thick  germanium  sample. 

10.  The  variation  of  the  probe  transmission  with  time  as  a 
result  of  diffusion,  for  an  average  carrier  denisty  of 
2.5xl019/cm3. 

11.  The  variation  of  the  probe  transmission  with  time  as  a 

result  of  diffusion,  for  an  average  carrier  density  of 
19  3 

4.75xlOA  /cm  . 

12.  Schematic  of  the  diffusion  process 

13.  Probe  pulse  transmission  of  a Spin-thick  HgQ  39CdQ  giTe 
sample  versus  time  delay  for  sample  temperatures  of  297K 
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Abstract 


The  nonlinear  transmission  of  ultrashort  optical  pulses 
(1  >1.06um)  through  a thin  germanium  sample  has  been  measured 
as  a function  of  pulsewidth.  These  measurements  provide  in- 
sight into  a recently  proposed  theoretical  model  that  describes 
the  enhanced  transparency  observed  at  high  intensities. 
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Recently  Elci  et  el.1  have  presented  a theoretical  model  that 
accounts  for  tie  ultrafast  transient  response  of  optically-created,  hot 
electron-hole  plasmas  in  germanium.  This  model  is  shown  to  be  in  sub- 
stantial agreement  with  the  observed  nonlinear  transmission  of  single, 
ultrashort  optical  pulses  through  germanium  as  a function  of  incident 
pulse  energy  at  temperatures  of  105  K and  297  K and  with  the  transmission 

of  a weak  probe  pulse  as  a function  of  picosecond  time  delays  after  an 

2 

intense  excitation  pulse.  In  addition,  the  theory  agrees  well  with 
experimental  nonlinear  transmission  data  obtained  as  a function  of  band 
gap  energy  achieved  by  pressure  tuning  a germanium  sample.3  One  feature 
of  this  model  is  the  predicted  strong  pulsewidth  dependence  of  the  non- 
linear transmission  of  these  intense#  ultrashort  optical  pulses  through 
germanium. 

This  pulsewidth  dependence  is  depicted  in  Fig.  1.  This  figure  is  a 
plot  of  the  theoretical  transmission  of  a 5.2pm  - thick  sample  of  ger- 
manium as  a function  of  the  number  of  incident  quanta  for  a sample  tem- 
perature of  77  K and  for  theoretical  optical  pulse  widths  of  2.5,  5, 
and  10  psec.  Elcl  et  al.1  have  suggested  that  experiments  with  variable 
duration  pulses  be  performed  to  provide  a further  test  for  the  theoretical 
model.  In  this  paper,  we  report  the  results  of  such  measurements. 

A schematic  of  the  experimental  apparatus  is  shown  in  Fig.  2.  A 
single  pulse,  having  a wavelength  of  1.06pm  and  a pulse  width  t of  5-10 
psec,  is  switched  from  a train  of  pulses  from  a mode-locked  Nd: glass  laser 
by  an  electro-optic  shutter.  To  produce  variable-width  pulses,  this 
single  optical  pulse  is  split  into  two  separate  parts  that  are  then  re- 
combined collinearly  after  one  part  has  been  delayed  with  respect  to  the 
other  with  a Michelson  interferometer.  The  recombined  pulse  is  then  focused 
on  the  sample  to  a spot  with  diameter  of  approximately  250pm.  The  sample 
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is  an  undoped,  high  purity  (Pm^n  “ 40.0  (lcm) , single  crystal  germanium 
wafer  that  is  cut  with  the  (111)  plane  as  the  face,  mounted  on  a KZF-2 
substrate,  and  polished  to  a thickness  of  7. Sum,  as  determined  from  an 
interferometric  technique. 

Figure  3 is  a plot  of  the  transmission  versus  incident  pulse  energy 
for  two  pulse  widths,  t and  2t , and  for  a sample  temperature  of  105  K. 

Each  point  shown  is  the  average  of  8 to  12  data  points,  and  the  error 
bars  indicate  one  standard  deviation.  Although  the  exact  pulse  width  x 
is  uncertain  from  one  optical  pulse  to  the  next,  it  can  be  readily  seen 
that  the  difference  between  the  two  curves  is  less  than  the  uncertainties 
in  the  data  and  substantially  less  than  that  predicted  by  theory. 

The  pulsewidth  dependence  of  the  transmission  can  be  understood  in 
terms  of  this  model  in  the  following  way.  Initially,  when  the  optical 
pulse  enters  the  Ge  sample  most  of  it  is  absorbed  creating  a large  number 
of  electrons  (holes)  in  the  central  valley  of  the  conduction  (valence) 
band.  These  electrons  are  rapidly  (compared  to  a psec)  scattered  to  the 
conduction  band  side  valleys  by  phonons.  Carrier-carrier  scattering 
events  ensure  that  the  carrier  distributions  are  Fermi-like.  Since  the 

energy  of  the  light  quanta  fuo  is  greater  than  the  direct  band  gap  E 

o o 

and  the  indirect  band  gap  E_ , a direct  absorption  event,  such  as  the  one 

V* 

described  above,  gives  an  excess  energy  of  approximately  -E„  to  elec- 
tronic  thermal  agitation.  Thus,  direct  absorption  followed  by  intervalley 
phonon  and  carrier-carrier  scatterings  causes  the  conduction  (valence) 
band  to  contain  a large  number  of  electrons  (holes)  with  a temperature  of 
approximately  1500  K. 

If  only  these  processes  are  considered,  the  Ge  transmission  will  rise 
as  the  optical  pulse  energy  is  increased  because  of  the  saturation  of  the 
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available  optically  coupled  state#  due  to  direct  absorption;  but  no 
dependence  on  pulsewidth  will  result.  However,  intravalley  electron- 
phonon  relaxation  and  plasmon-assisted  electron-hole  recombination,  when 
included,  will  cool  or  heat,  respectively,  the  carrier  distribution.  If 
the  phonon-assisted  cooling  dominates,  more  of  the  electrons  (holes)  will 
be  located  lower  (higher)  in  the  conduction  (valence)  band,  thus  more 
completely  filling  the  states  needed  for  absorption;  the  transmission 
will  further  increase.  However,  if  heating  due  to  plasmon-assisted 
recombination  is  the  stronger,  the  electrons  (holes)  will  be  located 
higher  (lower)  in  the  conduction  (valence)  band,  thus  freeing  the  states 
needed  for  optical  absorption;  the  transmission  will  increase  more  slowly 
as  a function  optical  pulse  energy. 

To  understand  the  pulsewidth  dependence  of  these  two  processes,  con- 
sider two  pulses  of  equal  energy  but  different  width  incident  on  a Ge 
sample.  While  the  longer  pulse  interacts  with  the  sample,  there  is  more 
time  for  both  phonon-assisted  relaxation  and  plasmon-assisted  recombination. 
If  heating  due  to  plasmon  emission  is  dominant,  the  longer  pulse  will  have 
the  higher  distribution  temperature  and,  consequently,  the  lower  trans- 
mission. If  the  phonon-assisted  relaxation  is  stronger,  similar  arguments 
show  the  transmission  will  have  the  opposite  pulsewidth  dependence.  Al- 
though we  have  simplified  the  situation  somewhat,  these  are  the  dominant 
effects.  The  coupling  constants  in  the  paper  by  Elci  et  al.1  were  such 
that  phonon-assisted  cooling  of  the  carrier  distribution  was  small  during 
the  time  the  pulse  was  present  in  the  sample,  thus  plasmon  heating  effects 
determined  the  pulsewidth  dependence. 

Thus,  the  present  data  represents  the  first  substantial  disagreement 
between  experiment  and  a heretofore  successful  model.  (From  our  discussion 
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of  the  pulsewidth  dependence  of  the  model  it  appears  that  for  some  reason 
either  the  plasmon-assisted  recombination  is  not  as  strong  as  we  originally 
calculated  or  that  the  phonon-assisted  relaxation  is  stronger) . Whether 
this  disagreement  is  due  to  the  simplifying  approximations  made  in  the  model, 
to  the  experimental  uncertainties  of  physical  constants  used,  or  to  the 
details  of  the  physical  processes  considered  in  the  model  is  currently 
being  investigated.  Because  the  calculations  of  Elci  et  al.*  have  pro- 
vided a theoretical  description  of  a substantial  body  of  experimental  data, 
we  still  believe  this  model  to  be  substantially  correct.  However,  the 
present  measurements  indicate  the  need  for  further  experimental  and  theo- 
retical study  in  order  to  place  our  understanding  of  these  ultrafast  pro- 
cesses on  a firm  basis. 

We  thank  W.  P.  Latham  and  J.  C.  Matter  for  their  comments  concerning 
the  manuscript. 
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Figure  Captions 


1.  Predicted  single  pulse  transmission  of  a thin  Ge  sample  vs.  incident 
pulse  onerqy  density  (number  of  quanta/interaction  volume)  for  pulse- 
widths  of  2.5,  5,  and  10  psec  for  an  initial  lattice  temperature  of 
77  K (from  Elci  et  ai.1). 

2.  Block  diagram  of  experimental  apparatus  for  producing  pulses  of 
variable  width. 

3.  Experimental  single  pulse  transmission  of  a 7.5pm  - thick  Ge  sample 
vs.  incident  pulse  energy  for  pulsewidths  of  t and  2t  for  a lattice 
temperature  of  105  K.  Corresponding  theoretical  curves  are  shown 
for  pulsewidths  of  2.5,  5,  and  10  psec. 
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ABSTRACT 


Recently,  Elci,  et  al.  developed  a first  principles 
model  that  accounts  for  the  generation  and  transient  be- 
havior of  dense  electron-hole  plasmas  produced  in  germanium 
by  picosecond  optical  pulses.  The  agreement  between  this 
model  and  early  experiments  is  substantial.  However,  recent 
transmission  measurements,  involving  optical  pulses  of  vary- 
ing width,  are  in  disagreement  with  initial  predictions  of 
this  model.  Here,  we  emphasize  the  dependence  of  the  optical 
properties  of  the  Ge  plasma  on  the  electron-phonon  coupling 
constant,  the  broadening  in  the  plasmon-assisted  recombina- 
tion, and  the  energy  band  structure,  and  we  suggest  adjust- 
ments to  the  original  model  that  should  produce  agreement 
with  these  more  recent  experiments.  The  pulsewidth  depend- 
ence in  the  transmission  is  shown  to  be  sensitive  to  the 
relative  strengths  of  the  electron-phonon  relaxation  and 
the  plasmon  recombination  during  the  temporal  evolution  of 
the  hot  electron-hole  plasma. 


I . INTRODUCTION 


I 
I 

I Recently,  Elci,  et  al.  have  proposed  a theoretical 

model  that  describes  the  generation  and  the  temporal  evolu- 
tion of  hot  electron-hole  plasmas  produced  in  the  elemental 
semiconductor  germanium  by  the  absorption  of  intense,  ultra- 
short  optical  (X  - 1.06pm)  pulses.  This  model  has  been  shown 
I to  be  in  agreement  with  the  experimentally-observed  enhanced 

transmission  of  single,  energetic  pulses  through  thin  ger- 
manium samples  as  a function  of  incident  pulse  energy  for 

2 3 

sample  temperatures  of  105  K and  300  K.  ' It  has  also  suc- 
cessfully described  the  transmission  of  a weak  probe  pulse 
through  the  germanium  as  a function  of  time  delay  after  an 
energetic  excitation  pulse  for  different  excitation  energies 
and  for  sample  temperatures  of  105  K and  300  K.3  Subse- 
quently, van  Driel,  et  al.*  have  conducted  further  trans- 
mission studies,  in  which  the  energy  band  gap  of  the  ger- 
manium sample  is  tuned  by  hydrostatic  pressure,  that  corrob- 
orate the  proposed  model. 

However,  Bessey,  et  al.^  have  recently  experimentally 
measured  the  optical  pulsewidth  dependence  of  the  transmission 
of  thin  germanium  samples,  as  suggested  by  Elci,  et  al.6, 
and  l:ave  found  these  measurements  to  be  in  substantial  disa- 
greement with  the  predictions  of  the  model.  In  addition, 

Elci,  et  al.  have  reported7  that,  although  normalized,  rela- 
tive transmissions  are  in  agreement  with  experiment,  absolute 
calculated  transmissions  are  somewhat  higher  than  those 

li 


ill 
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measured.  The  question  thus  arises  as  to  whether  these 
disagreements  with  the  proposed  model  are  due  to  the  assump- 
tions made  in  deriving  the  model,  the  experimental  uncer- 
tainties in  physical  constants  used  in  the  calculations, 
or  the  physical  processes  included  in  the  model. 

We  demonstrate  in  this  paper  that  variations  produced 
in  the  model  by  attempts  to  simplify  the  germanium  energy 
band  structure  and  by  the  experimental  uncertainties  in  the 
electron-phonon  coupling  constants  and  the  strength  of  the 
plasmon-assisted  electron-hole  recombination  are  sufficient 
to  allow  this  disagreement  between  model  and  experiment. 

In  Section  II,  we  begin  with  a brief  review  of  the 
theoretical  model  of  Elci,  et  al.^  Then,  in  Section  III, 
we  demonstrate  that,  in  this  model,  the  absolute  optical 
transmission  is  drastically  affected  by  replacing  the  non- 
parabolic energy  band  structure  of  germanium  by  a parabolic 
band  structure.  In  Section  IV,  we  show  that  the  optical 
pulsewidth  dependence  of  the  phonon-assisted  carrier  relax- 
ation processes  and  the  plasmon-assisted  recombination 
processes  is  opposite  and  that  the  pulsewidth  dependence 
of  the  transmission  can  be  varied  at  will  by  adjusting  the 
relative  strengths  of  the  two  processes. 


II. 


REVIEW  OF  THE  MODEL 


According  to  the  model  proposed  in  Ref.  1,  the  trans- 
mission of  a single  optical  pulse  through  a thin  germanium 
sample  as  a function  of  incident  pulse  energy  can  be  ac- 
counted for  as  follows.  When  an  excitation  pulse  is  inci- 
dent on  a Ge  sample,  the  unreflected  portion  of  the  pulse 
enters  the  sample  where  most  of  it  is  absorbed  by  direct 
transitions,  creating  a large  number  of  electrons  (holes) 

in  the  central  valley  of  the  conduction  (valence)  band. 

- _i  4 

The  electrons  are  rapidly  ( < 10  sec)  scattered  to  the 
conduction  band  side  valleys  by  long  wave  vector  phonons. 
Particle-particle  scattering  events,  which  occur  at  a rate 
faster  than  or  comparable  to  the  direct  absorption  rate, 
ensure  that  the  carrier  distributions  are  Fermi-like  and 
that  both  electron  and  hole  distributions  have  the  same 
temperature,  which  can  be  different  from  the  lattice  temp- 
erature. Since  the  photon  energy  fta)o  is  greater  than  either 
the  direct  energy  gap  EQ  or  the  indirect  band  gap  E^,  such 
a direct  absorption  event  followed  by  phonon-assisted 
scattering  of  an  electron  to  the  side  valleys  results  in 
the  photon  giving  an  excess  energy  of  tiwo-EG  to  thermal 
agitation.  This  excess  energy  results  in  an  initial  dis- 
tribution temperature  (approximately  1500  K for  a lattice 
temperature  of  300  K)  due  to  direct  absorption  that  is 
greater  than  the  lattice  temperature.  If  the  above  processes 
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wore  the  only  ones  considered,  the  single  pulse  optical 
transmission  would  begin  at  its  Beer's  law  value  and  in- 
crease as  a function  of  incident  pulse  energy  because  of 
the  saturation  of  the  optically  coupled  states  as  a result 
of  direct  absorption.  Note  that,  owing  to  the  rapid  scat- 
tering of  electrons  from  the  central  to  side  valleys  of 
the  conduction  band  by  phonons,  any  decrease  in  the  num- 
ber of  states  available  in  the  central  valley  for  direct 
absorption  is  ultimately  determined  by  the  buildup  of  the 
electron  population  in  the  X and  L-side  valleys. 

As  the  carrier  density  increases  because  of  direct 
absorption,  the  plasma  frequency  of  those  carriers  increases. 
At  sufficiently  large  plasma  frequencies,  an  electron  in 
the  T-valley  can  recombine  with  a hole  near  the  top  of  the 
valence  band  via  emission  of  a plasmon.  Normally,  plasmon- 
assisted  recombination  can  occur  only  if  the  plasma  fre- 
quency ii)p  is  larger  than  the  direct  band  gap  frequency  EQ/fi. 
However,  in  our  case,  the  plasma  resonance  is  considerably 
broadened  due  to  the  nature  of  the  direct  absorption  and  sub- 
sequent scattering  between  the  central  and  the  side  conduction 
band  valleys.  Therefore,  plasmon-assisted  recombination  can 
occur  at  plasma  frequencies  lower  than  E^/ft.  These  collec- 
tive plasma  oscillations  have  lifetimes  that  are  short  com- 
pared to  a picosecond.  The  energy  lost  in  the  decay  of  the 
plasmcns  is  rapidly  transferred  to  single  electron  and  hole 
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states  and,  ultimately,  increases  the  temperature  of  the 
carrier  distribution.  Therefore,  as  the  electron  and  hole 
populations  increase,  the  plasmon-assisted  recombination 
rate  increases,  which  limits  the  carrier  number  and  raises 
the  electron-hole  temperature.  Thus,  further  increase  in 
the  sample  transmission  is  prevented  as  the  electrons  (holes) 
are  heated  and  removed  from  the  optically  coupled  states. 

Free  carrier  absorption  events,  although  less  important, 
also  serve  to  increase  the  distribution  temperature  while 
the  pulse  is  incident  on  the  sample. 

Owing  to  direct  absorption  and  plasmon-assisted  recom- 
bination, the  interaction  region  of  the  sample  contains  a 
large  number  of  carriers  with  a high  electron-hole  tempera- 
ture. These  electrons  (holes)  located  high  in  a conduction 
(valence)  band  valley  can  relax  by  emitting  phonons.  The 
effect  of  this  relaxation  is  to  reduce  the  electron-hole 
temperature  and  increase  the  lattice  temperature.  As  the 
distribution  cools  by  means  of  these  intraband  phonon-assisted 
transitions  and  carriers  fill  the  states  needed  for  absorption, 
the  transmission  increases.  The  phonon-assisted  relaxation 
occurs  on  a time  scale  comparable  to  the  pulse  width. 

After  the  passage  of  the  optical  pulse,  the  interaction 
region  of  the  sample  contains  a large  number  of  carriers  with 
a high  distribution  temperature.  The  final  temperature  is 
determined  by  the  number  of  quanta  in  the  optical  pulse  and 


1 If. 

by  the  relative  strengths  of  the  plasmon-assisted  recombina- 
tion and  the  phonon-assisted  relaxation.  The  plasmon- 
assisted  recombination  is  essentially  turned  off  when  the 
excitation  pulse  has  passed  through  the  sample.  As  time 
progresses,  the  distribution  will  continue  to  cool  by  phonon- 
assisted  intravalley  relaxation.  Experimentally,  a weak 
probe  pulse  interrogates  the  evolution  of  the  distribution 
after  the  passage  of  the  first  pulse  and  is  a sensitive 
measure  of  whether  the  optically  coupled  states  are  avail- 
able for  absorption  or  are  occupied.  Immediately  after  the 
passage  of  the  excitation  pulse  the  probe  transmission  is 
small  since  the  electrons  (holes)  are  located  high  (low)  in 
the  bands  because  of  the  high  distribution  temperature, 
leaving  the  states  that  are  optically  coupled  available  for 
direct  absorption.  Later,  as  the  distribution  temperature 
cools  and  carriers  fill  the  states  needed  for  absorption, 
the  transmission  increases. 

The  position  of  the  optically  coupled  states  needed 
for  direct  absorption  depends  on  the  model  band  structure. 

A discussion  of  the  effects  of  the  simplified  band  structure 
on  the  germanium  transmission  is  given  in  the  next  section. 


III.  SIMPLIFICATIONS  TO  THE  BAND  STRUCTURE 


One  of  the  basic  simplifications  in  the  theoretical 
calculations  of  Ref.  1 is  to  replace  the  actual  germanium 
energy  band  structure  with  a parabolic  band  structure  having 
two  degenerate  valence  bands  and  a conduction  band  with  a 
direct  valley  having  energy  gap  Eq  and  ten  equivalent  side 
valleys  with  averaged  energy  gap  Eg.  Actually,  the  central 
valley  in  germanium  is  highly  non-parabolic,  and  the  six 
side  valleys  in  the  (1,0,0)  direction  are  located  higher  in 
energy  than  the  four  side  valleys  in  the  (1,1,1)  direction. 

In  this  section,  we  investigate  numerically  the  effect  of 
assuming  a parabolic  central  conduction  band  valley  on  the 
absolute  value  of  the  optical  transmission.  We  illustrate 
that  such  assumptions  are  sufficient  to  account  for  the  dis- 
crepancy between  the  calculated  and  observed  transmissions 
reported  in  Ref.  1.^ 

If  the  calculated  transmission  is  to  fit  accurately  the 

experimental  data,  the  number  of  available  states  in  the 

conduction  and  valence  bands  must  be  properly  estimated. 

However,  an  enormous  simplication  in  the  calculations  results 

if  the  bands  are  taken  to  be  parabolic.  At  the  band  edge, 

the  central  valley  of  the  conduction  band  can  be  taken  to 

be  parabolic  with  an  effective  mass  m = 0.04m  , where  m 

o e e 

is  the  electron  rest  mass.  However,  the  photon  energy 
is  greater  than  the  direct  band  gap,  and  the  optically  coupled 
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states  are  located,  in  general,  above  the  central  valley 
edge,  where  the  central  valley  is  broader  than  this  ef- 
fective mass  would  predict.  Given  that  one  wishes  to  make 
a parabolic  approximation  to  the  band  structure,  the  ques- 
tion then  arises  as  to  the  "best"  choice  for  the  central 
valley  effective  mass. 

One  solution,  which  has  already  been  mentioned,  is  to 
choose  the  effective  mass  equal  to  the  band  edge  effective 
mass,  mQ  = 0.04me.  However,  this  procedure  results  in  an 
erroneous  value  for  the  linear  absorption  coefficient  since 
it  does  not  accurately  reflect  the  curvature  of  the  bands 
f near  the  optically  coupled  states.  It  also  results  in  an 

erroneous  estimate  of  the  states  available  up  to  the  optically 
coupled  states. 

g 

Another  solution  is  to  equate  the  expression  for  the 
Beer's  law  absorption  coefficient  calculated  assuming  para- 
bolic bands  to  the  experimentally  measured  coefficient  at 
the  frequency  and  to  solve  for  mQ.  This  procedure  ensures 
an  accurate  linear  absorption  coefficient  but  still  results 

( 

in  an  error  in  estimating  the  number  of  states  available  up 
to,  and  including,  the  optically  coupled  states.  This  pro- 
cedure produces  a value  of  m * 0.06m  at  77  K and  m = 0.12m 

o e o e 

at  300  K.  A higher  value  of  mQ  is  obtained  for  the  higher 
temperature,  since  the  band  gap  is  less  and  the  optically 


119 


coupled  states  are  located  higher  in  the  actual  band  struc- 
ture . 

Figure  1 and  Table  I illustrate  the  magnitude  of  the 
error  that  can  be  introduced  by  making  a parabolic  approxi- 
mation to  the  direct  conduction  band  valley.  Figure  1 shows 
the  simplified  germanium  band  structure  used  in  the  calcu- 
lations of  Ref.  1 for  the  two  effective  masses,  m = 0.04m 

o e 

and  mQ  = 0.1me,  at  a temperature  of  300  K.  The  optically 
coupled  states  are  indicated  in  the  figure.  Table  I illus- 
trates the  effect  that  these  two  choices  of  effective  mass 

t 

can  have  on  the  number  of  available  states  in  the  conduction 

band  central  and  side  valleys  up  to  the  optically  coupled 

states.  Note  that,  as  the  central  valley  effective  mass  is 

varied  from  0.04m  to  0.10m  , the  number  of  electrons  required 

e e 

to  fill  the  conduction  band  up  to  and  including  the  optically 
coupled  states  is  not  appreciably  changed;  however,  the  .um- 
ber of  holes  required  to  fill  the  valence  band  is  tripled. 

The  choice  of  the  central  valley  effective  mass  can  have 
a dramatic  effect  on  the  germanium  transmission.  Figure  2 
is  a plot  of  the  single  pulse  transmission  of  a 5.2um-thick 
germanium  sample  as  a function  of  the  incident  number  of 
quanta  for  four  effective  masses.  The  sample  temperature  is 
iOO  K.  The  curves  are  generated  by  using  the  model  presented 
in  Ref.  1 and  varying  the  central  valley  effective  mass,  mQ. 
The  transmission  is  higher  for  the  smaller  effective  mass 
since  a given  number  of  holes  (created  by  direct  absorption 

r 

i 

r 
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of  the  optical  pulse)  will  more  nearly  fill  the  states  avail- 
able in  the  valence  band  up  to  the  optically  coupled  states 
(see  Fig.  1).  Figure  3,  which  shows  the  transmission  of  a 
weak  probe  pulse  as  a function  of  time  delay  after  the  ex- 
citation pulse,  illustrates  the  effect  of  mQ  on  the  sample 
transmission  as  the  carrier  distribution  cools. 

The  "best"  choice  for  mQ  probably  lies  between  the  value 
at  the  band  edge  and  that  calculated  to  give  the  correct  ex- 
perimental linear  absorption  coefficient  at  cjo,  since  the 
former  results  in  underestimating  the  available  valence  band 
states  needed  for  saturation  and  the  latter  overestimates 
this  nur..ber.  Notice  that,  if  a parabolic  approximation  is 
made  to  the  energy  band  structure,  a single  effective  mass 
cannot  be  chosen  that  accurately  predicts  both  the  linear 
absorption  coefficient  and  the  number  of  available  states 
up  to,  and  including  the  optically  coupled  states.  Also  note 

that  the  "best"  choice  for  m at  300  K will  be  different  from 

o 

the  "best"  choice  at  77  K since  the  direct  band  gap  is  narrower 
at  300  K,  and  the  optically  coupled  states  are  located  higher 
(lower)  in  the  conduction  (valence)  band. 

It  is  clear  then  that  assuming  a parabolic  central  con- 
duction band  valley  leads  to  uncertainty  in  the  value  of  the 
absolute  transmission.  Similar  uncertainties  are  introduced 
by  attempts  to  obtain  a simplified  model  for  the  valence 
band  and  for  the  conduction  band  side  valleys. 


IV.  PLASMONS,  PHONONS  AND  PULSEWIDTH  DEPENDENCE 

In  this  section,  we  use  the  model  of  Elci,  et  al.1 
to  investigate  numerically  the  roles  of  the  intravalley 
electron-phonon  relaxation  and  plasmon-assisted  electron- 
hole  recombination  in  determining  the  dependence  of  the 
optical  transmission  on  pulsewidth.  As  we  pointed  out  in 
the  introduction,  the  experimental  measurements5  are  in 
substantial  disagreement  with  the  original  calculations,1 
as  shown  in  Fig.  4.  Note  that  increasing  the  pulsewidth 
lowers  the  theoretical  transmission  for  larger  pulse  inten- 
sities. However,  there  is  little  if  any  pulsewidth  dependence 
in  the  data. 

In  the  model,  intravalley  phonon-assisted  relaxation 
and  plasmon-assisted  recombination  will  cool  and  heat,  respec- 
tively, the  carrier  distribution.  If  the  phonon-assisted 
cooling  dominates,  more  of  the  electrons  (holes)  created  by 
direct  absorption  of  the  optical  pulse  will  be  located  lower 
(higher)  in  the  conduction  (valence)  band.  Thus  the  states 
needed  for  absorption  will  be  more  completely  filled,  and 
the  transmission  will  be  enhanced.  However,  if  heating  due 
to  plasmon-assisted  recombination  is  the  stronger,  the  elec- 
trons (holes)  will  be  located  higher  (lower)  in  the  conduction 
(valence)  band.  The  states  needed  for  optical  absorption 
will  thus  be  freed,  and  this  process  will  tend  to  decrease  the 
transmission  as  a function  of  optical  pulse  energy. 
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To  understand  the  pulsewidth  dependence  of  these  two 
processes,  consider  two  pulses  of  equal  energy  but  different 
width  incident  on  a germanium  sample.  While  the  longer  pulse 
interacts  with  the  sample  there  is  more  time  for  both  phonon- 
assisted  relaxation  and  plasmon-assisted  recombination.  If 
heating  due  to  the  plasmon  emission  is  dominant,  the  longer 
pulse  will  have  a higher  distribution  temperature  and,  con- 
sequently, a lower  transmission.  This  trend  is  illustrated 
in  Fig.  5.  This  figure  is  a plot  of  the  single  pulse  trans- 
mission of  germanium  as  a function  of  incident  pulse  energy 
for  the  three  pulse-widths:  2.5,  5,  and  10  psec.  The  phonon- 
assisted  relaxation  terms  are  neglected  in  this  figure.  If, 
however,  the  phonon-assisted  relaxation  is  stronger,  similar 
arguments  show  that  the  transmission  will  have  the  opposite 
pulsewidth  dependence,  as  presented  in  Fig.  6.  Figure  6 is 
identical  to  Fig.  5 except  that  the  plasmon-assisted  recom- 
bination terms  have  been  neglected  while  retaining  the  phonon 
relaxation  terms. 

The  actual  pulsewidth  dependence  is  determined  by  which 
of  these  two  processes  is  stronger.  The  coupling  constants 
in  the  paper  by  Elci,  et  al.*  were  such  that  phonon-assisted 
cooling  of  the  carrier  distribution  was  small  during  the  time 
the  pulse  was  incident  on  the  sample,  thus  plasmon  heating 
effects  determined  the  pulsewidth  dependence.  The  experimental 
data  obtained  by  Bessey,  et  al.5  indicate  that  either  the  plas- 
mon effects  were  overestimated  in  the  original  calculations 
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or  the  phonon-assisted  relaxation  was  underestimated.  Indeed, 
an  uncertainty  does  exist  in  both  processes. 

Consider  for  example,  the  accuracy  with  which  the  electron- 
optical  phonon  coupling  constant  Qq  is  known.  Table  2 contains 
a list  of  several  experimentally  and  theoretically  determined 
values  for  Qq.  Evidently,  there  is  a substantial  uncertainty 
in  Qq.  This  uncertainty  has  a dramatic  effect  on  the  evolu- 
tion of  the  optical  transmission  both  while  the  excitation 
pulse  is  incident  on  the  sample  and  after  its  passage.  Figure  7 
illustrates  the  change  that  occurs  in  the  transmission  as  a 
function  of  incident  pulse  energy  when  Qq  is  varied.  The 
plasmon  terms  have  the  same  strength  as  in  Ref.  1.  For  the 
lower  values  of  Qq,  there  is  negligible  cooling  during  the 
period  the  pulse  interacts  with  the  sample;  the  final  values 
of  transmission  are  determined  primarily  by  direct  absorption 
followed  by  plasmon  heating  of  the  carrier  distribution.  For 
this  range  of  values,  the  dependence  on  Qq  is  small.  However, 
for  the  larger  values  of  Qq  (for  example,  Qq  = 2xl0~3  ^2-)  * 
the  phonon-assisted  cooling  competes  with  the  plasmon-assisted 
heating  causing  the  electrons  (holes)  to  be  located  lower 
(higher)  in  the  conduction  (valence)  band,  which  results  in 
a higher  transmission. 

The  uncertainty  in  Qq  also  affects  the  temporal  evolu- 
tion of  the  germanium  transmission  after  the  passage  of  the 
excitation  pulse  as  shown  in  Fig.  8.  The  rise  in  these 
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transmission  curves  is  directly  determined  by  the  time  it 
takes  the  electrons  (holes)  to  relax  to  the  bottom  (top)  of 
the  conduction  (valence)  band  and  more  completely  clog  the 
states  needed  for  absorption.  In  fact,  interrogation  of 
the  temporal  evolution  of  the  transmission  of  the  sample 
by  a weak  probe  pulse  would  represent  a sensitive  method  for 
determining  Qq,  if  all  other  parameters  associated  with  the 
model  were  precisely  known. 

There  is  also  an  uncertainty  in  the  plasmon-assisted 
recombination  that  is  a result  of  the  approximations  made 
in  the  calculation.  In  Ref.  1,  the  plasmon  recombination 
rate  was  calculated  using  the  Drude  form  of  the  dielectric 
function.  Thus,  the  calculation  of  the  plasmon  rate  must  be 
regarded  as  phenomenological  with  the  plasmon  broadening 
being  an  estimated  constant  in  this  model.  The  plasmon 
broadening  was  estimated  to  be  proportional  to  ten  times  the 
rate  for  scattering  of  an  electron  from  the  T-valley  to  the 
(111) -valley  accompanied  by  the  emission  of  a phonon. 

One  convenient  way  of  demonstrating  the  effect  of  the  un- 
certainty in  the  plasmon  recombination  rate  is  to  vary  the 
plasmon  broadening  by  varying  the  phonon-assisted  intervalley 
scattering  time,  tq,  which  is  inversely  proportional  to  the 
broadening.  Figure  9 illustrates  the  change  in  the  trans- 
mission as  a function  of  incident  pulse  energy  when  tq  is 
varied.  Increasing  tq  markedly  increases  the  transmission  for 
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larger  incident  energies.  For  larger  tq,  the  plasmon 
broadening  is  smaller  and  plasmon  recombinations  are  less 
effective  in  heating  the  carrier  distribution.  Thus,  the 
electrons  (holes)  are  located  lower  (higher)  in  the  conduc- 
tion (valence)  band;  the  transmission  is  higher. 

If  either  the  electron-optical  phonon  coupling  constant 
or  the  plasmon  broadening  were  known  with  certainty,  the 
other  constant  could  be  estimated  within  this  model  by  varying 
the  unknown  constant  to  obtain  the  proper  transmission  and 
pulsewidth  dependence.  However,  the  uncertainties  in  the 
phonon  relaxation  and  plasmon  recombination  rates  make  an 
accurate  evaluation  of  either  constant  impossible  within  the 


present  model. 


CONCLUSION 


In  the  previous  sections,  we  have  illustrated  the 
effect  of  the  simplified  band  structure,  the  phonon-assisted 
relaxation,  and  the  plasmon-assisted  recombination  on  the 
transmission  of  picosecond  pulses  through  germanium.  In 
particular,  we  have  shown  that  the  pulsewidth  dependence  of 
the  plasmon  effects  is  opposite  to  that  of  the  phonon- 
assisted  relaxation.  Because  of  the  uncertainties  in  the 
electron-optical  phonon  coupling  constant  and  the  plasmon 
broadening,  the  deviation  of  the  experimentally-measured 
pulsewidth  dependence  of  the  germanium  transmission  from 
the  initially  predicted  pulsewidth  dependence  can  be  accounted 
for  within  the  framework  of  the  model,  if  the  proper  balance 
is  maintained  between  plasmon  recombination  and  phonon  re- 
laxation during  the  passage  of  the  pulse  through  the  sample. 
Since  this  pulsewidth  dependence  can  be  adjusted  almost  at 
will  because  of  these  uncertainties,  the  pulsewidth  experi- 
ments do  not  provide  a definitive  test  of  the  model.  We 
believe  the  model  to  be  essentially  correct  since  it  agrees 
with  a substantial  body  of  experimental  data;  however,  further 
experimental  and  theoretical  study  is  required  to  completely 
determine  the  model  and  its  physical  constants. 
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LIST  OF  TABLES 

Table  1 The  optically  coupled  energies  and  number  of 
electrons  or  holes  up  to  and  including  the 
optically  coupled  energies  for  several  values 
of  the  central  valley  effective  mass  using  the 
model  band  structure  at  300  K. 

Table  II  Several  experimental  and  theoretical  values  of 
the  electron-optical  phonon  coupling  constant. 
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Fig.  7 The  single  pulse  transmission  of  germanium  as 

a function  of  incident  pulse  energy  for  several 
values  of  the  electron-optical  phonon  coupling 
constant,  Qq. 

Fig.  8 The  instantaneous  probe  transmission  as  a function 
of  time  delay  after  an  excitation  pulse  of  photon 
density  ST  (number  of  quanta/interaction  volume) 
for  several  values  of  the  electron-optical  phonon 
coupling  constant,  Qq. 

Fig  9 The  single  pulse  transmission  of  germanium  as  a 

function  of  incident  pulse  energy  for  several 

-14 

values  of  the  plasmon  broadening;  tq  (~10  sec.) 
is  the  intervalley  phonon-assisted  scattering  time 
used  in  the  broadening  in  Ref.  1. 
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ABSTRACT 

The  conventional  perturbation  theory  in  quantum 
mechanics  for  the  interaction  of  electromagnetic  radiation 
with  matter,  which  is  based  on  the  time-dependent  vector 
and  scalar  potentials,  is  shown  not  to  be  gauge  invar  iant- 
The  perturbation  problem  is  reformulated,  using  the  approach 
of  K.-H.  Yang,  in  a gauge- invar iant  way,  in  which  the  trans- 
itions between  states  are  induced  by  the  quantum  mechanical 
power  operator.  In  order  to  use  an  unperturbed  Hamiltonian 
which  does  not  involve  the  electromagnetic  field,  the  gauge 
is  chosen  to  give  the  electric  dipole  interaction.  This 
choice  can  be  made  if  the  magnetic  field  is  negligible  and 
the  electric  field  is  slowly  varying  over  atomic  dimensions. 
As  an  example,  the  Kramers-Heisenberg  formula  is  derived 
from  the  gauge  invariant  electric  dipole  interaction. 


I . INTRODUCTION 


Over  half  a century  has  elapsed  since  the  discovery 

of  quantum  mechanics  and  the  first  treatment  of  radiation 

problems  in  the  new  quantum  theory.^  Although  most  quantum 

mechanics  texts  treat  the  interaction  of  radiation  and 

matter,  few  books  do  more  than  mention  gauge  transformations 

2-5 

and  relegate  them  to  the  problems.  The  necessity  of  doing 

calculations  in  a gauge  invariant  way  is  often  not  emphasized. 
In  fact,  many  books  proceed  to  expand  out  the  terms  involving 

6 _ 9 

the  vector  potential,  and  treat  them  as  perturbations. 

Other  books  use  the  electric  dipole  interaction  from  the  start, 

without  necessarily  relating  it  to  the  vector  potential . ' ** 

A few  books  do  point  out  that  the  electric  dipole  interaction 

is  related  to  the  vector  potential  interaction  by  a gauge 
12  — 16 

transformation.  The  perturbing  question  arises,  however, 

as  to  which  interaction  to  use  - the  one  involving  the  vector 
potential  or  the  electric  dipole  interaction  - or  is  it 
merely  a matter  of  taste?  It  is  the  purpose  of  this  paper 
to  show  that  it  does  make  a difference  and  the  electric  dipole 
interaction  should  be  used  when  it  can  be.  It  may  be  necessary 
to  use  higher  order  electric  multipoles  and  magnetic  multipoles 
as  well. 
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In  order  to  substantiate  this  claim,  a gauge  invariant 

time-dependent  perturbation  theory  is  developed.^"7  Since 

the  Schrodinger  equation  can  not  usually  be  solved  exactly, 

1 8 

the  interaction  between  the  radiation  and  the  atom  is 
treated  by  perturbation  theory.  The  most  convenient  and 
natural  choice  for  the  unperturbed  Hamiltonian  is  the 
Hamiltonian  for  the  atom  in  the  absence  of  the  radiation 
field.  As  we  will  show,  this  choice  of  the  unperturbed 
Hamiltonian  implies  that  the  electric  dipole  interaction 
is  the  correct  gauge-invariant  interaction  to  use  in  these 
problems,  when  the  magnetic  multipole  and  higher-order 
electric  multipole  interactions  are  negligible. 

The  reason  the  electric  dipole  interaction  for  these 
problems  plays  such  a fundamental  role  is  that  it  is  the 
only  interaction  for  which  the  atomic  Hamiltonian  and  the 
electromagnetic  field  are  decoupled  in  a gauge-invariant 
way.  Thus  the  expansion  coefficients  of  the  wave  function 
in  terms  of  the  unperturbed  atomic  states  can  be  interpreted 
as  true,  gauge-independent  probability  amplitudes  for  finding 
the  atomic  system  in  an  unperturbed  state.  For  another 
choice  of  gauge,  the  gauge-invariant  unperturbed  Hamiltonian 
is  shown  to  couple  the  atom  and  the  electromagnetic  field. 

If  the  Hamiltonian  is  broken  up  in  a way  that  does  not 
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preserve  gauge  invariance,  the  expansion  coefficients  for 
the  wave  function  are  gauge  dependent  and  cannot  be  inter- 
preted as  probability  amplitudes.  The  fundamental  significance 

19  . 

of  the  electric  dipole  interaction  was  noted  by  Lamb  in 
1552  when  he  remarked  that  another  choice  of  gauge  gave  results 
that  were  not  in  agreement  with  experiment. 

The  subject  of  gauge  transformations  in  quantum  mechanics 
has  been  recently  discussed  in  an  excellent  article  by  Yang,^7 
who  emphasizes  the  importance  of  doing  calculations  in  a gauge 
invariant  way.  Some  of  his  points  are  reviewed  here  for  com- 
pleteness, and  to  make  the  fundamental  ideas  available  to  a 
wider  audience.  In  particular,  Yang  uses  the  correspondence 
principle  to  show  that  the  gauge-invariant  energy  operator 
for  the  system  is  not  the  Hamiltonian,  which  is  gauge  dependent 
but  the  Hamiltonian  without  the  time-dependent  scalar  potential 
From  this  fundamental  point,  a gauge-invariant  time-dependent 
perturbation  theory  is  developed,  in  which  the  transitions 
between  states  is  induced  by  the  gauge-invariant  quantum 
mechanical  power  operator.  We  have  extended  his  treatment  by 
writing  the  power  operator  in  terms  of  the  electric  dipole 
energy  and  derivatives  of  the  electric  field,  which  shows 
the  naturalness  and  convenience  of  the  electric  dipole  ap- 
proximation. The  gauge  invariant  perturbation  theory  is  then 

applied  to  the  scattering  of  light  by  an  atom  to  derive  the 

20 

Kramers-Heisenberg  formula.  These  ideas  can  be  used  to 
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supplement  the  usual  textbook  treatment  of  radiation  theory 

in  quantum  mechanics  and  quantum  optics. 

The  role  of  gauge  transformations  in  quantum  theory 

21  22 

was  apparently  first  pointed  out  by  Fock.  Weyl  developed 

a unified  theory  of  electromagnetism  and  matter  based  on  a 

principle  of  gauge  invariance . His  theory  was  the  harbinger 

23 

of  the  modern  gauge  theories  of  particle  physics.  The 
principle  of  gauge  invariance  states  that  the  Schrodinger 
equation  is  form  invariant  under  simultaneous  gauge  trans- 
formation of  the  wave  function  and  the  vector  and  scalar 

, 24 

potentials.  Subsequently,  Pauli  called  the  gauge  trans- 
formation of  the  wave  function  a "gauge  transformation  of 
the  first  kind,"  while  the  usual  gauge  transformation  of  the 
electromagnetic  potentials  he  called  a "gauge  transformation 
of  the  second  kind." 

The  goal  of  developing  a completely  gauge  invariant  formu- 
lation of  quantum  electrodynamics  has  been  pursued  by  many 

2 5 26 

workers,  and  has  been  beset  with  many  technical  difficulties. 

A resurgence  of  interest  in  the  role  of  the  electromagnetic 

27 

potentials  in  quantum  theory  occurred  when  Aharonov  and  Bohm 
suggested  that  they  had  some  significance  apart  from  the  electro- 
magnetic fields  they  described.  The  development  of  the  laser 
also  stimulated  interest  in  the  interaction  of  the  electro- 
magnetic field  with  atoms,  molecules,  and  solids.  The  gauge 
invariant  Hamiltonian  for  a system  of  nonrelativistic  particles 
interacting  with  the  radiation  field  was  derived  by  Power  and 
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Zienau  using  a canonical  transformation.  This  method 

29 

has  been  used  and  extended  by  many  others.  A good  review 
of  this  work  has  recently  been  given  by  Woolley."*® 

In  Sec. II  the  conventional  perturbation  treatment 
of  radiation  is  reviewed  to  point  out  some  common  misconceptions. 
In  Sec.  Ill  the  form  invariance  of  the  Schrodinger  equation 
under  gauge  transformations  is  shown.  The  gauge  invariance  of 
expectation  values  for  observables  and  their  equations  of 
motion  is  discussed  in  Sec.  IV,  which  leads  to  the  gauge  in- 
variant energy  operator.  The  formulation  of  a gauge-invariant 
time-dependent  perturbation  theory  is  given  in  Sec.  V . Then 
in  Sec.  VI  the  electric  dipole  interaction  is  derived  by  making 
a particular  choice  of  gauge  and  neglecting  magnetic  and  higher- 
order  electric  multipoles.  In  Sec.  VII  these  ideas  are  applied 
to  derive  the  Kramers-Heisenberg  formula  for  the  scattering 
of  light  by  an  atom.  Finally  the  paper  is  summarized  and  con- 
clusions given  in  Sec  VIII. 

II.  CONVENTIONAL  APPROACH 

The  conventional  treatment  of  radiation  problems  in  quantum 
6 — 8 

mechanics  is  reviewed  here  briefly  to  point  out  its  inade- 
qacies  and  some  common  misconceptions.  The  main  objection  is 
that  the  conventional  perturbation  theory  is  not  gauge  invariant. 
This  fact  leads  to  possible  spurious  interactions  and  expansion 
coefficients  that  are  not  probability  amplitudes.  For  an  un- 
perturbed Hamiltonian  which  does  not  depend  on  the  electromagnetic 
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field,  the  gauge  should  be  choBen  such  that  the  electric 
dipole  interaction  is  the  perturbation,  which  is  discussed 
in  Sec.  VI.  This  interaction  is  not  a weak  field  ap- 
proximation to  the  conventional  treatment,  but  is  the 
correct  gauge  invariant  interaction  when  the  magnetic  field 
is  negligible  and  the  electric  field  is  slowly  varying  over 
the  atom. 

The  Hamiltonian  for  a particle  of  mass  m,  charge  q,  and 
momentum  p = -i-ftV  in  an  electromagnetic  field  described  by 
the  time  varying  vector  potential  A and  scalar  potential  AQ 
is 

2 

| H = H(A,Aq)  = (P  - qA/c)  /2m  + V + qAQ  , (2.1) 

where  V is  an  external  static  potential  energy.  The  quadratic 
term  in  Eq.  (2.1)  can  be  expanded  to  give 

H(A,Aq)  = Hq  - (q/2mc) (A-p  + p-A)  + (q2/2mc2)A2  + qAQ,  (2.2) 

where  Hq  is  the  Hamiltonian  H(0,0)  in  the  absence  of  an  electro- 
magnetic *ield 

Hq  = p2/2m  4 V . (2.3) 


The  Coulomb  gauge 


V-A  = 0 


(2.4) 
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is  often  used  in  Eq.  (2.2)  so  that  p*A  * A*p.  The 
conventional  procedure  is  to  use  Hq  as  the  unperturbed 
Hamiltonian,  and  the  remainder  in  Eq.  (2.2)  as  the  pertur- 
bation. However,  this  approach  is  not  gauge  invariant. 

The  coefficients  in  the  expansion  of  the  wave  function  i); 
in  terms  of  the  eigenstates  of  H0  are  gauge  dependent,  and 
can  therefore  not  be  interpreted  as  probability  amplitudes. 
Appendix  A gives  a discussion  of  the  lack  of  gauge  invariance 
of  the  expansion  coefficients  in  the  conventional  perturbation 
theory.  Even  the  Hamiltonian  in  Eq.  (2.1)  or  (2.2)  is  not 
gauge  invariant,  which  is  discussed  in  Sec.  IV. 

That  there  are  difficulties  with  the  conventional  per- 
turbation theory  can  be  seen  by  considering  the  case  in  which 
there  is  no  electromagnetic  field.  The  situation  in  which 
the  electric  field  E = 0 and  the  magnetic  induction  B = 0 
can  be  described  by  the  potentials  A **  0 and  Aq  = 0 . Then 
the  perturbation  in  Eq.  (2.2)  is  zero,  as  it  should  be. 
However,  as  we  show  in  Sec.  Ill,  the  same  situation  of  no 
electromagnetic  field  can  also  be  described  by  the  vector 
potential  A * VA  and  the  scalar  potential  Aq  =*  -3A/3(ct)f 
where  A is  an  arbitrary  function  of  space  and  time.  Then 
the  perturbation  in  Eq.  (2.2)  is  not  zero,  but  is  a spurious 
interaction.  When  treated  by  perturbation  theory,  this 
interaction  would  give  incorrect  results.  It  is  horrifying 
to  think  of  the  time  that  could  be  spent  in  calculating  the 
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effects  of  this  spurious  interaction  to  nth  order.  This 
example  should  emphasize  the  importance  of  doing  perturba- 
tion theory  in  a gauge  invariant  way. 

As  we  will  see  in  Sec.  VI,  the  gauge  can  be  chosen, 
in  problems  where  the  magnetic  field  is  negligible  and  the 
wavelength  of  the  radiation  is  long  compared  to  atomic  di- 
mensions, such  that  Eq.  (2.1)  becomes 

H = Hq  - qE (0 , t ) • r , (2.5) 

where  E(0,t)  is  the  electric  field  at  the  origin  (chosen  to 
be  the  center  of  the  atom) . The  electric  dipole  interaction 
in  Eq.  (2.5)  is  obviously  gauge  invariant,  since  it  does  not 
depend  on  the  potentials.  It  is  also  the  correct  perturbation 
to  use  in  conjunction  with  Hq. 

It  has  sometimes  been  thought  that  the  Hamiltonian  in 

Eq.  (2.5)  is  only  an  approximation  to  Eq.  (2.2)  that  is  ap- 

2 31 

plicable  in  the  weak  field  case  when  the  A term  is  negligible, 
because  the  electric  dipole  interaction  is  linear  in  the 
electric  *ield.  This  erroneous  conclusion  stems  from  the 
similarity  in  form  between  the  -qE*r  term  and  the  -(q/mc)A*p 
term  in  Eq.  (2.2)  when  the  Coulomb  gauge  in  Eq.  (2.4)  is 
used . 

The  matrix  elements  of  the  two  interaction  terms  are  not 
equal,  but  they  are,  however,  related.  If,  for  example,  Aq 

« 
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is  chosen  to  be  zero  (as  it  can  be  in  the  absence  of  charge 
density) , then  the  electric  field  is 

E = - 9A/3 (ct ) for  A = 0 . (2.6) 

~ o 

For  a harmonically  varying  field  of  frequency  u>,  Eq.  (2.6) 
implies  that 

E ( 0 , t)  » iwA(0,t)/c  . (2.7) 

The  commutation  relations  between  the  components  of  r 
and  p , 

lxi'  pj]  = 6ij'  (2.8) 

where  i = 1,2,3,  implies  that  the  commutator  between  Hq  in 
Eq.  (2.3)  and  r is 

[H-  r ] = - iftp/m  (2.9) 

O'*  * 

The  eigenstates  <J>n  of  Hq  satisfy  the  eigenvalue  problem 

H * “ en  ♦„  » (2.10) 

on  n n 

where  e^  is  the  eigenenergy.  When  matrix  elements  of  Eq.  (2.9) 

are  taken  between  the  eigenstates  <t>  the  result  is 

n, 

<*>  (r)  = -i(p)  /m  , 

nm  ~ nm  nm 


(2.11) 
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where  w is  defined  as 
run 


-h  iii 


nm 


= e -e 
n m 


(2.12) 


and  the  matrix  elements  are  defined  as 

(r)  = <<t>n  I f <*>_>  • (2.13) 

nm  n m 

The  matrix  elements  of  the  electric  dipole  interaction  are 
thus 

-q  E ( 0 , t ) • (r)  ^ = (u)/<JJnm)  (_cl/mc)A  (0 , t)  . (p)  ^ (2.14) 

when  Eqs.  (2.7)  and  (2.11)  are  used. 

If  the  resonance  condition. 


■hw  = e -e  (2.15) 

n m 

is  satisfied,  it  appears  that  Eq.  (2.5)  is  equivalent  to 

2 32 

Eq.  (2.2)  without  the  A and  Aq  terms.  However,  for  virtual 
processes  in  perturbation  theory,  the  resonance  condition  in 
Eq.  (2.15)  is  not  satisfied.  This  point  will  be  discussed  in 
more  detail  in  Sec.  VII. 

Since  the  electric  dipole  interaction  in  Eq.  (2.5)  and 
the  -(q/mc)  A'p  interaction  in  Eq.  (2.2)  give  matrix  elements 
differing  by  a factor  of  w/w  in  Eq.  (2.14),  it  can  be 
decided  experimentally  which  interaction  to  use.  Lamb33  has 
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shown  experimentally  that  the  electric  dipole  interaction 
of  Eq.  (2.5)  is  the  proper  interaction. 

We  shall  show  in  Sec.  VI  on  theoretical  grounds  that 
the  electric  dipole  interaction  is  not  a weak  field  approxi- 
mation to  the  full  Hamiltonian  in  Eq.  (2.2),  but  is  fully 
equivalent  to  Eq.  (2.1)  under  the  approximations  that  (a) 
magnetic  effects  are  negligible  (i.e.,  the  magnetic  induction 
B does  not  appear)  and  (b)  the  field  E(r,t)  is  slowly  varying 
over  the  atom  (i.e.,  the  dipole  approximation  is  valid). 
Furthermore,  Eq.  (2.5)  is  the  correct  Hamiltonian  to  use  when 
Hq  is  taken  to  be  the  unperturbed  Hamiltonian.  It  is  inter- 
esting to  note  that  many  early  workers  in  quantum  theory  did 

20  34 

use  the  electric  dipole  approximation.  ' In  the  next 
section  the  form  invariance  of  the  Schrodinger  equation  under 
gauge  transformations  is  discussed. 

III.  GAUGE  TRANSFORMATIONS 

The  Schrodinger  equation  for  a particle  in  an  electro- 
magnetic field  described  by  the  vector  potential  A and  the 
scalar  potential  Aq  is  shown  in  this  section  to  be  form  in- 
variant if  the  wave  function  and  the  potentials  are  trans- 
formed in  an  appropriate  way.  The  electromagnetic  field  is 
unchanged  if  the  potentials  undergo  a gauge  transformation, 
while  the  probability  density  ♦ ^ is  unchanged  if  ^ is  multi- 
plied by  a phase  factor.  These  transformations  were  shown 
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by  Fock  as  early  as  1927  to  be  intimately  related.  Dis- 
cussions of  the  form  invariance  of  the  Schrodinger  equation 
under  gauge  transformations  are  available  in  several  books,35,37 
but  the  subject  is  reviewed  here  because  of  its  importance  for 
subsequent  sections. 

The  Schrodinger  equation  for  a particle  of  mass  m in  an 
electromagnetic  field  is 

(l/2m)  (p-q£/c)2  g/  + ViJ>  + qA  g<  = ifi  3i^/3t  (3.1) 

»•  o 

when  the  Hamiltonian  of  Eq.  (2.1)  is  used.  A new  wavefunction 
' can  be  obtained  by  multiplying  by  a phase  factor, 

g> ' = exp  ( i qA/cti)g>  , (3.2) 

where  A = A(r,t)  is  a real  function.  Equation  (3.2)  is  called 

2 4 

a local  gauge  transformation  of  the  first  kind  after  Pauli. 

If  Eq.  (3.2)  is  substituted  into  Eq.  (3.1)  the  Schrodinger 
equation  becomes 

( l/2m)  (p-qA'/c)2  g, ' + Vg, ' + qA^g,  ’ = i*3g»'/3t  , (3.3) 

which  has  the  same  mathematical  form  as  Eq.  (3.1).  In  order 
to  preserve  the  form,  the  new  vector  potential  A'  is  defined 
as 


A' 


A + VA 


9 


(3.4) 
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and  the  new  scalar  potential  A^  as 

A’  = A - 3A/3 (ct)  . (3.5) 

o o 

Equations  (3.4)  and  (3.5)  are  the  customary  gauge  trans- 

38 

formations  of  electromagnetic  theory,  which  are  called 

24 

gauge  transformations  of  the  second  kind  by  Pauli. 

The  nonuniqueness  of  the  vector  and  scalar  potentials 

is  due  to  the  way  in  which  the  observed  electric  and  magnetic 

38 

fields  are  expressed  in  terms  of  the  potentials.  The  mag- 
netic induction 

B = V*A  , (3.6) 

is  invariant  under  Eq.  (3.4).  Likewise,  the  electric  field 

E = -VAq  - 3A/3 (ct)  , (3.7) 

is  also  invariant  under  Eqs.  (3.4)  and  (3.5).  Therefore  any 
choice  of  gauge  function  A can  be  used  and  the  same  electric 
and  magnetic  fields  will  be  obtained.  However,  the  wave 
function  is  changed  according  to  Eq.  (3.2). 

Thus  the  Schrodinger  equation  is  form  invariant  under  the 
local  gauge  transformation  of  the  first  kind  in  Eq.  (3.2)  if 
the  scalar  and  vector  potentials  transform  according  to  the 
customary  gauge  transformations  of  the  second  kind  in  Eqs. 

(3.4)  and  (3.5).  In  the  next  section  the  gauge  invariance  of 
expectation  values  for  observables  is  discussed. 
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IV.  GAUGE  INVARIANCE  OF  EXPECTATION  VALUES 

The  expectation  value  of  a physical  observable  should 

not  depend  on  the  choice  of  gauge.  If  it  did,  then  its 

value  would  be  arbitrary,  since  the  gauge  function  A in 

Sec.  Ill  is  arbitrary.  In  this  section,  the  consequences 

of  this  reasonable  requirement  are  investigated.  Following 
17 

Yang,  the  correspondence  principle  is  used  to  determine 
the  proper  operators  for  observables. 

If  0 is  a Hermitian  operator  which  describes  an  observable, 
then  its  expectation  value  in  the  state  ^ must  be  equal  to  its 
expectation  value  in  the  state  ' in  Eq.  (3.2), 

<4/ 1 e«*»>  = <\p>  | 0 ' ' > , (4.1) 

in  order  to  ensure  gauge  invariance.  The  gauge  transformed 
operator  0'  is  defined  as 

0'  = exp(iqA/-hc)  0 exp  (-iqA/lic)  , (4.2) 

which  will  be  called  a gauge  transformation  of  the  first  kind 
on  the  operator.  A prime  denotes  a gauge  transformation  of 
the  first  kind  in  Eq.  (4.2)  on  all  operators  except  the  vector 
potential  A and  the  scalar  potential  Aq,  for  which  a prime 
denotes  a gauge  transformation  of  the  second  kind  defined  in 
Eqs.  (3.4)  and  (3.5),  respectively.  Under  a gauge  transformation 


of  the  first  kind  in  Eq.  (4.2)  A and  Aq  are  unchanged,  so 
the  prime  on  them  is  reserved  for  the  gauge  transformation 
of  the  second  kind. 

The  Schrodinger  equation  in  Eq.  (3.1)  is  form  invariant 
under  the  gauge  transformation  on  the  wave  function  in  Eq. 

(3.2),  when  Eqs.  (3.4)  and  (3.5)  are  used  for  the  new  poten- 
tials. The  expectation  values  of  all  operators  are  inde- 
pendent of  gauge  transformations  of  the  first  kind  when  Eqs. 
(4.1)  and  (4.2)  are  used.  However,  not  all  operators  are 
form  invariant  under  the  gauge  transformation  of  the  first 
kind  in  Eq.  (4.2).  For  example,  the  momentum  operator 
p = -irfiV  is  not  form  invariant  under  Eq.  (4.2),  but  the  com- 

V 
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bination  (p  - qA/c)  is.  A gauge  transformation  of  the 
first  kind  in  Eq.  (4.2)  on  (p  - qA/c)  gives 

(p  - qA/c)'  = exp(iqVhc)  (p  - qA/c)  exp(-iqA/ftc) 

= (p  - qVA/c)  - qA/c 

= (p  - qA'/c) , (4.3) 

where  Eq.  (3.4)  has  been  used  for  A'.  The  calculation  here 
has  been  given  explicitly  to  show  how  the  vector  potential 
becomes  transformed  by  a gauge  transformation  of  the  second 
kind  as  a consequence  of  the  gauge  transformation  of  the  first 
kind  on  p.  Likewise,  the  operator  rh  y 3t  is  not  form  invariant 
under  the  gauge  transformation  in  Eq.  (4.2),  but  the  combination 
(ihyst  - qAQ)  is.  A gauge  transformation  of  the  first  kind  on 


. 


U.0 

(ift3/3t-  qA  ) gives 
o 

( ifi  3 / 3 1 - qA0)  ' = exp  (iqA/ftc)  (i*3/3t  - qAQ)  exp  ( -iq A/f»c ) 

= (ift3/3t  + q3A/3 (ct) ) -qAQ 

= ifi  3 / 3 1 - qA^  , (4.4) 

where  Eq . (3.5)  has  Leen  used  for  A^.  The  scalar  potential 

becomes  transformed  by  a gauge  transformation  of  the  second 
kind  as  a consequence  of  the  gauge  transformation  of  the  first 
kind  on  rfi9/3t.  We  emphasize  again  that  the  prime  on  the  vector 
potential  A*  and  the  scalar  potential  A^  denotes  the  gauge 
transformation  of  the  second  kind  given  in  Eqs . (3.4)  and  (3.5), 

respectively.  The  prime  on  all  other  operators  9'  denotes  the 
gauge  transformation  of  the  first  kind  in  Eq.  (4.2).  The 
operator  9 may  involve  A or  Aq  as  in  Eqs.  (4.3)  or  (4.4),  but 
the  prime  on  the  operator  9 still  denotes  the  operation  shown 
in  Eq.  (4.2). 

Another  operator  of  particular  interest  that  is  not  form 
invariant  under  Eq.  (4.2)  is  the  Hamiltonian  H(A,Aq)  in  Eq. 
(2.1).  The  Hamiltonian  transforms  under  Eq.  (4.2)  as 

H'(A,Aq)  * exp(iqA/tic)  H(A,Aq)  exp  ( -iqA/fic) 

- H(A\Aq) 

* H (A ' ,A^)  + q3A/3 (ct)  , (4.5) 


where  Eqs.  (4.3),  (3.4),  and  (3.5)  have  been  used.  Al* hough 
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the  time  evolution  of  the  system  with  wave  function  is 
described  by  the  Hamiltonian  H(A,Aq)  through  the  Schrodinger 
equation,  it  need  not  represent  the  instantaneous  energy  oper- 
ator. In  fact  Eq.  (4.5)  shows  that  the  total  Hamiltonian  H(A,Aq) 
can  not  represent  the  energy  operator,  since  its  expectation 
value  is  not  gauge  invariant,  i.e., 

<*|H(A,A  ) *>  = <iT  | H * (A,Ao)*’> 

= <<*>’  |H  (A',Aq)  ♦ ' > 

y <*'  |H(A’  > . (4.6) 

The  question  thus  arises  as  to  the  proper  choice  for  the  energy 
operator  and  other  operators  representing  physical  observables. 

This  question  is  answered  in  the  remainder  of  this  section. 

An  operator  representing  a physical  observable  must,  of 
course,  have  an  expectation  value  which  is  independent  of  the 
choice  of  gauge.  Yang17  has  used  the  correspondence  principle 
to  select  the  appropriate  operators  to  represent  physical 
observables.  According  to  him,  an  operator  represents  a physical 
quantity  with  a classical  analogue  only  if  the  equation  of  motion 
for  the  expectation  value  of  the  operator  is  of  the  same  form 
as  the  equation  of  motion  for  the  corresponding  classical  quantity. 

The  equation  of  motion  for  the  expectation  value  of  an  oper- 
ator 6 can  be  obtained  by  differentiating  Eq.  (4.1)  with  respect 

••  •• 

to  the  time  t,  and  using  the  Schrodinger  equation.  The  Schro- 
dinger equation  in  Eq.  (3.3)  can  be  written  as 
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H ( A'  , A'  ) , (4.7) 

~ o 

where  H(A',A^)  is  the  Hamiltonian  in  Eq.  (2.1)  with  the  vector  and 
scalar  potentials  A’  and  A^,  respectively.  The  time  derivative  of 
Eq.  (4.1)  can  be  written  as 

iftd<*'  |6'*'>/dt  = -<4»'  |H(A' ,O)0’i/>'> 

+<K-'  | (i*9/3t  - qA^)0’^'>  , (4.8) 

where  the  scalar  potential  has  been  taken  out  of  the  Hamiltonian. 
From  Eqs.  (2.1)  and  (4.5)  it  can  be  seen  that 

H ( A 1 , 0 ) * (p  - qA/c) 2/2m  + V = H’(A,0)  (4.9) 

is  form  invariant  under  gauge  transformation.  Thus  when  Eqs.  (4.9), 
(4.4),  (4.2)  and  (4.1)  are  used,  it  can  be  seen  that  Eq.  (4.8)  is 

manifestly  gauge  invariant,  i.e.,  it  is  valid  if  the  primes  are 
removed . 

Even  though  Eq.  (4.8)  is  manifestly  gauge  invariant,  it  is 
not  the  mosu  convenient  form  in  which  to  write  the  time  development 
of  the  expectation  value.  The  time  derivative  in  the  last  term 
of  Eq.  (4.8)  can  be  carried  out  to  give 

i-ftd<4/t  | 0 '!)»•  >/dt  « <<!>'  | (6  ' ,H(A'  ,A^)  ) > 

+ <<T  | (if»3e,/3t)K»,> 


(4.10) 
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The  manifest  gauge  invariance  has  been  lost  in  this  traditional 
form,  since  the  right  hand  side  of  Eq.  (4.10)  is  the  sum  of 
two  terms,  neither  of  which  is  gauge  invariant.  However,  since 
the  left  hand  side  of  Eq.  (4.10)  is  gauge  invariant,  the  sum  of 
the  two  terms  on  the  right  hand  side  gives  a gauge  invariant 
result. 

Equation  (4.10)  can  now  be  applied  to  give  the  time  dependence 
of  some  observables,  for  which  the  explicit  gauge  invariance  can 
be  seen  by  inspection.  Since  Eq.  (4.10)  is  gauge  invariant,  the 

primes  can  be  removed.  When  it  is  applied  to  the  position  operator 

41 

r,  the  time  rate  of  change  of  the  average  particle  position  is 

md<ij>  | r<J;>/dt  - <4»|  (p  - qA/c)4<>  . (4.11) 

By  using  Eqs.  (4.1)  and  (4.3)  it  can  be  seen  that  Eq.  (4.11)  is 
manifestly  gauge  invariant,  i.e.,  it  is  valid  with  primes  on  \ , 
r,  and  A.  From  the  classical  analogy  the  operator  on  the  right 
hand  side  of  Eq.  (4.11)  is  used  to  define  the  velocity  operator 

v « m-1(p  - qA/c) , (4.12) 

which  is  form  invariant  under  gauge  transformation  of  the  first 
kind. 


The  time  rate  of  change  of  the  average  velocity  operator 
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v in  Eq.  (4.12)  can  also  be  obtained  by  applying  Eq.  (4.10) 

41 

Without  the  primes) . The  equation  of  motion  is 

<\l>\v\l)>/dt  ■ -<<l>  | ( TV)  + <ij/|qEl(/> 

+ (l/2c)  <\|>|  (V*B  - B*v)i)/>  , (4.13) 

where  the  electric  field  E and  the  magnetic  induction  B are 
defined  in  Eqs.  (3.7)  and  (3.6),  respectively.  Equation  (4.13) 
is  the  quantum  mechanical  analogue  of  Newton's  second  law,  since 
the  right  hand  side  is  the  average  of  the  forces.  The  last  term 
in  Eq.  (4.13)  is  the  average  value  of  the  Hermitian  Lorentz  force 
operator. 

The  energy  of  the  system  should  also  be  a gauge  invariant 
17 

quantity.  The  question  was  previously  raised  in  Eq.(4.6)  as  to  what 
operator  to  use  for  the  energy.  The  Hamiltonian  H(A,Aq)  is  not 
form  invariant  under  gauge  transformation  as  Eq.  (4.5)  shows. 
However,  Eq.  (4.9)  shows  that  H(A,0)  is  form  invariant  under  a 
gauge  transformation  of  the  first  kind.  If  Eq.  (4.10)  is  applied 
to  H (A , 0 ) , the  result  is 

d<i|>  | H (A,  0)  <J*>/dt  = <4»  | Pi)<>  . (4.14) 

The  right  hand  side  of  Eq . (4.14)  is  the  average  of  the  quantum 

mechanical  power  operator  P defined  as 
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P = (q/2)  (E*v  + yjS)  , (4.15) 

where  the  velocity  operator  y is  defined  in  Eq.  (4.12) . Equation 
(4.14)  is  also  manifestly  gauge  invariant,  i.e.,  it  is  valid 
with  primes  on  ^ , A,  and  P.  From  Eq.  (4.13)  it  can  be  seen 
that  the  electric  force  is  qE.  The  electric  force  qE  does 
work,  but  not  the  magnetic  force.  Equation  (4.15)  is  the 
quantum  mechanical  version  of  the  classical  statement  that  the 
power  put  into  the  system  is  the  time  rate  of  change  of  the 
energy  of  the  system.  We  can  therefore  identify  H(£,0)  with 
the  quantum  mechanical  energy  operator. 

V.  GAUGE  INVARIANCE  OF  PROBABILITY  AMPLITUDES 
In  this  section  a gauge  invariant  form  of  time-dependent 

17 

perturbation  theory  is  developed,  following  the  method  of  Yang. 

The  customary  perturbation  theory,  discussed  in  Sec.  II,  is 
based  on  using  Hq  in  Eq.  (2.3)  as  the  unperturbed  Hamiltonian, 
and  the  remainder  of  Eq.  (2.2)  as  the  perturbation.  As  shown 
in  Appendix  A,  the  customary  approach  is  not  gauge  invariant. 

For  an  arbitrary  gauge,  the  interpretation  of  the  expansion 
coefficients  in  the  wave  function  as  probability  amplitudes  for 
finding  the  system  in  an  energy  eigenstate  of  the  unperturbed 

Hamiltonian  H_  is  not  correct. 

o 

19 

As  Lamb  has  noted,  the  correct  gauge  to  choose  to  interpret 
the  expansion  coefficients  as  probability  amplitudes  of  finding 
the  system  in  eignstates  of  the  unperturbed  Hamiltonian  Hq  is 
the  one  for  which  the  Hamiltonian  has  the  electric  dipole  interaction 


given  in  Eq.  (2.5).  This  choice  of  gauge,  however,  is  not  ar- 
bitrary, but  is  seen  to  be  a consequence  of  the  gauge  invariant 
formulation  of  time -dependent  perturbation  theory. 

The  energy  operator  H(A,0),  shown  to  be  form  invariant 
under  gauge  transformations  in  the  last  section,  is  u3ed  as  the 
unperturbed  Hamiltonian.  It  satisfies  the  eigenvalue  problem 

H(A,0)*n(t)  = en(t)*n(t)  , (5.1) 

where  and  are  in  general  functions  of  the  time  t because 
A is  time  dependent.  The  wave  function  for  the  Schrodinger 
equation  in  Eq.  (3.1)  can  be  expanded  in  terms  of  the  complete 
set  of  states  {ty  } as 


= £ cn^n  • (5.2) 

When  Eq.  (5.2)  is  substituted  into  Eq.  (3.1)  and  Eq.  (5.1)  is 
used,  the  equation  of  motion  for  cn  is 

ifl*n  * cncn  +£<lt'nl  (<3Ao"ift9/3t)  Vcm  * (5*3) 

m 

In  contrast  to  Eq.  (A2) , Eq.  (5.3)  is  gauge  invariant.  It  should 
be  emphasized  that  in  Eq.  (5.3)  is  in  general  a function  of 
time  because  of  Eq.  (5.1). 
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It  is  now  shown  that  the  coefficients  cn  in  Eq.  (5.2) 
are  unchanged  under  a gauge  transformation.  If  a gauge  trans- 
formation of  the  first  kind  in  Eq.  (3.2)  is  made  on  the  wave 
function  <(<  in  the  Schrodinger  equation  in  Eq.  (3.1),  and  the 
gauge  transformations  of  Eqs.  (3.4)  and  (3.5)  are  used,  the 
transformed  Schrodinger  equation  in  Eq.  (3.3)  is  obtained.  The 
same  procedure  can  be  used  on  the  eigenvalue  problem  in  Eq.  (5.1) 
to  obtain 

H (A*  ,0)\p‘  = e **  , (5.4) 

— n n n 

where  ip ^ is  related  to  i(/n  by  Eq.  (3.2)  and  where  Eq.  (4.9)  has 
been  used.  The  unperturbed  energy  is  unchanged  under  the 
gauge  transformation  on  the  wavefunction.  Since  both  ' and  tp^ 
are  related  to  \p  and  4/n  by  Eq.  (3.2),  the  expansion  of  in  terms 
of  has  the  same  expansion  coefficients  as  Eq.  (5.2), 

* (5.5) 

Thus,  when  Eq.  (5.5)  is  substituted  into  the  Schrodinger  equation 

in  Eq.  (3.3),  the  equation  for  the  time-dependent  coefficients  c^ 
is 

ihc  - e c I (qA*  - i-ftd/3t)  ip ' >c  ,c 

n n n rn  1 o *m  m . (5.6) 

m 


From  Eqs.  (4.4),  (4.2),  and  (4.1)  it  can  be  seen  that  the  matrix 


elements  in  Eq.  (5.6)  are  equal  to  the  matrix  elements  in  Eq. 

(5.3).  Since  { ) are  independent  of  the  choice  of  gauge  they 
can  be  interpreted  as  the  probability  amplitude  of  finding  the 
system  in  energy  eigenstates  of  H(A,0),  the  physically  meaning- 
ful unperturbed  Hamiltonian. 

Although  Eq.  (5.3)  has  been  shown  to  be  gauqe  invariant, 
its  physical  significance  can  be  clarified  by  rewriting  the 
interaction  term.  If  Eq.  (5.1)  is  differentiated  with  respect 
to  time  it  can  be  shown  (Appendix  B)  that  for  n * m 

<4,J  (qAo~ifl3/3t)'*'m>  = itl(en  ~ em)  P4'm>  <5'7) 

where  P is  the  power  operator  in  Eq.  (4.15).  Therefore  Eq.  (5.3) 
can  be  written  as 

( ifi3/8t  - e ) c = V*  i-fiU  - e)  | P*  > c , (5.8) 

n n / ^ n in  n in  m 

where  the  "dressed"  eigenenergy  e is  defined  as 

"n  = Fn  + (qAo  ” i*a/at)<(/n>  • (5.9) 

Equation  (5.8)  remains  gauge  invariant  because  the  matrix  elements 
of  the  power  operator  P in  Eq.  (4.15)  are  gauge  invariant.  The 
dressed  eigenenergy  « in  Eq.  (5.9)  is  also  gauge  invariant  because 
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both  terms  on  the  right  hand  nj.de  are  gauge  invariant  by  Eqs . 

(5.1) , (4  4),  (4.2)  and  (4.1).  Equation  (5.8)  explicitly  shows 

that  the  transitions  between  states  are  induced  by  the  power 
operator.  In  contrast,  in  the  conventional  time-dependent  per- 
turbation theory  in  Eq.  (A2)  the  transitions  between  states  are 
induced  by  gauge  dependent  potentials. 

Although  Eq . (5.8)  shows  clearly  the  physical  significance 

of  the  interaction  which  induces  transitions,  it  can  be  rewritten 
in  a form  in  which  the  electric  dipole  approximation  becomes 
more  natural.  The  classical  power  can  be  written  as 

Pc ^ = qE*v  d(qE*r)/dt  - q(3E/3t)T 

-qOEj/ax.lXjV.  (classical),  (5.10) 

in  terms  of  the  components  x^  of  r,  etc.  The  summation  convention 
on  repeated  indices  ( i , j = 1,2,3)  is  used.  The  dipole  energy 
qE-r  explicitly  appears  in  Eq.  (5.10).  An  expression  similar  to 
Eq . (5.10)  also  holds  quantum  mechanically.  In  Appendix  C it  is 

shown  that  the  matrix  elements  of  the  power  operator  are 

‘:n:PV  = *nl  IqE-r,  H(A,0)]^m 

-(q/4)<«('n|{{vj,^Ei/8xj),XL}i(<m>  . (5.11) 

The  curly  brackets  denote  the  anticommutator 


1 a , b ) = ab  + ba 


(5.12) 


170 


which  is  used  in  Eq.  (5.11)  to  give  a Hermitian  operator. 

Equation  (5.11)  is  a quantum  mechanical  form  of  Eq.  (5.10). 

If  Eq.  (5.11)  is  substituted  into  Eq.  (5.8),  the  result 

is 

(ifid/3t  - c )c  = |qE*rii/  > 

n n t—j  I n 1 - • m 

m f n 

+ ( ifiq/4)  (tn-tm)"1<^n|{{vj,8Ei/3x.}/xi)^ni>]cm  # (5.13) 

This  equation  is  still  exact  and  is  also  manifestly  gauge  invar- 
iant, since  c , c and  the  matrix  elements  are  explicitly  gauge 
n n 

invariant.  Since  the  electric  field  E * E(r,t)  is  a function 
of  r,  it  cannot  be  taken  outside  of  the  inner  product  in  Eq . (5.17) . 
In  the  next  section  it  will  be  shown  how  Eq.  (5.13)  can  be  sim- 
plified by  making  some  approximations. 
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VI.  ELECTRIC  DIPOLE  APPROXIMATION 

If  the  magnetic  field  in  the  problem  of  interest  is  neg- 
ligible, then  we  can  simplify  Eq.  (5.13)  by  choosing  the  gauge 
such  that  the  vector  potential  is  zero.  If  we  can  further 
assume  that  the  electric  field  is  slowly  varying  over  the  di- 
mensions of  the  atom,  we  obtain  the  electric  dipole  approximation. 

Specifically,  we  choose  the  gauge  function  A in  Eq.  (3.2) 
such  that  the  new  vector  potential  A’  is, 

A*  = A + VA  = 0 . (6.1) 

Thus  the  gauge  function  A is 

r 

A(r,t)  = - /dr1  • Afr^t)  + A(0,t),  (6.2) 

o 

which  is  path  independent.  The  new  scalar  potential  A^  in 
Eq.  (3.5)  can  be  written  as 

r 

A^(r,t)  = Jdr1  • [V1Ao(r],t)  + 3a  (r  , t ) / 3 (ct ) ) 

0 

♦ IAo(0,t)  - 3A(0,t)/d(ct) 1 (6.3) 

If  Eq . (3.7)  is  substituted  into  Eq.  (6.3)  we  have 

r 

Ao(r,t)  * ■ /drl  ' E (Fj_  * t)  (6.4) 

0 

if  A{ 0#t)  is  also  chosen  such  that 
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3A(0,t)/3(ct)  » Ao(0,t) 


(6.5) 


The  Schrodinger  equation  in  Eq.  (3.3)  thus  becomes  in  the 
new  gauge 


H0*'  - q - t)i|/'  = ifi3i)/  '/3t  , (6.6) 

o 

where  Hq  in  Eq . (2.3)  is  the  Hamiltonian  in  the  absence  of  the 

field.  A multipole  expansion  of  Eq.  (6.6)  can  be  made,  but  if 
the  electric  field  is  only  slowly  varying  over  atomic  dimensions 
the  main  contribution  is  the  electric  dipole  interaction, 


[Ho  - q e ( o , t ) • rnp'  = iiiaui'/at  , (6.7) 

for  which  the  Hamiltonian  is  the  same  as  Eq.  (2.5).  Thus  the 
contention  that  the  two  Hamiltonians  in  Eqs.  (2.1)  and  (2.5) 
are  equivalent  if  there  are  no  magnetic  effects  and  the  dipole 
approximation  is  made  is  verified.  The  choice  of  gauge  in  Eqs. 
(6.1)  and  (6.4)  makes  Eq . (5.13)  more  tractable. 

Equation  (5.13)  simplifies  considerably  if  E varies  slowly 
over  the  extent  of  the  wave  functions  i^.  In  this  electric  dipole 
approximation  E can  be  replaced  by  its  value  at  the  origin  E(0,t), 
and  the  spatial  derivatives  of  E can  be  neglected.  Then  Eq . (5.13) 

can  be  written 
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(if»3/3t  ~ e )c 
n n 

= - E (0 , t)  (6'8) 

which  is  still  valid  for  any  choice  of  gauge. 

This  equation  can  be  simplified  still  further  if 
magnetic  effects  are  negligible  so  that  the  gauge  in  which 
A = 0 and  Aq  is  given  by  Eq.  (6.4)  can  be  used.  Then  the 
unperturbed  Hamiltonian  in  Eq.  (5.1)  is  H(0,0)  = Hq  given 
in  Eq.  (2.3),  for  which  the  eigenvalue  problem  is  given  in 
Eq.  (2.10)  with  en  = e^.  The  scalar  potential  Aq  in  Eq.  (6.4) 
becomes  in  the  dipole  approximation 

A *»  -q  E(0,t)  * r . (6.9) 

O ~ 

Thus,  Eq.  (6.8)  can  be  rewritten  as 

( ift  3/3  t - e)  c = V'-q  E(0,t)  • (r ) c , (6.10) 

n n ~ ~ nm  in 

n 

where  tl.c  matrix  elements  (r)  are  defined  in  Eq.  (2.13) 

~ nm 

and  are  now  time  independent.  Equation  (6.10)  is  completely 

equivalent  to  the  Schrodinger  equation  in  Eq.  (6.7). 

Thus  the  Hamiltonian  of  Eq.  (2.5)  is  the  physically 

meaningful  Hamiltonian  to  use  if  H is  chosen  as  the  un- 

o 

perturbed  Hamiltonian.  The  probability  amplitudes  and  en- 
ergies calculated  will  then  be  gauge  invariant,  in  contrast 
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to  the  ones  calculated  from  Eq.  (2.2)  using  the  Coulomb 
gauge  in  Eq.  (2.4).  In  the  next  section  these  ideas  will 
be  applied  to  the  scattering  of  a photon  by  an  electron 
bound  in  an  atom. 

VII.  KRAMERS-HEISENBEFG  FORMULA 

In  this  section  the  equation  for  the  time-dependent 

probability  amplitudes  derived  in  Eq.  (6.10)  is  applied  to 

calculate  the  scattering  cross  section  for  a photon  from 

an  atomic  electron,  which  is  called  the  Kramers-Heisenberg 

(KH)  formula.2®  In  the  conventional  perturbation  theory42 

based  on  the  Hamiltonian  in  Eq . (2.2),  the  three  Feynman 

diagrams  shown  in  Fig.  1 must  be  evaluated  to  obtain  the  KH 

2 

formula  to  first  order  in  the  fine  structure  constant  e /-fie. 
Even  though  the  conventional  perturbation  theory  is  not 
gauge  invariant,  the  KH  formula  is.  A simplification  can 
be  made  by  using  the  gauge  invariant  formulation  of  the  last 
section  to  calculate  the  cross  section.  Then  only  the  two 
Feynman  diagrams  in  Fig.  1(a)  and  (b)  need  to  be  evaluated 
if  the  - ,E-r  interaction  is  used,  where  the  charge  q = -e. 

In  order  to  calculate  the  cross  section  for  the  scat- 
tering of  photons,  the  quantized  electromagnetic  field  is 
4 3 

convenient.  The  vector  potential  for  the  quantized  electro- 

44 


magnetic  field  in  the  Coulomb  gauge  is 
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A (r , t ) = ^c(fi/2u)l/)1/2{akae(a)exp[i(k-r-wt)]  + h.c.}  , (7.1) 

k ,a 

where  h.c.  denotes  the  Hermitian  conjugate  of  the  preceding 

term,  V is  the  volume,  w = kc  is  the  angular  frequency,  and 

A =0.  The  Coulomb  gauge  in  Eq.  (2.4)  is  used,  and  the 
o 

a / \ 

polarization  vectors  e,  (a=l,2)  are  orthogonal  to  the  wave 
vector  k.  The  photon  creation  and  annihilation  operators, 
ak^  and  a^a  , respectively,  satisfy  the  boson  commutation 
relations , 


t . 

ak’a’ 1 

■ 6aa* 

9 

"V. 

V«'' 

= o , 

(7.2) 

, + 

aka' 

'Va'1 

= 0 . 

The  quantized  electric  field  E is  calculated  from  Eq.  (2.6), 
which  gives 

1 /?  a(®) 

E(r,t)  = (ftw/2(/  ) ' ( iak(J  e expl  i (k*r-wt)J  + h.c}  . (7. 3) 

k ,a 

The  electric  field  is,  of  course,  gauge  invariant  and  can  be 
used  in  Eq.  (6.10). 

In  the  initial  state  i the  electron  is  in  the  atomic  eigen- 
state $n  of  Ho  in  Eq.  (2.10)  with  energy  en , while  the  photon 
has  the  wave  vector  k and  polarization  a with  energy  ftto  . in 
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the  final  state  f the  electron  is  in  the  state  # with  energy 

c , while  the  photon  has  the  wave  vector  k ' and  polarization 
m 

a'  with  energy  -fW  . The  intermediate  state  in  Fig.  1 (a)  has 
one  electron  in  the  state  with  energy  r and  no  photons. 

In  contrast,  the  intermediate  state  in  Fig.  1 (b)  has  one 
electron  in  the  state  <f>j  and  two  photons  with  wave  vectors  and 
polarization  (k,a)  and  (k',a').  The  total  energy  of  the  inter- 
mediate state  in  Fig.  1 (b)  is  + tiw  + •fiu)'. 

Equation  (6.10)  can  be  solved  by  iteration  to  second  order 
to  obtain  the  probability  amplitude  c^  (t)  corresponding  to 
Fig.  1 (a) , 


c (2a) 


(t) 


~ (E ' • r ) . (E  • r ) . 

2 4 *->  -o  ~ mg  o - gn 


(ie  /fi)  £ 


, e . - (e  +fio)+i0) 
J D n 


(7.4) 


t 

* J dt,  exp(iAEtj/h)  exp  ( ienT/fi-iemt/fi) 

•T 


The  electric  field  amplitude  E is  defined  as 

~ o 


E = (hu>/2(/)1/2  e 
-o 


(7.5) 


and  E^  is  defined  analogously  with  w and  ^ replaced  by  w ’ and  c' 
respectively.  Equation  (7.3)  for  £ has  been  used  in  the  inter- 

A /V 

action  term  in  Eq.  (6.10),  and  the  polarization  vectors  e and  e' 
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i I 

: I 

I 

correspond  to  polarizations  a and  a'  , respectively.  The  matrix 

elements  in  Eq.  (7.4)  are  defined  in  Eq.  (2.13).  It  is  assumed 

in  Eq.  (7.4)  that  the  interaction  has  been  switched  on  at  time 

-T  in  the  remote  past.  At  the  end  of  the  calculation  the  limit 

as  — T+-«®  will  be  taken.  The  quantity  +i0  in  the  denominator  is 

interpreted  as  an  imaginary  infinitesimal  whose  imaginary  part 

goes  to  zero  through  positive  values,  and  is  required  for  the 

convergence  of  the  integrals.  The  total  energy  change  in  the 

process  of  scattering  from  the  initial  state  with  energy  + -fito 

to  the  final  state  with  energy  + fiw'  is 

m 

( AE  = e + ftw ' - e - lioj  , (7.6) 

m n 

which  turns  out  to  be  zero  when  the  limit  -T-*--°°  is  taken. 

The  probability  amplitude  c^  ' (t)  corresponding  to  Fig.  1 
(b)  is  similar  to  Eq.  (7.4).  The  electric  field  amplitudes 
Eq  and  E^  are  interchanged,  since  the  photon  with  wave  vector  k' 
and  polarization  a'  is  now  emitted  first.  The  energy  in  the 
denominator  of  Eq.  (7.4)  is  replaced  by  e ^ + -fto  + hoo ' , since 
there  are  now  two  photons  in  the  intermediate  state. 

To  calculate  the  transition  probability  per  unit  time,  and 
then  the  cross  section,  it  is  extremely  important  that  the  am- 
plitudes cm  are  true  gauge- invariant  probability  amplitudes.  We 

i 

\ have  established  this  fact  in  Sec.  V.  The  probability  per  unit 


i 

i 
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time  wfi  of  making  a transition  from  the  initial  state  i to  the 
final  state  f is 

wfi  = lim  | c^2a) (t)  + c^2b)(t)|2  /(t  + T) , (7.7) 

in  second  order  of  perturbation  theory,  since  the  numerator  is 
the  probability  of  making  a transition  during  the  time  t + T. 

The  only  part  of  c^2a^  and  c^2b^  that  depends  on  the  time  is  the 
integral  in  Eq.  (7.4),  so  that  the  limit  in  Eq.  (7.7)  gives 

• 2 1 
> lim  | f dt1exp(iAEt1/f»)  | (t  + T) 

T-*-oo  -1 

= 2TTfi  6(AE)  (7.8) 

The  6 function  whose  argument  is  the  change  in  energy  for  the 
scattering  process  in  Eq.  (7.6)  insures  that  energy  is  conserved 
for  the  overall  scattering  process.  In  contrast,  the  energy  of 
the  intermediate  state  j in  Eq.  (7.4)  is  not  equal  to  the  total 
energy  entering 

e . * e + -ft oo  , (7.9) 

J n 

since  the  denominator  would  then  vanish.  The  state  j is  thus 
spoken  of  as  being  a virtual  state,  and  for  these  states  the 
resonance  condition  in  Eq.  (2.15)  is  not  valid.  Thus  in  this 


1 7'J 


case  the  -qE-r  perturbation  is  not  equivalent  to  the  A*p  term 
of  the  conventional  perturbation  theory  in  Eq.  (2.2). 

The  differential  scattering  cross  section  is  defined  as 


da  = ^ wf  j/flux  , (7.10) 

i 

where  flux  = c/V  for  a single  incoming  photon  in  the  volume  V , 
and  the  sum  is  over  all  final  states  in  which  a photon  is  scat- 
tered into  an  element  of  solid  angle  dft  in  the  direction  of  k'. 

If  Eq.  (7.4)  for  c^2a^(t)  and  the  similar  expression  for  cj;21^  (t) 

in  m 

are  substituted  into  Eq.  (7.7),  and  Eq.  (7.8)  is  used,  the  result 
from  Eq.  (7.10)  for  the  cross  section  is 


do/dft  = r2  m2  w (u >')^ft  ^ 
o 


E 

1 


(e ' -r)  . (e  *r)  . 

Vic  - mj  * ~ jn 

oj  . - U) 

jn 


oj . ♦ a) ' 

Dn 


(7.11) 


2 2 

where  r - e /4irmc  is  the  classical  electron  radius  and  w. 

o j n 

is  defined  as  in  Eq.  (2.12).  Equation  (7.11)  is  the  Kramers- 
Heisenberg  (KH)  formula  for  the  scattering  of  a photon  by  an 
atomic  electron.  The  two  terms  within  the  absolute  value  bars 
correspond  to  the  two  diagrams  of  Fig.  1 (a)  and  (b)  with  inter- 
action vertices  eE(0,t)*r.  This  form  with  the  electric  dipole 


180 


matrix  elements  is  the  same  form  as  originally  obtained  by  Kramers 
20 

and  Heisenberg.  A derivation  along  the  lines  given  here  is 

. 45 

given  by  Rower. 

42 

On  the  other  hand,  most  derivations  of  the  KH  formula  use 

the  conventional  perturbation  theory  based  on  the  Hamiltonian 

2 2 2 

in  Eq.  (2.2).  The  (e  /2mc  )A  term  which  is  represented  by  the 
Feynman  diagram  in  Fig.  1 (c) , must  also  be  included  as  a per- 
turbation. The  vertices  in  Fig.  1 (a)  and  (b)  are  given  in  this 

case  by  (e/mc)A*p.  A calculation  similar  to  the  above  gives  the 

42 

cross  section 


do/dn 

- (mtt)'1  2 


(u*/w) 


A At 

6 e^e 
ran  ' ~ 


(£,-F)mj  (£,F)m  + (£-F>m1(<i,-F>in 


“jn“  “ 


01 J + 6)' 


(7.12) 


The  first  term  within  the  absolute  value  bars  is  the  contribution 
2 

from  the  A term  in  Fig.  1 (c) , while  the  next  two  terms  come 
from  Fig.  1 (a)  and  (b) , respectively,  with  the  A*p  vertex.  It 

is  not  obvious  that  Eqs.  (7.11)  and  (7.12)  are  the  same,46  but 

47  48 

Dirac  and  ‘laym  have  shown  that  they  are.  For  the  sake  of 
completeness  the  derivation  is  given  in  Appendix  D using  our 
notation. 

It  is  at  first  surprising  that  the  gauge  dependent  conven- 
tional perturbation  theory  would  give  the  same  result  as  the 
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gauge  independent  electric  dipole  interaction.  As  shown  in 
Appendix  A the  amplitudes  a^  calculated  in  the  conventional 
perturbation  theory  are  gauge  dependent.  Therefore,  the  trans- 
ition probability  calculated  from  Eq.  (7.7)  on  which  the  cross 
section  in  Eq.  (7.10)  is  based  would  be  incorrect.  The  reso- 
lution of  this  problem  is  obtained  in  Appendix  E,  where  it  is 
shown  that  the  difference  between  the  gauge  independent  proba- 
bility amplitudes  cn  and  the  conventional  amplitudes  a^  is 
of  order  e,  the  charge  on  the  electron.  Therefore  the  discrepancy 
between  the  correct  electric  dipole  interaction  and  the  incorrect 
conventional  perturbation  theory  would  show  up  only  in  higher- 
order  perturbation  theory. 


VIII.  CONCLUSION 

In  this  paper  gauge  transformations  in  the  quantum  theory  of 
the  interaction  of  radiation  and  matter  have  been  studied.  The 
special  role  of  the  electric  dipole  interaction  -qE(0,t)*r  in 
problems  in  which  there  is  a negligible  magnetic  field  and  the 
wavelength  of  the  radiation  is  large  compared  to  the  size  of  the 
system  has  been  emphasized.  Calculations  involving  the  electro- 
magnetic field  should  always  be  done  in  a gauge  invariant  way. 

In  order  to  use  an  unperturbed  Hamiltonian  which  does  not  couple 
the  matter  to  the  radiation,  it  is  necessary  to  use  the  electric 
dipole  interaction.  Otherwise,  the  coefficients  in  the  expansion 

of  the  wave  function  in  terms  of  the  unperturbed  states  would  be 
gauge  dependent,  and  cannot  be  interpreted  as  probability  amplitudes. 


The  common  impression  that  the  expansion  coefficients  for 
the  wave  function  can  always  be  interpreted  as  probability  amp- 
litudes is  erroneous,  since  the  probability  that  the  system  be 
in  a given  energy  eigenstate  must  be  gauge  invariant.  Therefore 
the  conventional  perturbation  theory  based  on  an  unperturbed 

Hamiltonian  for  the  matter  uncoupled  to  the  radiation  field 

2 . . 

and  using  the  A-p  and  A interaction  terms  as  perturbations  is 

incorrect,  since  the  expansion  coefficients  are  gauge  dependent. 

The  electric  dipole  interaction  was  illustrated  by  calcu- 
lating the  cross  section  for  the  scattering  of  a photon  by  an 
atomic  electron.  The  cross  section  obtained  in  this  way  is  equi- 
valent to  the  Kramers-Heisenberg  formula  calculated  from  the  con- 
ventional perturbation  theory  using  the  potentials.  The  reason 
for  the  agreement  is  that  only  terms  to  second  order  in  the  fine 
structure  constant  are  included.  In  higher  orders  a discrepancy 
would  appear  between  the  two  approaches,  and  the  correct  result 
would  be  the  one  in  which  the  electric  dipole  interaction  is  used. 

The  gauge  invariant  formulation  of  the  perturbation  theory 
discussed  in  Sec.  V is  general,  and  can  be  extended  to  problems 
in  which  the  magnetic  field  can  not  be  neglected.  A gauge  in- 
variant formulation  in  terms  of  both  electric  and  magnetic  multi- 

17,30 


poles  can  be  given. 
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APPENDIX  A:  Conventional  Perturbation  Theory 

The  customary  procedure  in  solving  the  time -dependent 
Schrodinger  equation  is  to  use  the  Hamiltonian  in  Eq.  (2.2) 
with  Hq , defined  in  Eq.  (2.3),  as  the  unperturbed  Hamil- 
tonian and  the  remainder  in  Eq.  (2.2)  as  the  perturbation. 

The  unperturbed  Hamiltonian  satisfies  the  Schrodinger  equation 
in  Eq.  (2.10).  The  wave  function  ip  in  the  Schrodinger  equation 
in  Eq.  (3.1)  can  be  expanded  in  terms  of  {<J)n}  as 

* - 2>„  ♦„  (Al> 

n 

If  Eq.  (Al)  is  substituted  into  the  Schrodinger  equation  in 
Eq.  (3.1),  the  resulting  equation  for  the  coefficients  is 

ifi^n  = enan  + 2 I (-q/2mc)  (A-p  + p-£) 

m 

+ (q2/2mc2)A2  + qAQ|  4>m>am  . (A2) 

The  transitions  between  states  are  thus  induced  by  the  gauge 
dependent  potentials  in  Eq.  (A2) . 

If  a gauge  transformation  of  the  first  kind  in  Eq.(3.2) 
is  made  on  the  wave  function,  and  Eqs.  (3.4)  and  (3.5)  are 
used,  the  new  Schrodinger  equation  obtained  is  given  in 
Eq.  (3.3).  The  wave  function  \p'  can  be  expanded  in  terms  of 


v 
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the  eigenfunctions  in  Eq.  (2.10)  as 

♦ ' = £ an  4>n  • (A3) 

When  Eq.  (A3)  is  substituted  into  tne  Schrodingei  equation 
in  Eq.  (3.3),  the  equation  for  the  coefficients  a^  obtained 
is 


ifia^  = enan  + ^ <4>n  I I (~q/2mc)  (A'-p  + p-A') 

m 

+ (q2/2mc2)  (A1)2  + qA^  | ^>3^  . (A4) 


The  amplitude  a'  is  in  general  not  the  same  as  a . The 
n n 

question  then  arises  as  to  what  is  the  correct  probability 

amplitude  of  finding  the  system  in  the  state  with  energy 

19 

e^.  This  question  was  answered  by  Lamb  when  he  said  that 

on  the  basis  of  experiment  the  correct  gauge  to  use  is  the 

one  for  which  A = 0 as  given  in  Section  VI. 

That  the  coefficients  a and  a'  are  not  in  general  the 

n n 

same  can  be  seen  by  substituting  Eqs.  (Al)  and  (A3)  into  Eq. 
(3.2)  which  gives 

an  * L I exp  (iqA/flc)  $m>am  . (A5) 

Since  A = A(r,t),  the  coefficients  in  Eq.  (A5)  will  in  general 
be  complex  numbers  not  equal  to  6 


APPENDTX  B:  Matrix  Elements  of  the  Power  Operator 


In  this  appendix  Eq.  (5.7) 
method  of  Yang.17  If  Eq.  (5.1) 
to  time,  we  obtain 


is  derived  following  the 
is  differentiated  with  respect 


H(£,0)i|;  + H(£,0)i  = e*  + ei  . (Bl) 

m **  m mm  mm 


If  the  inner  product  of  Eq.  (Bl)  is  taken  with  respect  to 

ip  , the  result  is 
n 


e 6 + 

m nm 


( e -e  ) <>p  I Ip  > 

m n rn 1 ym 


<«;n|H(A,0)  U/m> 


(B2 ) 


The  time  derivative  of  H(A,0)  in  Eq.  (4.9)  is 


H ( A , 0 ) 


(q/2 ) ( (E  + VAo)-v  + V* (E+VAq) ) , 


(B3) 


when  Eqs . (3.7)  and  (4.12)  are  used.  However,  the  commutator 

of  H (A, 0)  with  q Aq  is 


' H(A,0) ,qAQ) 


(-ifiq/2)  (v*VAq  + VAq*v) 


(B4 ) 


Therefore  Eq.  (B3)  can  be  written  as 

H(A,0)  = P + (ifi)  "1[qAo,H(A,0)  I 


(B5) 


where  P is  the  power  operator  defined  in  Eq.  (4.15). 
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Equation  (B2)  then  becomes 

f-m  Snm  + (ifi)'1(cn-cm)<*nl(qA0-ih3/3t)*,m> 

= ^rJPV  • (B6) 

For  n 4 m Eq.  (5.7)  is  obtained,  while  for  n = m,  we  have 
a special  case  of  Eq.  (4.14). 
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APPENDIX  C:  Power  Operator  and  Electric  Dipole  Energy 

In  this  appendix  Eq.  (5.11)  is  derived.  A straight- 
forward calculation  shows  that  the  commutator  of  r and  H(A,0) 
is 


I r,H(A,0)l  * ihv  , (Cl) 

where  the  velocity  operator  v is  defined  in  Eq.  (4.12).  If 
Eq.  (Cl)  is  substituted  into  the  matrix  element  of  the  power 
operator  in  Eq.  (4.15),  the  result  is 

<4'Jp>J'm>  = (q/2ib)  f (tm-cn)<  «»  lE-r<J>  > 

n m m n n m 

-<4»nl  (E.H(A,0)xi  - xiH(A,0)Ei)ij»m>  1 (C2) 

where  Eq.  (5.1)  has  been  used. 

If  the  commutator  of  E^  and  H(£,0)  is  calculated  the 
result  is 

|£.,H(A,0)]  = (ift/2)  (Vj  OE./Sx.)  + OE^aXjlVj)  (C3) 

for  i * 1,2,3,  where  v^(j  = 1#2,3)  are  the  components  of 
the  velocity  operator  v in  Eq.  (4.12).  When  Eq.  (C3)  is 
substituted  into  the  last  term  in  Eq.  (C2) , the  result 
is 
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<*J  (E1H(A,0)x1  - xiH(f,0)Ei))|/m> 

= (r.  -e  )<<!'_  |E*ri)/  > 

n m n - m 

+ (ift/2)<^nl  ((Vj  f3Ei/3xj}  , <C4) 

where  the  curly  brackets  are  anticommutator  brackets  defined 
in  Eq.  (5.12).  If  Eq.  (C4)  is  substituted  into  Eq.  (C2) , 
we  obtain  Eq.  (5.11). 
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APPENDIX  D:  Equivalence  of  the  Two  Forms  of  the  Kramers-Heisenberg 
Formula 

The  two  forms  of  the  KH  formula  given  in  Eqs.  (7.11) 

and  (7.12)  have  been  derived  using  the  eE-r  perturbation  in 

2 2 2 

Eq.  (2.5)  and  the  (e/mc)A*p  + (e  /2mc  )A  perturbation  in 

* ^ <v 

Eq.  (2.2),  respectively.  Their  equivalence  has  been  shorn 

47  48 

by  Dirac  and  Baym,  but  for  the  sake  of  completeness  the 
derivation  is  given  here  using  our  notation. 

The  commutation  relation  between  the  components  of 
r and  p in  Eq.  (2.8)  implies  that  the  Hamiltonian  HQ  in  Eq. 
(2.3)  and  r have  the  commutator  given  in  Eq.  (2.9).  These 
two  commutation  relations  can  be  used  to  show  that  the  first 
term  on  the  right  hand  side  of  Eq.  (7.12)  is 

»>'1 

- “Wj  (D1) 

If  Eq.  (Dl)  is  substituted  into  Eq.  (7.12)  and  Eq.  (2.9) 
is  used,  the  cross  section  is 

do/dfl  - (r2  m2  w'/uh2)  * 


Using  Eq.  (2.12)  and  ths  energy  conservation  implied  by 
substituting  Eq.  (7.6)  into  Eq.  (7.8)  we  can  prove  that 


' '“jn 


-w>  1) 


Ufa) ' (u , — u>)  ^ +W 

jn 


(D3) 


and 

“mj“jn  (“wmj  " (“jn+‘“')"1^  “ * (D4) 

Equations  (D3)  and  (D4)  can  be  substituted  into  Eq.  (D2) , 
which  gives  Eq.  (7.11),  except  for  an  extra  tens.  That 
extra  term  vanishes. 


since  the  atomic  states  are  complete,  and  the  coordinates 
commute.  Therefore,  we  have  proved  the  equivalence  of  the 
two  forms  of  the  KH  formula  in  Eqs.  (7.11)  and  (7.12). 


0,  (D5) 


I 
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APPENDIX  E:  Relation  Between  Gauge  Dependent  and  Independent 
Amplitudes 

In  this  appendix  it  is  shown  why  the  conventional  per- 
turbation theory  discussed  in  Sec.  II  and  Appendix  A gives 
the  same  KH  cross  section  in  Eq.  (7.12)  as  the  KH  formula 
calculated  in  Eq.  (7.11).  The  equality  is  established  in 
Appendix  D.  It  is  shorn  that  the  difference  between  the 

coefficients  is  proportional  to  the  square  root  of  the  fine 

2 

structure  constant  e /lie,  so  that  the  discrepancy  between 
perturbation  theory  with  the  electric  dipole  interaction 
and  the  conventional  perturbation  theory  will  show  up  only 
in  higher  orders. 

The  gauge  function  A(r,t)  relating  the  electric  dipole 
interaction  in  Eq.  (2.5)  with  the  Hamiltonian  in  Eqs.  (2.1) 
or  (2.2)  is 

A(r ,t)  - -A(0,t)-r  + A(0,t),  (El) 

from  Eq.  (6.2).  The  coefficients  a^  in  Eqs.  (A3)  and  (A5) 
are,  in  this  case,  the  coefficients  cn  of  Eq.  (5.5).  Thus 
the  gauge  invariant  coefficients  cn  are  related  to  the  gauge 
dependent  coefficients  an  of  the  conventional  perturbation 
theory  by 

cn  " «*P(iqA(0,t)/ftc)£<*  |exp(-i(q/ftc)A(0,t) -rU  >a  (E2) 

m 11  v ' ■ B ' 

I 


f 

I 

I 
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from  Eq.  (A5) . If  the  exponential  in  Eq.  (E2)  is 
expanded,  we  obtain 

cn  * an  - i (q/hc)  £ A(0,t) • (E)nmam  + • (E3) 

Thus  the  difference  between  the  gauge  independent  proba- 
bility amplitude  cn  and  the  gauge  dependent  expansion  co- 
efficient aR  is  of  order  q = -e,  the  charge  on  the  electron. 
Therefore,  the  conventional  perturbation  theory  gives  the 
same  result  for  the  KH  formula  in  Eq.  (7.12)  as  the  gauge 

independent  form  in  Eq.  (7.11)  to  second  order  in  the  fine 

2 

structure  constant  o ■ e /ft c. 
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LIST  OF  FIGURES 

Fig.  1.  The  Feynman  diagrams  used  to  calculate  the 

Kramers-Heisenberg  formula.  When  the  Hamiltonian 
in  Eq.  (2.5)  is  used  only  (a)  and  (b)  are  needed. 
When  the  Hamiltonian  of  Eq.  (2.2)  is  used  all 
three  diagrams  are  needed. 


